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Abstract 


Title  of  Dissertation:  Fault  Tolerant  Rerouting  in  Broadband 

Multiclass  Networks 

Andrew  I.  Vakhutinsky,  Doctor  of  Philosophy,  1996 

Dissertation  directed  by:  Professor  Michael  O.  Ball 

Institute  for  Systems  Research  and 
College  of  Business  and  Management 

Modern  broadband  integrated  service  digital  networks  (B-ISDN)  must  han¬ 
dle  multiclass  traffic  with  diverse  quality  of  service  (QOS)  requirements.  The 
main  purpose  of  our  research  is  to  design  call  rerouting  mechanisms  which  pro¬ 
vide  rapid  restoration  of  network  services  in  case  of  link  failures.  We  suggest 
two  approaches:  virtual  circuit  (VC)  and  virtual  path  (VP)  reroutings.  The 
first  approach  is  more  reactive  while  the  latter  is  more  proactive.  The  appli¬ 
cability  conditions  for  the  first  approach  include  the  availability  of  a  layered 
network  structure  similar  to  VC/VP  architecture  which  is  widely  accepted  in 
asynchronous  transfer  mode  (ATM)  networks.  Another  applicability  condition 
is  the  extent  of  network  failure:  VP  level  restoration  is  designed  for  single  link 
failures  -  the  most  common  in  the  telecommunication  networks.  On  the  other 
hand,  in  case  of  less  predictable  multiple  link  failures,  VC-level  rerouting  is  ap¬ 
propriate.  These  two  rerouting  approaches  vary  in  the  amount  of  time  required 
to  carry  them  out.  Though  both  schemes  are  designed  to  work  in  real  time, 
VP-level  rerouting  tends  to  be  faster  and  can  be  performed  in  an  on-line  mode 
using  pre-computed  paths.  VC-level  rerouting  requires  real-time  computation  of 
routes  which  may  result  in  a  noticeable  impact  on  some  services.  On  the  other 


hand,  VP-level  rerouting  requires  a  substantial  amount  of  off-line  computation 
to  design  the  VP  layout  and  the  backup  routes. 

In  this  dissertation  we  propose  a  new  model  and  associated  algorithms  to 
solve  a  VC-rerouting  problem  in  real  time.  This  model  takes  advantage  of  the 
distributed  network  data  and  computational  resources  by  decomposing  the  prob¬ 
lem  at  an  early  stage  and  then  performing  the  computations  in  a  decentralized 
mode. 

In  order  to  solve  the  fault  tolerant  VP  layout  problem,  we  formulate  a  bi¬ 
criteria  optimization  model  reflecting  the  tradeoff  between  throughput  and  cer¬ 
tain  QOS  requirements.  The  model  involves  a  piece-wise  linear  approximation 
to  the  capacity  allocation  rule  for  variable  bit  rate  connections  statistically  mul¬ 
tiplexed  over  a  VP. 

Both  models  are  formulated  as  integer  programs.  The  solution  method  devel¬ 
oped  employ  relaxation  and  aggregation  of  variables,  feasible  solution  heuristics 
and  valid  inequalities.  The  results  of  the  computational  experiments  presented 
indicate  that  the  methods  developed  are  efficient  and  produce  accurate  solutions. 
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Chapter  1 


Introduction 


The  research  described  in  this  dissertation  involves  modeling  a  broadband  multi- 
class  telecommunication  network  which  is  tolerant  to  the  failures  of  its  elements. 
Fault  tolerance  is  the  ability  of  a  telecommunication  network  to  restore  its  op¬ 
erations  with  minimal  impact  on  the  provided  services  after  a  network  element, 
such  as  link  or  node,  fails. 

The  networks  under  consideration  are  broadband  integrated  services  digital 
networks  (B-ISDN)  which  are  designed  to  support  a  wide  variety  of  service  classes 
including  but  not  limited  to  voice,  data  and  video  in  the  digital  form.  These 
classes  are  different  in  their  required  bandwidth  (varying  from  64Kb/s  for  voice 
service  to  several  Mb/s  for  some  video  applications  and  data  transmission),  traffic 
speed  variation  (they  can  be  constant  bit  rate  (CBR)  or  variable  bit  rate  (VBR), 
the  latter  may  have  extremely  high  peak  to  mean  bit  rate  ratio)  and  delay 
tolerance. 

Modern  telecommunication  networks  providing  these  services  represent  a  va¬ 
riety  of  systems.  Among  this  variety  we  specifically  consider  large,  geographi¬ 
cally  dispersed  networks.  These  networks  generally  use  fiber-optic  high  capacity 
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links.  In  fact,  a  physical  link  in  the  modern  telecommunication  network  may 
be  capable  of  transmitting  at  speeds  in  the  Gb/s  range.  Several  types  of  net¬ 
work  architectures  have  been  proposed,  or  are  under  a  development,  to  support 
B-ISDN.  Probably  the  most  prominent  is  Asynchronous  Transfer  Mode  (ATM). 
We  use  some  ATM  specific  assumptions  in  our  work  but  our  models  apply  to 
other  cases  as  well.  The  first  assumption  is  a  connection-oriented  fixed  length 
packet  switching  mode  of  operation  which  means  that  the  information  transmis¬ 
sion  flow  between  the  end  users  is  implemented  in  the  form  of  packets,  or  cells, 
having  the  same  size  and  following  the  same  route,  or  sequence  of  links.  The 
second  assumption  is  the  availability  of  a  layered  network  structure  similar  to 
Virtual  Channel/Virtual  Path  (VC/VP)  architecture  which  is  widely  accepted 
in  ATM  networks.  This  type  of  architecture  allows  for  the  construction  of  what 
can  be  considered  as  additional  logical  links  in  the  network.  Our  modeling  as¬ 
sumptions  concerning  telecommunication  networks  are  described  in  more  detail 
in  chapter  2. 

In  order  to  create  a  testbed  for  the  proposed  call  rerouting  algorithms  and 
compare  performance  measures  of  different  VP  layouts,  we  perform  an  on-line 
call  routing  simulation  which  provides  us  with  a  meaningful  distribution  of  call 
routes  in  a  network.  This  simulation  method  is  described  in  chapter  3.  During 
the  simulation  process,  incoming  calls  of  a  certain  type,  origin  and  destination 
are  generated  according  to  Poisson  distribution.  Call  holding  times  are  generated 
as  exponentially  distributed.  We  consider  calls  of  different  types  (data,  voice  and 
video)  arriving  with  different  rates  between  different  node  pairs.  Each  call  type 
may  have  different  bandwidth  demand  and  mean  holding  time.  This  simulation 
stops  when  the  network  reaches  a  steady  state.  An  embedded  routing  algorithm 
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computes  a  call  path  between  origin  and  destination  taking  into  account  call 
parameters  and  current  network  status.  We  ran  our  simulation  for  different 
parameter  settings  varying  the  steady  state  network  utilization  level  and  call 
type  mixtures.  The  distribution  of  the  calls  in  the  network  is  used  as  input  data 
for  computational  experiments  conducted  to  evaluate  the  quality  of  the  VC- 
level  rerouting  algorithms  proposed  in  chapter  4.  It  is  also  used  in  chapter  6  to 
measure  an  adaptability  of  a  VP  layout  by  calculating  call  blocking  probability 
when  the  offered  traffic  is  increased  as  compared  with  the  projected  level. 

The  layered  architectures  of  ATM  networks  make  it  possible  to  reroute  calls 
on  different  levels.  We  consider  two  alternatives:  virtual  circuit  (VC)  level  and 
VP  level  rerouting.  In  the  first  case  each  call  is  treated  as  a  separate  entity,  in 
the  second  case  we  reroute  entire  VPs  to  pre-established  backup  routes  bypassing 
the  failure.  Generally,  VP  level  rerouting  is  preferred  since  it  can  be  carried  out 
much  faster  than  VC  level  rerouting.  This  can  be  explained  by  the  fact  that 
the  number  of  VPs  is  normally  much  lower  than  the  number  of  VCs  and  once 
a  backup  VP  is  pre-established,  a  required  amount  of  capacity  can  be  allocated 
to  it  in  real  time.  However,  VP  level  rerouting  requires  that  certain  conditions 
are  satisfied.  First,  the  backup  VPs  with  zero  capacities  are  to  be  established 
together  with  the  primaries  before  a  failure  occurs.  Second,  backups  should  not 
be  affected  by  the  same  failure  which  disconnects  the  corresponding  primaries. 
The  latter  fact  makes  the  VP  rerouting  non-applicable  in  case  of  multiple  link 
failures  affecting  large  parts  of  the  network.  Also,  VP  level  rerouting  is  less 
adaptive  and  uses  bandwidth  less  efficiently.  Thus  it  is  less  advantageous  when 
the  call  demand  is  not  steady  or  not  predictable.  Finally,  VP  level  rerouting 
can  be  performed  only  in  those  networks  which  support  the  layered  architecture. 
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On  the  other  hand,  VC  level  rerouting  has  broader  applicability  but  the  delays 
required  to  find  the  VCs  and  set  them  up  will  be  greater  than  the  delays  for  VP 
level  rerouting. 

In  chapter  4  we  present  an  approach  to  real-time  VC  rerouting.  We  assume 
that  there  are  several  priority  classes  of  calls  and  lower  priority  services  can  be 
preempted  to  provide  bandwidth  for  the  higher  priority  traffic.  We  formulate 
the  problem  as  an  integer  program  (IP),  considering  it  as  a  multicommodity  flow 
problem  with  side  constraints  and  zero-one  variables,  and  develop  a  real  time 
approach  to  solve  it.  Our  model  takes  advantage  of  the  distributed  network  data 
and  computational  resources  by  decomposing  the  problem  at  an  early  stage  and 
then  performing  the  computations  in  a  decentralized  mode.  Specifically,  a  two- 
phase  process  is  proposed.  In  the  first  phase,  a  centralized  algorithm  calculates 
a  capacity  allocation  solving  a  linear  relaxation  of  the  problem.  This  allocation 
is  then  distributed  to  each  node.  In  the  second  phase,  each  node  calculates 
paths,  in  parallel,  for  the  calls  it  originated.  In  the  process  of  rerouting  calls, 
our  approach  allows  for  the  possible  preemption  of  existing  calls  based  on  call 
priorities.  Computational  experiments  demonstrate  the  validity  of  our  approach. 

Chapter  5  is  devoted  to  VP  level  rerouting.  There  we  consider  four  different 
scenarios  for  a  reliable,  fault  tolerant  VP  layout.  These  scenarios  vary  in  two 
characteristics:  (1)  whether  the  backup  VP  is  end-to-end  link-disjoint  with  the 
primary  VP  or,  instead,  a  link-specific  backup  VP  is  created  for  each  link  on 
the  primary  VP,  and  (2)  whether  the  backup  has  the  same  or  lower  capacity 
relative  to  the  primary  VP.  For  each  scenario  we  give  an  IP  formulation.  We 
also  consider  the  VP  layout  problem  as  bi-criteria  optimization  problem  which 
involves  trading  off  between  two  performance  parameters: 
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1.  delay  and  control  overhead  minimization 
and 

2.  throughput  maximization. 

In  Chapter  6  we  develop  various  IP  techniques  to  solve  the  problems  formu¬ 
lated  in  Chapter  5.  Since  these  IP  problems  are  quite  large  and  extensive,  it 
appears  impractical  to  attempt  to  solve  them  to  optimality.  Instead,  we  develop 
procedure  to  find  high  quality  feasible  solutions.  Our  approach  is  similar  to  that 
described  in  Chapter  4  and  [6]  and  involves  aggregation  of  variables,  linear  re¬ 
laxation,  valid  inequalities  and  rounding  heuristics.  However,  as  opposed  to  the 
VC  rerouting  problem  applications  which  are  essentially  real-time,  VP  layout  is 
normally  performed  off-line.  Consequently,  when  using  this  approach  in  prac¬ 
tice,  therewould  be  more  computation  time  available.  The  output  to  the  solution 
process  is  a  Pareto-optimal  solution  set,  i.e.  the  set  of  solutions  in  which  neither 
of  the  parameters  can  be  improved  without  worsening  the  other. 

We  performed  a  set  of  computational  experiments  having  a  two-fold  goal: 

1.  to  demonstrate  a  trade-off  between  VP  layout  parameters; 

2.  to  estimate  the  quality  of  our  approximate  IP  solving  methods. 

In  the  first  case,  we  described  approximate  Pareto-optimal  solution  sets  for  model 
networks  by  providing  upper  and  lower  curves  for  the  set  of  exact  optimal  so¬ 
lutions.  In  the  second  case,  we  presented  an  integrality  gap  (i.e.  gap  between 
feasible  and  linear  relaxation  solutions)  and  computation  times. 

Another  performance  measure  which  we  consider  in  our  research  is  adapt¬ 
ability  of  a  particular  VP  layout  to  changes  in  traffic  demand.  Quantitatively, 
we  measure  adaptability  in  terms  of  call  blocking  probability.  Our  conjecture 
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was  that  longer,  lower  capacity  VPs  have  lower  adaptability  than  shorter,  higher 
capacity  VPs.  In  order  to  validate  this  conjecture,  we  simulated  random  genera¬ 
tion  of  calls  which  were  routed  through  the  network  using  previously  computed 
VPs.  The  results  validated  our  conjecture  and  quantified  the  extent  of  this  effect. 
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Chapter  2 


Foundations  of  Broadband  Multiclass 
Networks  and  Modeling  Assumptions 


In  this  chapter  we  give  a  brief  review  of  the  concepts  underlying  Broadband 
Multiclass  networks,  or  using  more  common  terminology  integrated  services  dig¬ 
ital  networks  (B-ISDN).  Recent  monographes  and  reviews  [35,  25,  41]  give  an 
excellent  introduction  into  this  topic. 

The  main  purpose  of  B-ISDN  design  is  to  support  a  wide  range  of  voice  and 
non-voice  applications  in  the  same  network.  These  networks  extend  the  concepts 
of  telephone  networks  by  incorporating  additional  functions  and  features  of  cur¬ 
rent  circuit  and  packet-switching  networks  for  data  to  provide  both  existing  and 
new  services  in  an  integrated  manner.  The  first  set  of  ISDN  recommendations 
was  adopted  by  CCITT  in  1984.  The  primary  rate  access  interfaces  were  defined 
with  total  bit  rates  of  1.544Mb/s  (i.e.  T1  bandwidth)  and  2.048Mb/s  (i.e.  El 
bandwidth)  including  a  64Kb/s  signaling  channel.  However,  it  was  soon  realized 
that  higher  bit  rates  are  required  for  applications  such  as  interconnection  of  local 
area  networks  (LAN),  video,  image,  etc.,  bringing  the  standardization  process 
to  the  introduction  of  of  broadband  concepts.  B-ISDN  is  conceived  as  an  all- 
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purpose  digital  network.  Activities  currently  under  way  are  leading  to  creating 
a  worldwide  network  facilitating  multimedia  information  exchange  between  any 
two  subscribers.  CCITT  Recommendation  1.113  defines  “broadband”  as  “a  ser¬ 
vice  or  system  requiring  transmission  channels  capable  of  supporting  rates  that 
are  greater  than  the  primary  access  rate” .  Currently,  B-ISDN  interfaces  support 
up  to  622  Mb/s  with  the  possibility  of  higher  rates  in  the  future. 

The  network  transport  mode  defines  how  information  supplied  by  the  network 
users  is  eventually  mapped  onto  the  physical  network.  Asynchronous  Transfer 
Mode  (ATM)  is  the  transport  mode  of  choice  for  B-ISDN.  Though  our  research 
network  model  is  not  limited  to  ATM,  we  find  it  useful  to  highlight  the  main 
features  of  this  mode  since  this  is  a  typical  and  probably  the  best  defined  example 
of  the  multiclass  telecommunication  network.  Below  we  highlight  some  of  the 
most  prominent  features  of  ATM  networks. 

•  Connection-oriented  network.  This  means  that  once  a  connection  is 
established,  the  same  route  through  the  network  is  used  for  the  entire 
holding  time.  However,  the  path  used  by  the  connection  depends  on  the 
state  of  the  network  and  usually  has  to  be  determined  before  the  connection 
is  established. 

•  Packet-switched  network.  The  call  traffic  is  transmitted  as  a  stream  of 
fixed-size  packets,  referred  to  as  ATM  cells.  An  ATM  cell  is  53  bytes  long, 
consisting  of  a  48-byte  information  field  and  a  5-bvte  header.  Each  header 
contains  routing  information  applied  to  switching  at  the  cross-connect 
nodes. 

•  Hierarchical  routing.  Two  types  of  routing  concepts  are  used  in  ATM 
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To  Terminals 


Figure  2.1:  Layered  Architecture  of  ATM  Networks  (as  it  appears  in  [9]) 

networks:  virtual  circuits  (VC)  and  virtual  paths  (VP).  Before  transmis¬ 
sion  can  start,  a  VC  has  to  be  established  between  end  nodes.  The  VP 
is  a  more  permanent  structure.  It  can  be  viewed  as  a  logical  link  in  the 
network  in  the  sense  that  it  provides  a  direct  connection  between  a  pair  of 
nodes  using  a  pre-established  route  and  has  its  own  fixed  bandwidth.  VPs 
normally  carry  several  VCs  multiplexed  together  whose  total  bandwidth 
should  not  exceed  the  total  bandwidth  of  the  VP.  Figure  2.1  depicts  this 
hierarchy.  Notice  that  a  VC  can  generally  traverse  more  than  one  VP.  The 
VP  and  VC  being  traversed  by  a,  cell  are  defined  by  fields  in  the  header 
called  virtual  path  identifier  (VPI)  and  virtual  channel  identifier  (VCI), 
respectively.  These  fields  are  used  for  switching  each  cell  between  links 
upon  the  cell’s  arrival  at  an  intermediate  switching  node. 

We  now  consider  how  features  of  ATM  networks  affect  our  model. 
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2.1  Quality  of  Service  and  Statistical  Multiplex¬ 
ing 

The  networks  under  consideration  carry  different  classes  of  calls.  Each  call  by 
its  nature  may  be  constant  bit  rate  (CBR)  or  variable  bit  rate  (VBR).  As  an 
example,  consider  three  main  call  types:  voice,  data  and  video.  Voice  calls  are 
generally  CBR  (unless  sound  compression  is  involved),  may  tolerate  compara¬ 
tively  high  cell  loss  (usually  10-6)  and  require  rather  low  bandwidth  (64 Kbit/s), 
but  delay  more  than  0.1s  would  make  necessary  installation  of  sophisticated  echo 
suppression  equipment  and  it  is  not  desirable.  On  the  other  hand,  data  calls  can 
tolerate  rather  high  (10s  or  more)  delays  but  cell  losses  require  retransmission 
of  the  cells.  Also,  data  calls  are  usually  VBR  with  sometimes  very  high  (106 
or  greater)  peak/mean  rate  ratio;  bandwidth  requirements  can  vary  substan¬ 
tially  for  different  applications.  Video  calls  are,  similar  to  voice  calls,  real-time 
services  which  implies  low  delay  tolerance  and  no  retransmission  of  lost  or  cor¬ 
rupted  cells,  but  unlike  voice  calls  they  are  VBR  and  have  very  high  bandwidth 
requirements  (up  to  8 Mbits/ s).  Each  VBR  call  is  characterized  by  its  mean 
and  peak  transmission  rates.  These  two  parameters  are  negotiated  by  the  call 
admission  mechanism  and  can  not  be  exceeded  for  the  entire  holding  time  of  the 
call. 

ATM  networks  provide  certain  Quality  of  Service  (QOS)  parameters  for  each 
type  of  calls.  We  consider  two  most  important  QOS  parameters: 

1.  Cell  Loss  Ratio  (CLR).  Each  fixed  capacity  VP  normally  carries  several 
VBR  VCs  statistically  multiplexed  together  (i.e.  sharing  the  VP  capacity). 
In  the  event  when  total  VC  transmission  rate  exceeds  VP  capacity  some 
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cells  are  dropped  (lost).  CLR  is  defined  as  a  ratio  of  the  lost  cells  number 
to  the  total  number  of  cells  transmitted.  In  ATM  networks  CLR  is  usually 
a  very  small  parameter  ranging  from  10-12  to  10-6. 

2.  Cell  Transmission  Delay  consists  of  propagation,  and  switching/queuing 
delay.  The  first  is  negligibly  low  unless  satellite  links  are  involved.  The 
second  part  can  be  reduced  by  eliminating  switching  in  intermediate  nodes. 
It  can  be  accomplished  by  low  hop  routing  and/or  using  end-to-end  VPs. 

Let  us  consider  how  the  two  QOS  parameters  are  influenced  by  the  buffer  size 
and  output  link  rate.  CLR  is  approximately  equal  to  the  probability  that  the 
buffer  is  full  upon  arrival  of  a  cell.  Call  transmission  delay  is  basically  the  sum 
of  a  relatively  small  propagation  delay  and  a  queuing  delay  which  is  determined 
by  the  average  length  of  the  buffer  queue  found  upon  the  cell  arrival.  Thus 
increasing  the  buffer  size  would  reduce  the  CLR  but  could  result  in  the  queuing- 
delay  becoming  prohibitively  large.  On  the  other  hand,  increasing  the  output 
link  rate  reduces  both  CLR  and  call  delay  but  may  lead  to  inefficient  use  of  the 
network  resources.  Clearly,  the  output  link  rate  can  not  be  less  than  total  mean 
arrival  rate.  Consider  two  marginal  cases.  If  the  output  link  rate  is  set  exactly 
to  mean  rate  and  buffer  size  is  extremely  large,  then  CLR  is  close  to  zero  but 
transmission  delay  can  be  rather  large,  sometimes  reaching  hours  as  is  the  case 
with  the  Internet.  If  the  output  link  rate  is  set  to  the  total  peak  rate,  then 
even  a  zero  length  buffer  would  provide  zero  CLR  but  the  link  would  be  idle 
most  of  the  time,  since  for  majority  of  the  services,  peak  to  mean  rate  ratio  is 
high.  In  order  to  optimize  ATM  network  design,  satisfying  both  low  CLR  and 
a  transmission  delay  criteria,  and  giving  an  acceptable  utilization  of  network 
resources,  the  choice  is  made  for  sufficiently  small  buffer  size  and  an  output  link 
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rate  which  is  between  total  mean  and  peak  arrival  rates  taking  advantage  of 
statistical  multiplexing  gain.  This  approach  results  in  low  CLR  (less  than  10-6) 
and  provides  real-time  services  with  transmission  delay  of  about  0.1  s. 

Statistical  multiplexing  is  one  of  the  main  advantages  of  ATM  networks  over 
more  traditional  networks  like  circuit-switched  and  synchronous  transfer  mode 
(STM)  networks.  The  result  of  statistical  multiplexing  is  lower  required  band¬ 
width  per  call  without  QOS  degradation  when  several  non-correlated  VBR  calls 
are  multiplexed  together  on  a  link.  Since  the  transmission  rate  distribution  of  a 
general  real  source  is  not  well  understood,  there  is  no  complete  analytical  model 
for  the  statistical  multiplexing.  Several  methods  were  designed  to  determine 
the  amount  of  capacity  to  be  allocated  in  order  to  satisfy  QOS  requirements. 
Among  them  we  mention  equivalent  capacity  approach  by  Gun  et  al  [16,  17]  and 
application  of  the  Shannon  noiseless  coding  theorem  to  its  calculation  in  [44], 

In  chapter  4,  we  develop  a  call  rerouting  model  where  the  total  capacity 
allocated  to  a  group  of  calls  is  computed  by  adding  the  calls’  effective  bandwidths 
([16]).  Later,  in  chapters  5  and  6  we  consider  a  more  complex  model  which 
involves  a  capacity  allocation  function  that  models  statistical  multiplexing  effect 
using  a  piece-wise  linear  approximation. 

2.2  Rerouting  Mechanisms 

There  are  various  rerouting  protocols  in  modern  telecommunication  networks. 
In  this  section  we  outline  the  assumptions  we  make  about  a  rerouting  procedure 
on  both  the  VP  and  VC  levels.  Table  2.1  gives  a  brief  comparison  between 
the  VC-  and  VP-level  in  terms  of  various  issues  such  as  approach  to  rerouting, 
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applicability,  performance  time,  amount  of  spare  capacity  required  and  likelihood 
of  call  refusal. 


Issue 

VC-level 

VP-level 

Type  of  approach 

reactive 

proactive 

Applicability 

wide 

two-layer  architecture 

multiple  failures 

single  link  failures 

Performance  time: 

quasi-real 

real-time 

on-line  computations 

large 

minimal 

off-line  computations 

low 

large 

Spare  capacity 

lower 

higher 

Call  refusal 

possible 

unlikely 

Table  2.1:  Comparison  between  VC-  and  VP-level  Rerouting 


Chapter  5  addresses  the  problem  of  laying  out  primary  and  backup  VPs  such 
that  the  failure  restoration  is  restricted  to  the  VP  level.  It  is  known  [21,  22] 
that  in  ATM  networks  it  is  possible  to  establish  a  zero  bandwidth  VP  between 
a  pair  of  nodes.  This  VP  is  used  as  a  backup  VP  for  a  corresponding  primary 
VP.  We  assume  that  each  end  node  of  the  backup  VP  can  be  a  terminating  node 
of  the  primary  VP  as  well  as  its  intermediate  node.  The  key  assumption  of  the 
VP  rerouting  scheme  is  a  simple  but  effective  restoration  procedure  which  can 
be  performed  with  a  high  degree  of  rapidness  and  reliability  as  soon  as  a  failure 
occurs.  In  this  case  the  failed  primary  VP  is  switched  to  its  backup  as  described 
below. 

Consider  a  link  failure  in  a  network  as  shown  in  Figure  2.2a.  We  assume 
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(a) 


(b) 


Figure  2.2:  Failure  Restoration  Process  (as  it  appears  in  [21]) 
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that  each  network  element  (NE)  maintains  a  local  data  base  for  managing  the 
information  of  spare  resource  utilization  of  the  links  connected  to  the  NE,  and 
that  terminating  NEs  hold  the  pairing  of  occupied  VPs  and  their  backup  VPs. 
If  a  node  receives  a  restoration  message,  it  performs  the  following  procedure  as 
a  transit  node  (Figure  2.2b): 

if  spare  bandwidth  >  required  bandwidth 

then  { 

capture  appropriate  bandwidth  in  the  data  base: 
retransmit  the  message  to  the  next  NE; 

} 

else  { 

return  “Insufficient  Bandwidth”  message 
to  terminate  the  process. 

} 

When  the  upstream  side  node  of  the  restoration  pair  receives  the  restoration 
message,  it  switches  traffic  from  the  failed  VP  to  the  backup  VP;  this  completes 
failure  restoration  for  this  VP  (Figure  2.2c).  This  restoration  algorithm  should 
be  considered  as  very  basic.  In  chapter  5  we  propose  two  enhancements: 

1.  If  there  is  not  sufficient  bandwidth  to  reroute  an  entire  VP,  some  protocols 
may  be  capable  of  taking  down  some  of  the  VCs  they  are  currently  carrying. 

2.  Additional  routing  flexibility  can  be  achieved  by  initially  routing  traffic- 
on  both  primary  and  backup  VPs.  Then,  when  failure  occurs,  traffic  is 
switched  from  one  VP  to  the  other.  In  this  case  there  is  practically  no 
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distinction  between  primal  and  backup  VPs  and  they  can  be  thought  of  as 
“mirror”  VPs. 

In  any  of  the  preceding  variants,  once  the  backup  VP  is  established  with 
appropriate  capacity,  it  is  a  relatively  simple  matter  to  switch  traffic  to  the 
backup  VP.  There  would  be  a  simple  entry  modification  to  the  routing  table  at 
the  beginning  VP  node  so  that  all  cells  whose  outgoing  VPI  was  the  primary 
VP  would  have  the  outgoing  VPI  be  the  backup  VP  after  the  primary  VP  fails. 
The  first  proposed  enhancement  would  require  the  installation  of  a  sort  of  filter 
which  dropped  all  cells  with  particular  VCIs  (or  setting  their  cell  loss  priority 
bit  to  one)  after  switching  to  lower  capacity  backup  VP.  The  second  proposed 
enhancement  can  be  considered  as  part  of  the  basic  algorithm  where  the  new 
outgoing  VPI  of  the  cells  using  failed  VP  would  be  the  same  as  VPI  already  in 
use  by  the  cells  routed  into  the  “mirror”  VP. 

In  lower  speed  private  networks  (our  primary  focus  in  chapter  4)  which  might 
use  ATM  or  fast  packet  switching  (where  VPs  do  not  exist),  we  envision  an  en¬ 
vironment  in  which  backup  bandwidth  will  not  be  reserved  necessitating  the 
preemption  of  lower  priority  calls  in  the  event  of  element  failures.  To  carry  out 
the  appropriate  restoration,  VC/call  level  rerouting  together  with  call  preemp¬ 
tion  will  be  required.  However,  even  in  the  case  of  ATM  networks,  where  full 
use  is  made  of  the  VP  concept,  this  routing  scheme  makes  sense  also  (for  VP 
rerouting),  since  a  single  trunk  failure  can  disrupt  up  to  4096  VPs,  all  of  which 
would  require  rerouting  with  very  high  probability,  and  all  of  which  will  be  dif¬ 
ficult  to  place  on  other  trunks  because  of  the  (in  all  probability)  large  capacity 
of  each.  In  such  cases,  dynamic  rerouting,  rather  than  a  simple  backup  scheme, 
could  produce  substantially  improved  results. 
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Chapter  3 


Testbed  for  the  Evaluation  of  Rerouting 
Algorithms  and  Strategies 


In  this  chapter  we  describe  an  aggregate  level  model  for  generating  multi-media 
calls.  This  model  serves  two  intended  purposes.  The  first  is  to  generate  call  data 
to  be  used  to  compare  routing  algorithms  and  routing  strategies.  The  second  is 
to  generate  a  sample  list  of  calls  affected  by  a  link  failure  which  is  used  as  an 
input  to  VC-level  rerouting  algorithm. 

3.1  Model  Description 

We  consider  a  connection-oriented  telecommunication  network  where  call  re¬ 
quests  arrive  according  to  a  certain  distribution.  The  call  is  admitted  if  there 
is  a  route  between  its  origin  and  destination  nodes  satisfying  its  bandwidth  re¬ 
quirements.  Each  session  (or  call)  uses  the  same  route  for  its  entire  holding- 
time. 

The  simulation  procedure  consists  of  generating  random  call  request  param¬ 
eters  and  applying  a  routing  algorithm  to  establish  a  connection.  A  call  request 
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may  be  rejected  or  blocked  if  the  routing  algorithm  can  not  find  an  available, 
feasible  route. 

Arriving  calls  belong  to  different  classes  (e.g.  voice,  data,  or  video)  and  each 
given  node  pair  may  have  its  own  mixture  of  different  classes.  For  each  call  the 
following  parameters  are  randomly  generated: 

1.  origin  and  destination; 

2.  required  bandwidth; 

3.  arrival  and  holding  times. 

We  assume  there  is  a  set  H  of  different  call  types.  These  types  may  differ  in 
quality  of  service  (QOS)  requirements  including  bandwidth,  maximum  tolerable 
delay,  maximum  cell  loss  probability.  For  the  sake  of  simplicity,  we  assume 
that  both  service  request  interarrival  times  and  holding  times  are  exponentially 
distributed,  though  it  is  not  a  binding  factor  in  our  model.  Different  pairs  of 
nodes  have  different  arrival  rates  for  each  of  the  services,  denoted  by  for  node 
pair  (i,j)  and  call  type  h.  Each  call  type  h  has  mean  holding  time  1  / //.*.,  the 
same  for  each  node  pair. 

Each  network  node  i  is  assumed  to  have  an  associated  parameter  called  node 
population  denoted  by  I).  The  origin  and  destination  of  an  arriving  call  are 
chosen  randomly  based  on  the  following  probabilities: 

Tj 

Prob{  Origin  =  >}  =  - 

A  all  kTk 
T 

Pro6{Destination  =  j  |  Origin  =  i\  =  — - - — — 

^all  k,k^i^k 

This  produces  a  mix  of  calls  such  that  number  of  calls  outgoing/incoming  from/to 
a  particular  node  proportional  to  its  population. 
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For  each  call,  its  origin  and  destination  and  required  bandwidth  are  used  as 
an  input  into  the  routing  algorithm.  It  is  assumed  that  the  routing  algorithm 
has  access  to  network  status  information  such  as  the  load  on  each  link.  Based 
on  that  information  a  decision  on  call  routing  is  made.  If  the  call  is  routed, 
the  network  information  is  updated  by  subtracting  the  call  bandwidth  from  the 
capacity  available  on  links  carrying  the  call.  If  there  is  not  enough  capacity  to 
route  the  call,  the  call  is  blocked.  A  routed  call  stays  in  the  network  for  its 
holding  time.  After  that  it  is  released  from  the  network  and  the  network  status 
is  updated. 

Each  call  is  routed  through  the  network  by  an  on-line  algorithm.  The  on¬ 
line  algorithm  deals  with  the  call  requests  one-bv-one,  accesses  only  current 
network  status  information,  and  can  not  use  any  a  priori  knowledge  of  the  entire 
request  sequence.  The  theory  of  the  on-line  routing  algorithms  is  well  studied 
in  the  literature  (c.f.  [40,  36]).  We  adopt  one  of  the  algorithms  proposed  in 
[13].  Essentially,  this  algorithm  is  a  type  of  shortest  path  algorithm  where  link 
length  is  chosen  by  exponential  scaling  of  its  congestion.  The  algorithm  can  be 
described  briefly  as  follows: 

A  call  between  origin  O  and  destination  D  is  routed  over  O — D  path  P,  if 

1.  P  is  a  min-hop  path; 

2.  ue  +  r /Be  <  1  for  all  e  €  P; 

3.  J2eeP  K'Ue  —  K- 

where  ue  is  a  relative  congestion  of  link  e,  Be  -  its  capacity  and  1  <  k  <  2 
is  a  parameter. 
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The  idea,  behind  this  approach  is  to  introduce  exponential  penalties  for  high  link 
utilization  and  not  to  allow  an  excess  of  the  link  capacity. 

It  was  proved  in  [4]  that  this  algorithm  performance  (in  terms  of  the  through¬ 
put)  is  at  least  l/0(log  n)  of  the  best  off-line  algorithm  performance,  where  n  is 
the  network  size.  Though  this  ratio  does  not  seem  sufficiently  small  from  a  prac¬ 
tical  point  of  view,  simulations  show  that  this  strategy  outperforms  traditional 
minimum-hop  and  reservation-based  strategies. 

We  evaluate  the  network  performance  by  the  call  blocking  ratio.  Thus  the 
goals  mentioned  in  the  introduction  to  this  chapter  are  accomplished  in  the 
following  ways:  In  order  to  compare  different  routing  strategies  and  algorithms, 
the  same  stream  of  calls  is  generated  for  all  simulation  runs  but  on  each  run 
a  different  routing  algorithm  is  used.  The  quality  of  the  routing  algorithm  is 
evaluated  by  the  call  blocking  probability.  In  order  to  obtain  a  sample  set  of 
VCs  affected  by  a  link  failure  for  the  VC-level  call  rerouting  algorithm,  the 
simulation  is  run  until  it  reaches  a  steady  state  and  then  we  extract  a  list  of  calls 
carried  by  a  particular  link  which  is  considered  as  failed. 

3.2  Implementation 

The  simulation  is  event  driven.  That  is,  we  maintain  a  list  of  events  and  keep 
extracting  the  minimum  time  event  from  the  list.  We  also  keep  updating  the 
list.  The  only  events  stored  in  the  list  are  calls  departures,  since  each  arriving 
call  is  either  routed  immediately  or  blocked.  Thus,  there  is  no  queue  of  arrivals 
and  we  have  to  keep  track  only  of  the  next  call  arrival.  If  the  next  event  is  call 
arrival,  we  try  to  route  the  call.  If  the  call  is  blocked,  we  do  not  consider  it 
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anymore;  otherwise,  we  keep  its  route  and  store  its  departure  time  in  the  list  of 
events.  If  the  next  event  is  call  departure,  the  call  is  extracted  from  the  list  and 
its  route  is  removed  from  the  network.  The  simulation  stops  after  a  prespecified 
number  of  calls  is  generated. 

Below,  we  give  more  formal  description  of  the  algorithm  as  a  pseudocode. 

Inputs  to  the  Algorithm 

1.  simulation  length 

2.  network  information 

(a)  node  list 

for  each  node: 

i.  node  population 

ii.  call  type  mix 

(b)  link  list 

for  each  link: 

i.  end  nodes 

ii.  link  capacity 

3.  call  type  list 

for  each  call  type: 

(a)  required  bandwidth 

(b)  mean  holding  times 
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Outputs  from  the  Algorithm 

1.  call  list 
for  each  call: 

(a)  nocle-or  (origin  node)  and  nodejdest  (destination  node) 

(b)  call  type 

(c)  traff  (mean  traffic  rate  for  call) 

(d)  VC  -  given  by  linki,  link 2, . . . ,  linkn 

The  Algorithm 

1.  CallArrival  =  0: 

2.  for  Counter  =  0  to  SimLen 

3.  if  (EventList  is  empty  OR  CallArrival<MinTime(EventList)) 

4.  GenerateCallQ 

5.  if  (FindRoute(CallParameters)  ==  SUCCESS) 

6.  CallDeparture  =  CallArrival  +  HoldTime(^); 

7.  AddCallToList(EventList); 

8.  UpdateLinklnf o(Network,  CallParameters); 

9.  else  block  the  call; 

10.  CallArrival  =  CallArrival  +  Interarrival (A) ; 

11.  else 

12.  CallRelease (Network ,  CallParameters); 

13.  RemoveCallFromList (EventList); 
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VARIABLES: 

Counter 

event  counter; 

SimLen 

simulation  length  (number  of  events); 

CallArrival 

arrival  time  of  the  next  call; 

CallDeparture 

departure  time  of  the  next  call; 

CallParameters 

call  parameters  (type,  bandwidth); 

Network 

network  data; 

EventList 

list  of  currently  active  calls. 

FUNCTIONS: 

InterarrivalQ 

generates  interarrival  time; 

GenerateCall() 

generates  call  request; 

HoldTime() 

generates  holding  time  of  the  call; 

FindRoute() 

finds  a  route  for  the  call; 

AddCallToList() 

adds  call  to  the  list  of  currently  active  calls; 

RemoveCallFromLi st  ( ) 

removes  call  from  the  event  list; 

MinTime() 

finds  the  minimum  departure  time; 

CallReleaseQ 

releases  the  bandwidth  capacity  on  the  links 

UpdateLinklnf o() 

adds  call  to  each  carrying  link. 
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Call  arrival  times  are  computed  by  adding  a  randomly  generated  interarrival 
time  to  the  last  arrival  time  (line  10).  The  first  arrival  is  set  to  0  (line  1).  Inter¬ 
arrival  times  are  generated  as  exponentially  distributed  random  variates  using 
the  inverse-transform  method  within  function  Interarrival(A).  The  departure 
time  of  a  call  is  computed  as  its  arrival  time  plus  a  randomly  generated  holding 
time.  The  holding  time  is  generated  as  an  exponentially  distributed  random 
variate  using  the  inverse-transform  method  within  function  HoldTime(^). 

In  the  beginning  of  any  simulation  loop  (line  2),  we  determine  the  next  event 
by  comparing  the  next  arrival  time  and  the  minimum  departure  time  from  the 
event  list  (line  3).  If  the  next  event  is  arrival,  the  call  request  is  generated  (line 
4)  and  the  routing  function  is  called  to  compute  its  route  (line  5).  If  the  call 
gets  through  (i.e.  a  feasible  route  is  found),  the  departure  time  of  the  call  is 
computed  (line  6),  the  call  departure  time  is  added  to  the  list  of  events  (line 
7),  and  its  parameters  are  added  to  the  lists  of  calls  routed  through  the  links 
carrying  the  call  (line  8).  If  the  call  is  blocked,  it  is  disregarded  by  the  algorithm. 

If  the  next  event  turns  out  to  be  call  departure,  then  the  call  is  taken  down 
(line  11).  That  means  releasing  capacity  bandwidth  on  the  links  used  by  the  call 
(line  12)  and  removing  the  call  from  the  event  list  (line  13). 

The  simulation  ends  when  the  number  of  events  defined  by  the  event  counter 
exceeds  the  simulation  length. 


24 


Chapter  4 


Distributed  VC-level  Call  Rerouting  in 
Multiclass  Broadband  Networks 


In  this  chapter  we  present  mechanisms  for  rerouting  around  faults  in  multiclass 
traffic  environments  which  will  be  encountered  in  most  of  the  future  ubiquitous 
broadband  integrated  service  digital  networks  (B-ISDN).  These  networks  should 
deliver  service  for  calls  with  widely  varying  bandwidth  and  quality  of  service 
(QOS)  requirements.  In  order  to  provide  high  service  availability,  it  is  necessary 
for  the  networks  to  react  in  real-time  to  network  faults,  the  topic  of  our  research. 

The  rerouting  problems  treated  in  this  chapter  involve  the  simultaneous  cal¬ 
culation  of  routes  for  several  calls.  Our  approach  differs  from  classical  which 
is  routing  in  circuit-switched  networks  [23]  in  that  we  are  most  concerned  with 
satisfying  a  large  group  of  calls  and  the  variety  of  bandwidth  requirements  giving 
the  problem  a  “packing”  dimension.  That  is,  calls  with  varying  bandwidth  must 
be  packed  into  bandwidth  limited  links.  The  “bandwidth  packing  problem”  has 
been  studied  in  the  literature  [24,  37,  38].  Previous  research  has  concentrated  on 
finding  an  optimum  or  near-optimum  solution  by  applying  off-line  algorithms. 
In  this  chapter  we  develop  a  model  and  associated  algorithms  designed  to  oper- 
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ate  in  a  distributed  environment.  Specifically,  a  two-phase  process  is  described. 
In  the  first  phase,  a  centralized  algorithm  calculates  a  capacity  allocation.  This 
allocation  is  then  distributed  to  each  node.  In  the  second  phase,  each  node  cal¬ 
culates  paths,  in  parallel,  for  the  calls  it  originated.  In  the  process  of  rerouting 
calls,  our  approach  allows  for  the  possible  preemption  of  existing  calls  based  on 
call  priorities. 

Routing  standards  and  schemes  typically  address  large-network  issues  by 
network  aggregation  techniques.  The  two-pliase  routing  scheme  proposed  in  this 
chapter  (which  does  not  explicitly  consider  network  aggregation)  has  immediate 
applicability  within  subnetworks.  These  “subnetworks”  themselves  may  be  quite 
large-for  example,  one  could  view  a  corporate  network  (such  as  IBM’s)  as  a 
subnetwork.  In  such  cases,  the  “local  graph”  of  such  a  scheme  as  the  Private 
Network-to-Network  Interface  (PNNI)  is  large  enough  to  present  the  challenging 
problems  that  we  are  trying  to  solve  here.  For  interworking  between  subnetworks 
(i.e.  where  other  subnetworks  are  seen  only  in  an  aggregated  sense)  this  scheme 
would  have  to  be  adapted  in  two  ways:  First,  it  would  have  to  be  adapted  to 
multilayered  rerouting,  since  as  it  stands  in  this  chapter  it  addresses  only  a  single 
network  layer.  Second,  we  would  have  to  find  an  appropriate  way  of  representing 
the  allocations  in  an  aggregate  form  for  higher  layer  clusters.  While  we  have  not 
yet  done  this  work,  we  expect  that  multilayer  schemes,  like  the  PNNI,  can  be 
accommodated  with  reasonable  effort. 

Another  unique  characteristic  of  the  problem  we  consider  is  allowing  the 
possibility  of  preempting  existing  calls.  Preemption  is  used  in  private  networks 
in  order  to  ensure  that  service  to  the  highest  paying  customers  is  disrupted 
least  and  established  or  restored  most  quickly.  For  virtual  path  networks  (VPN) 
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by  public  carriers,  preemption  can  be  an  effective  mechanism  for  sharing  links 
between  casual  customers  and  customers  to  whom  service  is  guaranteed  as  part 
of  a  VPN  contract.  Such  a  priority  mechanism  can  ensure  that  the  capacity 
guaranteed  to  a  VPN  customer  is  always  available. 

The  chapter  is  structured  as  follows:  section  4.1  gives  a  model  description  and 
general  approach.  The  algorithms  are  described  in  section  4.2,  computational 
experiments  are  presented  in  section  4.3  and  section  4.4  summarizes  the  chapter 
contributions. 


4.1  Problem  Formulation 

4.1.1  Assumptions  and  Fault  Scenario 

The  basic  scenario  we  are  presented  with  is  the  failure  of  a  network  link,  e, 
in  an  environment  in  which  we  may  assume  all  network  elements,  including  e, 
are  currently  handling  significant  amounts  of  traffic.  In  general,  at  the  time 
of  its  failure,  several  calls  will  be  routed  over  e.  We  denote  the  set  of  these 
calls  by  1C.  The  problem  to  be  addressed  is  the  determination  of  new  routes 
for  all  calls  in  1C.  Due  to  the  dynamics  of  the  network  information  used  by  the 
originating  node  when  setting  up  the  call,  it  is  possible  that  when  the  call  setup 
is  actually  carried  out,  it  will  be  refused  by  an  intermediate  node  due  to  the 
unavailability  of  sufficient  resources  on  a  link.  We  call  this  phenomenon  a  call 
refusal.  Alternatively,  it  is  possible  that  in  order  to  set  up  the  incoming  call  it 
is  necessary  to  take  down  the  route  of  an  existing  lower  priority  call.  We  call 
this  phenomenon  call  preemption.  In  such  cases,  the  originating  node  of  the 
preempted  call  would  try  to  find  a  new  route  for  that  call.  We  note  that  each 
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intermediate  node  along  a  call’s  path  contains  a  limited  amount  of  information 
concerning  the  call,  including  a  call  id,  the  originating  node,  call  bandwidth  and 
outgoing  link.  Certain  rerouting  strategies  might  require  that  this  information 
be  augmented. 

We  consider  two  failure  scenarios  for  a  network  link,  e.  Under  the  first 
scenario,  the  network  management  system  anticipates  the  failure  of  e  and  the 
system  has  sufficient  time  to  plan  route  adjustments  accordingly.  Under  the 
second  scenario,  link  e  experiences  a  sudden  failure.  In  the  first  case,  we  could 
assume  that  the  network  management  system  has  a  “reasonable  amount”  of 
time  to  react  to  the  failure  and  in  the  second,  the  network  must  react  much 
more  quickly.  In  either  case,  we  assume  that  several  calls  are  routed  over  e. 

The  call  model  is  as  follows:  Each  (single  connection)  call  is  composed  of  two 
paths:  A  forward  path  (relative  to  the  call  originator)  and  a  return  path.  When 
the  call  is  initially  set  up,  both  the  paths  are  established  by  the  originator, 
though  the  two  paths  may  indeed  be  different.  The  two  paths  then  may  be 
rerouted  independently  of  one  another  and  may  be  preempted  independently  of 
one  another.  If  one  of  the  paths,  however,  fails  to  be  rerouted  when  necessary, 
the  call  fails  (unless,  of  course,  the  communication  is  only  in  one  direction).  We 
think  that  this  gives  the  network  great  flexibility  because  it  allows  the  network 
to  find  feasible  routes  for  connections  that  might  otherwise  not  be  found,  thereby 
enhancing  perceived  network  availability  for  customers.  In  addition,  in  case  of 
link  failure,  communication  is  disrupted  only  in  one  direction.  Finally,  in  case 
of  dynamic  bandwidth  adjustment,  failure  to  find  the  necessary  bandwidth  for 
a  path  in  one  direction  does  not  necessitate  the  rerouting  of  the  entire  call  (i.e. 
both  directions).  To  simplify  the  formulation  and  analysis  given,  we  do  not 
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explicitly  address  the  correlation  of  the  forward  and  reverse  paths  in  our  model, 
i.e.  we  proceed  as  if  each  call  requires  a  path  in  one  direction  only.  Our  overall 
scheme  can  address  this  correlation  with  some  modest  modifications  which  we 
do  not  include  here. 

We  assume  routing  computations  are  distributed  throughout  the  network 
so  that  the  node  at  which  the  call  originates  calculates  a  route  for  the  call. 
Specifically,  under  normal  conditions,  the  call-origination  node  selects  a  route 
for  each  call  as  each  individual  call  request  is  received  by  the  node.  The  route 
calculations  are  based  on  link  load  information  passed  to  the  node  periodically 
by  other  nodes.  We  associate  a  fixed  bandwidth  requirement  with  each  call. 
Thus  our  model  most  directly  applies  to  circuit  switched  networks  where  each 
call  may  have  an  arbitrary  bandwidth  requirement.  However,  under  appropriate 
conditions,  it  is  also  applicable  to  ATM  networks  and  fast  packet  switching 
networks  with  variable  length  packets.  We  make  a  key  assumption  with  respect 
to  additivity  of  bandwidth;  that  is,  the  total  bandwidth  used  by  the  calls  on 
each  link  is  sum  of  all  bandwidths  associated  with  calls  routed  over  the  link. 
In  order  to  apply  the  model  to  ATM  and  fast  packet  switching  networks  at 
the  call/virtual  circuit  (VC)  level,  one  would  use  the  call’s  effective  bandwidth 
([16])  as  its  bandwidth  requirement.  The  additivity  assumption  would  imply 
that  effective  bandwidths  over  a  single  link  can  be  added  to  test  the  capacity 
constraint.  There  is  currently  a  large  amount  of  active  research  (e.g.  [17])  related 
to  this  issue  and  the  potential  error  that  might  be  introduced  by  doing  so.  We 
will  not  discuss  details  here  but  only  point  out  that  the  applicability  of  our 
model  in  this  context  clearly  is  impacted  by  the  results  of  this  research.  It  is 
also  possible  to  apply  our  model  to  the  determination  of  routes  for  virtual  paths 


29 


(VPs)  in  ATM  networks.  Since  fixed  bandwidths  are  associated  with  virtual 
paths,  the  additivity  assumption  appears  to  be  valid. 

Our  approach  makes  use  of  a  certain  capacity  allocation  scheme:  Each  call 
origination  node  is  given  its  own  capacity  allocation  for  each  link  in  the  network. 
The  capacity  allocated  to  a  particular  call  origination  node  will  typically  be 
substantially  less  then  the  total  capacity  available  on  the  link.  The  path  selection 
algorithms  that  the  call  origination  node  uses  will  penalize  any  violation  of  the 
capacity  allocation.  When  a  link  fails,  a  designated  node  will  compute  the 
capacity  allocation  for  all  call  origination  nodes  that  have  routed  calls  over  the 
failed  link.  This  capacity  allocation  will  be  based  on  an  aggregate  rerouting 
of  all  affected  traffic.  The  capacity  allocation  will  then  be  transmitted  to  each 
origination  node  which  will  use  the  allocation  when  determining  new  routes  for 
all  calls  previously  routed  over  the  failed  link.  The  motivation  for  this  approach 
is  to  reduce  the  incidence  of  call  preemption  and  call  refusal  that  would  normally 
take  place  when  failures  occur.  Specifically,  it  is  anticipated  that,  when  failures 
occur,  since  the  origination  nodes  are  trying  to  carry  out  “bulk”  call  setup,  the 
information  they  would  have  on  the  traffic  load  on  links  would  be  inaccurate. 
Thus,  it  would  be  much  more  likely  that  call  setups  would  be  attempted  on  links 
that  did  not  have  sufficient  capacity.  The  intent  is  to  reduce  the  occurrence  of 
this  phenomenon  by  effectively  allocating  the  link  capacity  ahead  of  time. 

In  this  chapter,  we  describe  the  underlying  optimization  problems,  propose 
basic  solution  approaches  for  each  and  evaluate  the  merits  of  the  two-phase 
approach  in  determining  high  quality  routes.  We  understand  that  several  addi¬ 
tional  issues  must  be  addressed  in  order  to  make  this  approach  practical.  These 
include: 
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•  determination  of  what  additional  information  must  be  stored  by  each  node; 

•  thorough  analysis  of  capacity  allocation  step,  specifically  -  evaluation  of 
alternative  locations  for  computing  the  allocation,  protocols  for  distribut¬ 
ing  the  allocation  information  and  determination  of  time  required  to  carry 
out  distribution; 

•  consideration  of  pre-calculation  and  storage-for-fast-retrieval  of  the  band¬ 
width  allocation; 

•  the  development  of  a  hybrid  approach  which  reroutes  as  many  calls  as 
possible  by  switching  to  backup  VP’s,  and  then  executes  the  algorithms 
presented  here  to  carry  out  VC  level  rerouting. 

In  sections  4.1.3  and  4.1.2  we  present  formulations  of  the  rerouting  prob¬ 
lems.  The  purpose  of  these  is  to  describe  the  fundamental  problem  structure. 
These  formulations  can  not  be  used  directly,  since  for  the  target  environment 
distributed,  real-time  algorithms  are  required.  In  Section  4.2,  when  specific  algo¬ 
rithms  are  presented,  the  telecommunications  environment  is  taken  into  account. 

4.1.2  Notation 

Let  Q  =  (A/”,  C)  be  the  undirected  graph  representing  a  telecommunication  net¬ 
work;  A f  is  the  set  of  nodes,  C  is  the  set  of  links.  Since  each  link  carries  different 
calls  in  each  of  its  two  directions,  we  consider  each  link  [i.  j]  G  C  consisting  of  two 
directed  arcs:  (?',  j)  and  (j,  i).  The  set  of  all  arcs  is  denoted  by  A,  A  C  Af  x  Af. 
Residual  (free)  capacity  of  the  arc  is  defined  as  the  difference  between  the  total 
link  capacity  and  the  total  bandwidth  of  the  calls  carried  by  the  arc.  In  gen¬ 
eral,  the  arcs  (t.  j)  and  (j,  t)  have  different  residual  capacities  Bj.j  and  Bji  since 
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different  calls  are  routed  over  these  arcs.  We  will  also  introduce  the  following 
notation  for  the  network:  IN(i)  =  set  of  arcs  directed  into  node  i;  OUT(i)  = 
set  of  arcs  directed  out  of  node  i.  Everywhere  further  the  network  is  considered 
consisting  of  directed  arcs. 

As  indicated  previously,  we  let  1C  denote  the  set  of  calls  to  be  rerouted.  For 
each  call  A:  G  /C,  the  following  set  of  parameters  is  defined:  origin  0 ( k )  G  A/", 
destination  D(k )  G  J\f,  bandwidth  bk  and  priority  class  h(k)  G  H.  Additionally, 
we  associate  a  revenue  c,h(k)bk  with  call  k.  Thus  q  can  be  considered  as  a  rev¬ 
enue  obtained  from  routing  each  unit  bandwidth  of  class  l  (We  note  that  the 
requirement  that  revenue  is  a  linear  function  of  bandwidth  is  not  essential  to 
the  overall  approach).  Notice,  that  the  revenue  is  received  only  after  routing  the 
“entire”  call. 

For  the  purposes  of  the  exact  problem  formulation,  we  assume  the  set  of  all 
active  calls  M.  is  known  explicitly.  This  assumption  will  be  replaced  with  a 
reasonable,  practical  approximation  in  section  4.2.  Let 

{1  if  call  m  is  carried  by  arc  e 

(4.1) 

0  otherwise 

and  let  Ch(m)bm  be  the  cost  of  preempting  active  call  m  with  priority  class  h(m) 
and  bandwidth  bm.  Similarly  to  routing  revenue,  q  can  be  considered  as  a  cost 
of  preempting  each  unit  of  class  i  bandwidth. 

A  path  through  the  network  for  the  call  A:  is  defined  as  a  sequence  of  arcs 
(*o,  «i),  (hd-2),  •  •  • ,  (V-i Am)  where  i0  =  O(k)  and  im  =  D(k).  We  will  denote 
the  path  by  Pk-  The  set  of  all  feasible  paths  for  call  A:  G  1C  will  be  denoted  as 
V(k).  In  order  to  “penalize”  long  call  routes,  we  introduce  heuristic  coefficients 
ske  representing  the  penalty  we  pay  by  rerouting  call  k  over  arc  e. 
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4.1.3  Problem  Formulation 


In  this  section,  we  present  an  integer  programming  (IP)  formulation  of  the  prob¬ 
lem.  This  formulation  is  probably  the  most  natural.  It  includes  embedded 
“flow”  problems  for  each  call  to  be  routed.  The  formulation  uses  the  following 
0/1  variables: 

{1  if  call  k  is  rerouted 
0  otherwise 

1  if  call  k  is  rerouted  over  the  arc  e 

for  all  k  E  1C  and  e  G  A  :  xke  =  < 

0  otherwise 

{1  if  call  m  is  preempted 
0  otherwise 

The  exact  IP  formulation  is  given  by: 

(IPR): 


Maximize 

'y  '  Ch(k)bkUk 

y  '  Ch(m)bm,Zm 

Y  Y  44 

(4.2) 

keic 

m.£  A4 

keic  eeA 

subject  to 

f 

-yk  if  D(k)  = 

i 

E  < 

-  E  xe  =  < 

II 

S' 

O 

t+H 

i  for  all  i  G  M  and  k  G  A(4.3) 

c  or /•(•<: 

e£lN(i) 

0  otherwise 

E  bkXke  <  Be 

T  y  '  ni^iii  ■-■ni 

for  all  e  G  A, 

(4,4) 

keK 

m£M 

0<yk<i 

and  integer 

for  all  k  G  JCf 

(4,5) 

0<zm<l 

and  integer 

for  all  rn  G  At . 

(4,6) 

xk  >  0  and  integer  for  all  e  G  A  and  k  G  K[. 4.7) 


Objective  function  (4.2)  maximizes  the  revenue  received  by  routing  the  calls 
(the  first  term),  minimizes  the  capacity  violation  (the  second  term)  and  a  cost 
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which  encourages  shorter  routes  to  be  chosen  (the  last  term).  Coefficients  s\  are 
chosen  to  be  relatively  small  in  comparison  to  cost/revenue  coefficients.  Con¬ 
straints  (4.3)  are  imposed  upon  each  call  A:  are  the  usual  network  flow  conserva¬ 
tion  equations.  Since  variables  x\  are  integral  and  the  flow  value  for  each  call  is 
0  or  1  (the  y  variable  restriction),  the  call  can  be  rerouted  over  at  most  one,  non¬ 
split  path.  Hence  xke  G  {0, 1}.  Constraints  (4.4)  enforce  the  capacity  restriction 
on  each  arc  where  the  upper  bound  on  the  capacity  is  a  “soft”  constraint  -  it 
may  be  violated  by  preempting  the  calls. 

We  now  present  an  alternate  formulation,  which  uses  a  0/1  variable  for  each 
feasible  path.  Thus,  conceptually  it  potentially  requires  a  very  large  number  of 
variables/data.  In  general,  this  should  not  deter  one  from  considering  it  since 
formulations  of  this  type  can  be  quite  practical  if  one  generates  variables  and 
their  associated  columns  dynamically.  Specifically,  for  problems  of  this  type  “col¬ 
umn  generation”  approaches  are  used,  which  generate  columns  of  the  constraint 
matrix  only  as  they  are  required.  It  is  typical  that  the  number  of  columns 
generated  is  of  a  very  manageable  size  and  is  much  less  then  the  total  number 
which  appear  in  the  problem  statement.  Although,  initially  we  do  not  intend  to 
use  formal  column  generation  approaches,  our  approximate  algorithms  will  be 
similar  in  spirit  to  column  generation  and  will  use  many  of  the  same  concepts. 

The  new  variables  for  this  formulation  are  path  variables  defined  as  follows: 


1  if  call  A:  is  rerouted  over  the  path  p 

for  all  k  E  JC  and  p  G  V{k) 


^  0  otherwise 

Flow-Path  Formulation  with  Aggregate  Call  Preemption  Costs: 
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(FPAB): 


Maximize:  JZ  ^h(k)bkyk  - 

^  '-in 

-E 

E  • 

4  fkp 

(4.8) 

k(z)Cr 

mEM. 

ketCpev(k) 

subject  to 

E 

fkp  —  Uk 

for 

all  k 

G  /C, 

(4.9) 

pev(k) 

E  bk  E  fkp 

A  Be  + 

E 

bm  An 

for 

all  e 

G  A 

(4.10) 

k£lC  p:eEpEV(k) 

m£A4 

VI 

o 

yk  <  1 

for 

all  A: 

G  /C, 

(4.11) 

fkp  >  0  and  integer  for  all  k  G  JC  and  p  G  V(k)  (4.12) 
0  <  z,m  <  1  and  integer  for  all  m  G  A4,(4.13) 

Constraints  (4.9)  state  there  is  exactly  one  path  for  the  call  if  it  is  rerouted 
and  no  paths  if  it  is  not.  Constraints  (4.10)  are  similar  to  those  of  (4.4). 

The  path  formulation  and  associated  column  generation  approaches  provide 
a  general  framework  for  the  development  of  a  class  of  approximate  and  exact 
algorithms.  Specifically,  the  process  can  be  viewed  as  starting  with  a  set  of  calls 
to  be  rerouted.  A  list  of  paths  is  built  up  dynamically  over  time.  At  any  given 
time  there  exists: 

•  a  list  of  calls 

•  a  list  of  paths 

•  a  “matrix”  of  patli/call  compatibilities  and  assignment  costs  (e.g.  spk ) 

•  a  (possibly  partial)  assignment  of  calls  to  paths 

An  any  given  iteration  the  algorithm  could  improve  the  assignment  of  calls  to 
paths  or  alternatively,  generate  one  or  more  new  paths.  Of  course,  these  two  op¬ 
tions  are  clearly  linked.  If  it  did  not  seem  possible  to  generate  a  good  assignment 
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of  calls  to  paths  using  the  current  set  of  paths,  then  this  would  signal  to  need  to 
generate  additional  paths.  Furthermore,  paths  could  be  generated  specifically  to 
try  to  accommodate  certain  unassigned  or  poorly  assigned  calls. 


4.2  Algorithms 

In  the  previous  section,  we  presented  exact  IP  formulations  of  the  problem. 
Theoretically,  we  could  have  plugged  those  formulations  into  an  IP  solver  and 
obtained  an  optimal  solution.  Such  an  approach  would  not  be  appropriate  in 
a  real-time  distributed  setting.  We  propose  a  two-phase  approach  in  which  the 
final  decision  on,  and  calculation  of,  a  call’s  path  remains  with  the  call  origination 
node. 

In  this  section  we  outline  the  design  of  the  algorithms  and  overall  problem 
structure  and  derive  some  theoretical  running  time  bounds  for  the  algorithm. 
The  next  section  gives  an  experimental  evaluation  of  the  quality  of  these  algo¬ 
rithms. 

Subsection  4.2.2  gives  a  formulation  of  the  Phase  I  capacity  allocation  prob¬ 
lem  in  which  aggregated  link  capacity  is  allocated  to  the  call  originating  nodes. 
An  example  solution  to  the  capacity  allocation  problem  is  sketched  in  Figure  4.1 
where  two  different  shading  colors  show  capacities  allocated  to  call  originating 
nodes  A  and  G.  The  capacity  allocation  problem  is  a  linear  program  which  can 
be  solved  using  standard  commercial  software.  The  solution  provides  a  global 
capacity  allocation.  In  calculating  this  allocation,  the  total  bandwidth  allocated 
to  all  call  origination  nodes  on  a  given  link  might  include  both  bandwidth  cur¬ 
rently  available  and  bandwidth  that  must  be  freed  by  preempting  existing  calls. 
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capacity  allocated 
from  Ato  B,  C,  D  and  E 

capacity  allocated 
from  G  to  B,  E,  D  and  F 

free  capacity 


Destination 


nodes 


Figure  4.1:  Phase  I:  Capacity  Allocation  by  Solving  a  Linear  Program 
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In  the  second  phase,  each  call  origination  node  sets  up  calls  “assuming”  all  allo¬ 
cated  bandwidth  is  available.  When  each  individual  call  is  actually  set  up,  the 
nodes  along  the  call’s  path  will  preempt  existing  calls  if  necessary.  The  Phase  I 
problem  must  be  solved  once  by  a  single  central  node.  Possible  candidate  loca¬ 
tions  for  solving  this  problem  include  one  of  the  two  nodes  adjacent  to  the  failed 
link  and  the  network  control/management  center. 

Subsection  4.2.3  presents  an  approximate  algorithm  for  the  bulk  path  selec¬ 
tion  problem  which  is  solved,  in  parallel,  by  each  call  origination  node.  Figure  4.2 
demonstrates  the  way  by  which  an  aggregated  link  capacity  allocated  to  nodes  A 
and  G  is  assigned  to  individual  calls.  The  bulk  path  selection  problem  is  a  variant 
of  the  bandwidth  packing  problem.  The  algorithm  presented  involves  the  iter¬ 
ative  solution  of  shortest  path  problems  and  is  motivated  by  first-fit-decreasing 
bin  packing  heuristics. 

4.2.1  Network  Environment 

One  of  the  basic  notions  underlying  this  chapter  is  to  have  network  manage¬ 
ment  functions  and  network  control  functions  cooperate  on  this  problem.  Our 
paradigm  is  that  each  node  has  a  copy  of  the  network  topology  database  which 
has  such  information  as:  the  set  of  links  and  nodes  in  the  network,  current  load¬ 
ing  of  the  links  (or  a  moving  average  of  this)  perhaps  broken  down  by  preemption 
priority  or  by  QOS  class.  This  topology  database  is  updated  via  network  control 
protocols  and  is  used  as  input  to  the  route  computation  algorithms. 

In  addition,  each  link  tracks  the  bandwidth  currently  associated  with  the 
connections  crossing  it.  When  a  link  fails,  the  node  on  each  side  of  the  failing 
link  may  access  this  information  and  may  use  network  management  (or  other) 
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delimeters  of 
the  call  channels 


Originating 

nodes 


Figure  4.2:  Phase  II:  Bulk  Rerouting 


protocols  to  send  the  list  of  affected  calls  to  the  node  that  is  performing  the 
Phase  I  computation  (as  mentioned  previously,  the  node  performing  the  Phase 
I  calculation  could  be  one  of  the  nodes  adjacent  to  the  failed  link).  Note  that 
some  of  this  information  may  have  to  be  kept  for  network  accounting  and  charge- 
back  anyway,  but  it  may  also  have  some  real-time  network  control  uses.  It  is 
also  possible  to  aggregate  this  information  somewhat  as  suggested  later  in  this 
chapter.  In  any  case,  it  is  perfectly  reasonable  for  a  node  to  keep  such  information 
(IBM’s  NBBS  does  keep  much  of  this  information  in  the  switches  [1,  43]). 

The  time  required  to  manipulate  and  possibly  transmit  the  information  for 
the  Phase  I  computation  should  be  reasonable.  Specifically,  if  a  link  carrying- 
4096  connections  in  each  direction  fails,  then  the  nodes  adjacent  to  it  each  send 
4096  x  (calling  party  id,  preemption  priority  class,  and  traffic  parameters).  Cer¬ 
tainly  this  is  not  an  outrageous  amount  of  information  and  probably  less  than 
lOOIv  bytes.  The  node  performing  the  Phase  I  computation  distributes  (pos¬ 
sibly  via  multicast  directly,  or  possibly  through  the  topology  database  update 
mechanism)  the  results  of  Phase  I  to  the  rest  of  the  network  for  the  Phase  II 
computation. 

4.2.2  Phase  I:  Capacity  Allocation 

Below,  we  give  a  linear  programming  formulation  of  the  capacity  allocation  prob¬ 
lem.  We  introduce  the  following  variables: 

for  all  (i,  j)  e  Af  x  A/",  t  €  H  :  yfj  is  total  amount  of  priority  t  capacity 

assigned  to  calls  whose  origin,  destination 
pair  is  (/../): 
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for  all  e  G  A,  i  G  A/" :  i*  is  total  amount  of  capacity  assigned  on  arc  e 

to  calls  whose  origination  node  is  i. 


The  x'e  variable  values  are  the  capacity  allocations  that  are  passed  from  Phase  I 
to  Phase  II.  As  before  s\  is  the  cost  coefficient  which  discourages  long  paths. 

The  LP  problem  is  obtained  from  problem  (IPR)  by  aggregating  x  and  y 
variables  and  relaxing  the  integrality  constraints. 

(CA): 


Maximize 
subject  to 

E  xl  ~ 

eeOUT(j) 


E  A  E  Vij  ~  E  Ch(m)bwZm  -EE  (4-14) 

ieh  (ij)eJJ'xjJ'  meM  ieM  eeA 


E 

eelN(j) 


■e  = 


H  Vij 

A!  m  EX'  Vim 


Ex'1  <c  B  N 7  A  b  z 

^ e  —  e  i  /  ,  uemum/°m 


ieAf 


m£A4 


if  j  A  i 

for  all  i,j  G  M ,  (4.15) 

if  j  =  i 

for  all  e  G  A ,  (4.16) 


^  <  E  bk  for  all  i,j  G  J\f,  l  G  H,  (4.17) 

k:0(k)=i,V(k)=j,h(k)=£ 

x\,  yfj  >  0  for  all  i,j  G  A/",  e  G  A,  i  G  H ,  (4.18) 
0  <  2m  <  1  for  all  m  G  A4.  (4.19) 


It  can  be  easily  seen  that  the  formulation  (CA)  immediately  follows  from  (IPR) 
if  we  perform  the  substitutions: 

Xe  =  E  bkXe  t  (4-20) 

k:0(k)=i 

Vij  =  E  bkVk  (4-21) 

k:0(k)=i,D(k)=j,h(k)=£ 

and  relax  the  integrality  constraints. 

Because  of  the  discrete  nature  of  calls  (they  are  not  allowed  to  be  split),  it  is 
likely  that  the  Phase  II  algorithm  will  not  be  able  to  “pack”  all  calls  within  the 
capacity  allocated,  if  the  amount  of  the  capacity  allocated  is  exactly  equal  to  the 
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total  bandwidth  of  the  calls  to  be  packed.  To  compensate  for  this  effect,  at  least 
partially,  we  introduce  an  additional,  “virtual”  priority  class  Thus,  we  update 
the  set  of  priority  classes:  H'  =  H\ J{f},  and  define  Q/  <$;  min{Q  :  £  G  H}.  Then 
the  objective  function  (4.14)  becomes: 

Maximize  J2  ^  E  Vij  ~  Y  $ h(m)bmzm  YY  (4-22) 

'•  II'  (i,j) G,Vx,V  meM  /•  A' <  •  .  I 

A  new  heuristic  upper  bound  for  the  y- vari abl es  is  added: 

yfj  <  £  Y  bk  for  a]1  h3  e  N.  (4.23) 

k:0(k)=i,D(k)—j 

(e  =  .1  is  a  heuristically  set  parameter), 


E  £i~  Y  ■ 

eeOUT(j)  eelN(j) 


for  all  i,  j  e  A/",  (4.24) 


and  the  flow  conservation  constraints  (4.15)  are  changed  to: 

-  E eeH'  Vij  if  3  / 

Emm  EeeH'  Vim  if  3  =  ?- 
(i.e.  H  is  replaced  by  H')  and  we  obtain  a  new  formulation  (CA')  which  is  a 
slight  modification  of  (CA). 

Till  now,  it  was  assumed  the  set  of  currently  active  calls  to  be  known  explic¬ 
itly.  It  might  not  be  the  case  in  the  real-life  situations  since 


1.  this  set  of  calls  is  enormously  large: 

2.  the  algorithm  will  be  executed  by  each  individual  source  node  and  these 
nodes  generally  will  not  have  global  call  information  available. 


Instead,  we  will  assume  that  for  each  link,  in  addition  to  the  residual  capac¬ 
ity,  the  amount  of  bandwidth  allocated  to  each  priority  class  of  calls  is  known. 
Specifically,  we  have  for  each  link,  e  e  A,  and  priority  class  l  G  H, 

Ble  =  bandwidth  amount  along  link  e  occupied  by  calls  of  priority  class  L 
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We  then  define  variables  for  each  e  E  A  and  i  E  H, 


5^  =  the  amount  of  bandwidth  preempted  from  priority  class  t  on  link  e, 

Note  that  we  have  made  a  rather  significant  simplifying  assumption  that  the 
cost  of  preempted  traffic  is  a  linear  function  of  the  total  bandwidth  freed  on  each 
link.  This  assumption  represents  an  approximation  for  two  reasons: 

1.  bandwidth  is  freed  in  discrete  increments  corresponding  to  the  bandwidth 
associated  with  the  calls  that  must  be  preempted. 

2.  preempting  of  a  single  call  will,  in  general,  free  bandwidth  on  several  links. 


While  there  are  certainly  circumstances  under  which  this  assumption  clearly 
leads  to  significant  cost  distortions,  we  feel  that  this  simplification  (or  another 
similar  one)  is  necessary  in  order  to  derive  a  practical  approach  to  routing. 
Using  definitions  for  2  variables  and  Sem ,  one  can  express  B[  and  zle  as  follows: 


B[  =  Y.  Sembm 

m:  mEA4,h(m)=£ 

^  '  bemhm^m 

m: 


&  = 
~ e 


Next,  we  use  5  variables  to  formulate  the  LP  which  can  be  used  in  practical 
computations. 
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(CAP): 


Maximize  y  ce  y  4  /'IE  “EE  SIK 

t eH  (i,j)eNxN  eeAleH  ic-.W  .1 

subject  to 

(4.25) 

M 

Hi 

1 

M 

Hi 

II 

—  YUeH  Vij  if'./  /  ;  ,  n  .  . 

for  all  j  e  N , 

(4.26) 

eeOUT(j) 

eelN(j) 

Y^rneU  YUeH  Vim  if  3  =  « 

< 

y  x\  <  Be  +  y  z£e  for  all  e  G  A, 
ieAf  ten 

(4.27) 

iiii  <  E  bk  for  a11  e  A':  £  e  //. 

k:0(k)=i,D(k)=jMk)  =  t 

(4.28) 

£  <  Be 

for  all  e  €  A,  £  €  H, 

(4.29) 

;f|l  y-j,  zle>0  for  all  i,  j  G  A/”,  e  G  A,  P  G  H. 

(4.30) 

Notice  that  the  cost  of  call  preemption  in  this  formulation  may  be  over- 
counted  since  it  is  summed  up  over  several  links  the  call’s  bandwidth  was  freed. 
To  compensate  this  effect,  a  corrective  coefficient  p  was  applied.  Theoretically, 

1 

number  of  hops  in  the  max  hop  route  of  the  existing  call  ^ 

After  some  experimentation,  we  choose  p  =  0.5. 

Similar  to  defining  problem  (CA'),  substitution  of  priority  set  H  with  H' 
can  be  performed  and  then  we  obtain  formulation  (CAP') 

Problems  (CA)  (or  (CA'))  and  (CAP)  (or  (CAP'))  were  solved  with  a 
linear  programming  solver  for  the  variety  of  instances.  The  solutions  obtained 
indicated  that  the  approximation  was  sufficiently  accurate. 

Observe  that,  since  problems  formulated  in  this  section  are  essentially  multi- 
commodity  flow  problems,  specialized  algorithms  could  have  been  used.  Explo¬ 
ration  of  such  approaches  lies  beyond  the  scope  of  our  current  study. 

As  a  result  of  solving  the  problem  (CAP'),  we  assign  arc  capacities  to  each 
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origination  node.  Those  capacities  are  input  into  the  algorithm  described  in  the 
next  section. 

4.2.3  Phase  II:  Bulk  Routing 

In  this  section,  we  consider  an  algorithm  which  works  independently  for  all 
origination  nodes.  Each  node  must  receive  a  capacity  allocation  from  the  Phase 
I  algorithm.  That  capacity  allocation  for  call  origination  node  ?'0,  for  each  link 
e,  is  denoted  by  Be(i0).  This  link  capacity  is  set  equal  to  the  value  of  variable  xle 
computed  in  Phase  I.  The  problem  that  must  be  solved  in  Phase  II  is  a  variant 
of  FPAB  (given  in  Section  4.1.3)  in  which  the  zle  (call  preemption)  variables 
are  not  present.  That  is,  each  call  origination  must  carry  out  path  selection 
using  only  the  capacity  allocated  in  Phase  I.  The  paths  calculated  by  each  call 
origination  node  will  automatically  lead  to  appropriate  call  preemption  since  this 
was  planned  in  determining  the  Phase  I  allocations. 

Two  considerations  are  critical  in  designing  the  routing  algorithm:  first,  calls 
should  be  routed  over  shortest  paths  and,  second,  care  should  be  taken  in  “pack¬ 
ing”  the  calls  into  the  constructed  routes  as  tightly  as  possible.  In  addition,  the 
running  time  should  be  small  enough  to  allow  for  real-time  use.  To  address  these 
issues,  the  algorithm  designed  makes  use  of  shortest  path  algorithms  (see  e.g. 
[2])  and  handles  the  calls  in  the  spirit  of  the  so-called  first  fit  decreasing  (FFD) 
heuristic  which  is  known  to  be  quite  effective  for  the  bin  packing  problem  (see 
e.g.  [32]). 

We  arrange  the  calls  within  the  same  priority  class  in  descending  order  of 
their  bandwidth.  As  the  algorithm  proceeds,  we  keep  a  pointer  to  the  first 
unrouted  call,  k.  Focusing  on  A:,  we  execute  a  special  shortest  path  algorithm 
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which  constructs  a  shortest  path  tree  (SPT)  rooted  at  node  io  using  only  those 
arcs  which  have  enough  capacity  to  handle  k.  Afterwards,  we  try  to  route  as 
many  calls  as  possible  over  that  SPT  taking  them  in  the  order  they  appear  in 
the  list  (descending  by  bandwidth).  The  process  iterates  using  the  new  first 
unrouted  call.  When  we  can  not  route  anymore  calls  within  the  same  priority 
class,  we  proceed  to  the  next  priority  class.  Clearly,  this  approach  is  similar 
to  the  FFD  algorithm  for  the  bin  packing  problem.  Indeed,  the  computational 
results  presented  in  the  next  section  indicate  that  it  is  quite  effective  in  practice. 

Input  parameters  of  the  algorithm  are  network  (N,A),  the  origination  node 
i o,  allocated  capacities  on  the  arcs  Be(i0 ),  e  G  A  and  set  of  calls  /Q0  =  {A:  G 
JC\0(k)  =  ?'o},  output  is  the  set  of  routes.  We  suppose  the  priority  levels  are 
ordered  in  the  descending  order:  C  >~  h  >~  ■  ■  ■  >~  1\h\  and  denote  A(w)  =  {e  |  e  G 
A.  B,  >  w}.  To  better  access  the  largest  bandwidth  calls,  all  calls  are  grouped 
by  their  destinations  and  stored  at  corresponding  binary  trees.  Another  set  of 
internal  data  consists  of  capacity  availability  labels  w.j  stored  at  each  node  j 
which  show  the  maximum  amount  of  capacity  that  can  be  used  to  route  a  call  to 
the  node  or,  equivalently,  it  is  the  capacity  of  the  “narrowest”  arc  on  the  path 
from  i0  to  j  in  the  SPT.  In  general  the  arc  lengths  defining  the  shortest  path 
are  call  delays.  However,  in  this  case  we  consider  the  min-hop  path  and  hence 
we  have  unit  arc  lengths. 

1. for  priority  level  t  =  t\,  £2,  ■  ■ . ,  1\h\  do 

2.  Initialize  binary  trees  7),  Vj  G  A f: 

3.  Tj  contains  items  {k  G  /Q0  |  h(k)  =  t,  D(k )  =  j }  with  keys 

4.  while  there  is  a  non-empty  tree  Tj  do 

5.  find  call  k  :  b~k  =  maxjG_,v'  maxj£7-  b 


46 


6.  Construct  an  SPT  rooted  at  iq  for  the  network  (JV,  A{bk))] 

7.  if  there  is  a  path  from  i0  to  D(k ),  route  the  call  k  and  update  the  labels  Wj 

8.  Remove  the  call  k  from  its  tree:  Tk  <—  To(k) 

9.  do 

10.  find  call  k  :  b~k  =  maxjeA/- max^^. {bk  \  bk  <  Wj }; 

11.  if  A:  was  found, 

12.  route  the  call  k  in  the  current  SPT; 

13.  remove  the  call  k  from  its  tree:  Tk  <-  ^D(k) 

14.  update  the  labels  Wj  Vj  G  J\f; 

15.  while  there  is  a  call  A:; 

The  next  issue  we  want  to  discuss  here  is  the  running  time  of  the  algorithm. 
All  operations  we  perform  with  binary  trees  require  time  logarithmic  on  the 
size  of  the  tree.  These  operations  are:  removing  the  items  (lines  8  and  13)  and 
finding  the  max  key  item  (lines  5  and  10).  Clearly,  each  call  is  considered  only 
once.  Indeed,  if  there  is  enough  bandwidth  in  the  existing  SPT  for  the  call  (line 
11),  it  is  routed  and  removed  (lines  12-13):  if  there  are  no  calls  which  can  be 
routed  over  the  existing  SPT  (line  15),  a  shortest  path  algorithm  is  run  to  route 
the  largest  bandwidth  unrouted  call  (line  6).  If  there  is  a  path  for  that  call,  it  is 
routed  and  removed,  otherwise,  it  is  removed  anyway,  since  there  is  no  chance  for 
the  call  to  be  routed  later.  Capacity  availability  labels  Wj  are  updated  in  time 
0( |A/"| )  (lines  7  and  14).  The  time  required  to  extract  the  largest  bandwidth 
call  is  0(  |A/"|  log  |/C|).  If  we  denote  the  time  required  to  construct  an  SPT  by 
Tspt ,  the  worst  case  running  time  of  the  algorithm  is  0(|/C|(|A/"|  log  |/C|  +Tspt))- 
The  average  should  be  significantly  better  since,  in  general,  several  calls  will  be 
routed  over  each  SPT  constructed.  The  running  time  for  the  shortest  path 
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Call 

Bandwidth 

Priority 

Percentage 

Type 

Distribution 

Distribution 

Mix  1 

Mix  2 

Mix  3 

voice 

all  are  .064  Mb/s 

10%  are  1,  90%  are  2 

54% 

47% 

41% 

data 

50%  are  .0096  Mb/s,  50%  are 

uniformly,  1-4 

43% 

47% 

53% 

uniformly  from  .01  to  1.  Mb/s 

video 

uniformly,  1  to  3  Mb/s 

uniformly,  1-3 

3% 

6  % 

6% 

Table  4.1:  Call  Types 

algorithm  can  vary  depending  on  the  approach  used  and  the  exact  nature  of  the 
cost  function.  In  the  specific  context  we  considered,  since  the  objective  was  a 
pure  “min-hop”,  a  breadth  first  search  algorithm  could  be  used  which  runs  in 
time  0(\A\). 

4.3  Experimental  Results 

The  purposes  of  our  experiments  were: 

1.  estimate  the  quality  of  our  solution 

2.  estimate  the  time  required  to  perform  the  computations 

Our  sample  data  were  obtained  by  following  a  general  simulation  routine  de¬ 
scribed  in  chapter  3.  In  our  experiments  we  generated  calls  of  three  different 
types:  voice,  data  and  video.  Table  4.1  gives  bandwidth  and  priority  charac¬ 
teristics  for  the  calls.  Each  call  of  priority  class  1,  2,  3  and  4  was  assumed  to 
have  revenue  1.0,  0.1,  0.01  and  0.001  respectively.  Table  4.1  also  presents  three 
different  mixes  of  calls  which  were  considered. 
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The  network  used  for  our  simulation  had  a  topology  similar  to  nation-wide 
ARPA  net  with  59  nodes,  71  undirected  links  (respectively,  142  directed  arcs), 
as  it  appears  in  [10],  p.  135.  Each  link  had  a  full  duplex  capacity  of  154  Mb/s. 
The  overall  parameter  setting  was  chosen  so  that  after  a  link  failure  the  total 
number  of  calls  to  be  rerouted  would  be  between  300  and  1,600. 

We  evaluate  our  approach  in  two  general  areas.  First,  there  are  the  rev¬ 
enue/cost  components  reflected  in  the  objective  functions  previously  described, 
i.e.  the  total  revenue  of  all  calls  rerouted  minus  lost  revenue  of  calls  preempted. 
The  second  criterion  to  be  considered  involves  call  refusal.  A  call  refusal  occurs 
when  the  initial  setup  attempt  fails.  A  call  refusal  causes  a  substantial  cost 
in  overall  network  overhead  as  well  as  delay  associated  with  reestablishing  the 
call.  To  accurately  evaluate  these  concerns,  a  detailed  simulation  is  required.  In 
lieu  of  such  an  analysis,  we  carried  out  certain  experiments  which  indicate  how 
alternate  approaches  would  perform  under  “extreme”  conditions.  These  were 
sufficient  to  derive  definitive  conclusions. 

The  experimental  evaluation  consisted  of  two  parts.  In  the  first  part  the 
quality  of  the  call  routing  algorithm  was  estimated  without  considering  call  pre¬ 
emption.  In  the  second  part  the  evaluation  of  call  preempting  capabilities  was 
performed.  These  parts  of  the  experiment  are  described  below  in  sections  4.3.1 
and  4.3.2  respectively. 

4.3.1  Evaluation  of  the  Call  Rerouting  Algorithm 

In  lieu  of  simulation,  we  use  two  methods  to  establish  a  baseline  solution  for 
which  to  compare  the  quality  of  our  two-phase  approach  to  call  rerouting  with¬ 
out  preemption.  In  both  approaches,  each  call  is  rerouted  by  an  independent 
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execution  of  a  path  selection  (shortest  path)  algorithm.  For  the  first  approach, 
which  can  be  considered  a  “worst  case”  scenario,  we  assume  that  each  call  orig¬ 
ination  node  reroutes  its  calls  but  receives  no  updated  information  about  link 
status  during  the  rerouting  process.  Thus,  all  nodes  assume  that  they  have  full 
access  to  all  available  bandwidth  capacity.  This  is  called  “fast”  rerouting.  This 
approach  will  necessarily  lead  to  high  levels  of  call  refusal  in  that  competing  call 
origination  nodes  will  try  to  setup  calls  over  the  same  sets  of  links.  If  a  call  is 
refused,  we  do  not  attempt  to  reroute  it  -  another  very  pessimistic  assumption. 

Under  the  second  approach,  we  assume  that  a  central  node  with  perfect 
information  reroutes  all  calls.  That  is,  whenever  a  route  is  calculated,  the  precise 
status  of  all  links  is  known  and  as  a  result,  there  is  no  call  refusal.  The  central 
node  uses  a  first  fit  decreasing  heuristic  approach,  by  ordering  calls  in  decreasing 
order  of  bandwidth.  We  call  this  method  “slow”  rerouting  since  time  would  be 
required  to  propagate  the  information  about  the  state  of  the  network  and  since 
the  actual  implementation  of  this  approach  would  require  the  distribution  of  all 
routes  from  the  central  node.  The  two  approaches  constitute  two  extremes  in 
that 


•  the  fast  approach  requires  the  minimal  amount  of  information  exchange, 
but  generally  would  achieve  low  quality  routes; 

•  the  slow  approach  should  achieve  the  highest  quality  results  (subject  to  the 
limitations  of  the  first  fit  decreasing  approach  to  packing)  but  is  impractical 
due  to  the  extreme  information  exchange  requirements. 

Our  approach  can  be  viewed  as  a  compromise  which  executes  an  approximate 
global  rerouting  with  modest  information  exchange  requirements  while  keeping 
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final  path  selection  in  the  hands  of  the  call  origination  nodes. 

Two  sets  of  results  are  presented.  The  first  is  aimed  at  evaluating  the  effec¬ 
tiveness  of  the  overall  two  phase  approach  and  the  second  specifically  addresses 
the  Phase  II  packing  algorithm.  Table  4.2  gives  the  characteristics  of  the  calls 
to  be  rerouted.  For  each  utilization  and  call  mix  combination  the  number  of 
calls  to  be  rerouted  and  their  total  revenue  are  given.  The  results  presented  in 
Table  4.3  address  the  overall  approach.  The  first  four  rows  give  the  results  of 
the  fast  and  slow  rerouting  algorithms  respectively.  The  fifth  and  sixth  rows 
give  the  results  of  the  two  phase  approach  applied  to  the  same  problems.  Note 
that  it  shows  significant  improvement  over  the  fast  approach  and  comes  close 
to  achieving  the  quality  levels  of  the  slow  approach.  Note  that  the  results  do 
degrade  somewhat  for  the  highest  utilization  level.  This  can  be  expected  since 
for  these  levels  packing  becomes  more  difficult.  In  the  next  two  rows  results  are 
given  when  preemption  is  allowed.  This  gives  an  indication  of  the  rather  sig¬ 
nificant  potential  improvements  provided  by  preemption.  The  summary  of  the 
results  is  plotted  in  Figure  4.3.  Each  of  the  three  parts  in  this  figure  corresponds 
to  a  particular  distribution  of  the  call  types.  Pie  diagrams  on  the  right  show  the 
fractions  of  the  number  of  calls  in  each  class.  Plots  on  the  left  show  the  fraction 
of  the  rerouted  calls. 

Evaluation  of  the  call  preemption  algorithm  will  be  presented  in  section  4.3.2. 

The  final  row  in  Table  4.3  gives  the  running  times  in  seconds.  We  temper  the 
discussion  given  below  by  noting  that  the  experiments  were  run  on  a  Sun  Sparc 
10  which  may  or  may  not  be  similar  to  the  hardware  environment  anticipated 
in  an  actual  telecommunications  network.  The  running  time  of  the  first  phase  is 
always  lies  between  20  seconds  and  1  minute.  This  would  certainly  not  be  fast 
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Mix  1 

Mix  2 

Mix  3 

Utilization 

30% 

50% 

65% 

30% 

50% 

65% 

30% 

50% 

65% 

Total  call  revenue 

32.46 

79.50 

115.28 

49.79 

97.65 

124.99 

38.29 

114.40 

139.38 

Number  of  calls 

357 

897 

1566 

305 

646 

1414 

287 

670 

1375 

Table  4.2:  Call  Characteristics 


Total  Rev.  and  Number 

Mix  1 

Mix  2 

Mix  3 

of  the  Rerouted  Calls 

30% 

50% 

65% 

30% 

50% 

65% 

30% 

50% 

65% 

“Fast” 

revenue 

13.76 

33.73 

13.67 

32.37 

42.44 

10.66 

16.67 

52.78 

15.61 

Scenario 

number 

237 

503 

227 

229 

360 

162 

191 

359 

227 

“Slow” 

revenue 

32.46 

79.48 

56.25 

49.79 

97.64 

32.72 

38.29 

114.26 

48.72 

Scenario 

number 

357 

869 

148 

305 

630 

171 

287 

636 

179 

2-Phase  w/out 

revenue 

32.18 

78.03 

47.76 

49.56 

95.83 

32.07 

38.03 

109.25 

40.54 

preemption 

number 

355 

882 

306 

304 

620 

254 

284 

648 

280 

2-Phase  with 

revenue 

32.18 

79.16 

110.05 

49.56 

95.83 

119.53 

38.28 

110.80 

135.02 

preemption 

number 

355 

882 

1287 

304 

620 

1118 

286 

653 

1076 

Running 

Phase  I 

26.8 

28.5 

30.6 

25.1 

67.1 

34.7 

24.5 

55.2 

39.7 

Time 

Phase  II 

0.25 

0.60 

0.96 

0.25 

0.53 

0.88 

0.25 

0.48 

0.88 

Table  4.3:  Evaluation  of  the  Two-Phase  Approach 
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enough  for  real  time  operation  but  is  not  unreasonable  for  the  scenario  where 
the  network  management  system  had  provided  advance  warning  of  a  failure.  In 
order  to  achieve  real-time  response,  alternate  approaches,  including  heuristics, 
for  solving  the  LP  could  be  considered.  The  running  times  for  the  second  phase 
algorithm  were  always  less  than  1  second.  The  times  reported  consisted  of  the 
time  used  to  execute  sequentially  the  Phase  II  algorithm  at  each  call  origination 
node.  Of  course,  in  an  actual  implementation  these  calculations  would  be  carried 
out  in  parallel  leading  to  much  smaller  times.  We  feel  these  times  are  quite 
encouraging  in  that  they  would  seem  to  allow  for  real-time  operation. 

Table  4.4  presents  results  concerning  the  quality  of  the  Phase  II  algorithm. 
Three  sets  of  experiments  were  carried  out.  In  the  first  two,  the  fast  and  slow 
rerouting  algorithms  were  each  run.  The  capacity  taken  up  on  each  link  for  calls 
associated  with  each  call  origination  node  were  then  allocated  to  that  node.  The 
Phase  II  algorithm  was  then  given  the  job  of  routing  calls  using  only  the  capacity 
allocated.  This  tested  the  ability  of  the  phase  II  algorithm  to  choose  a  good  set 
of  calls  and  also  to  pack  them  into  the  bandwidth  restricted  links.  The  results 
for  this  experiment  are  given  in  the  first  8  rows  of  Table  4.4.  Note  that  the 
phase  II  algorithm  generally  did  significantly  better  than  the  fast  algorithm  and 
generally  came  close  to  achieving  the  results  of  the  phase  II  approach.  In  the 
9th  and  10th  rows,  the  LP  value  is  compared  with  the  value  achieved  in  Phase 
II  (the  Phase  I  value  is  an  upper  bound  on  the  Phase  II  value).  Note  that  for 
the  medium  and  low  utilization  cases,  the  Phase  II  algorithm  came  very  close 
to  achieving  the  LP  value.  For  the  high  utilization  it  was  not  as  close.  This, 
again,  can  be  explained  by  the  fact  that  for  these  levels  packing  becomes  more 
difficult.  Of  course,  further  investigations  would  be  required  to  determine  how 
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Total  Revenue  and  Number 

Mix  1 

Mix  2 

Mix  3 

of  the  Rerouted  Calls 

30% 

50% 

65% 

30% 

50% 

65% 

30% 

50% 

65% 

“Fast” 

revenue 

13.76 

33.73 

13.67 

32.37 

42.44 

10.66 

16.67 

52.78 

15.61 

Scenario 

number 

237 

503 

227 

229 

360 

162 

191 

359 

227 

Phase  II 

revenue 

14.90 

37.67 

16.54 

32.19 

46.35 

10.68 

18.14 

60.64 

20.53 

Algorithm 

number 

228 

496 

234 

220 

353 

167 

172 

388 

226 

“Slow” 

revenue 

32.46 

79.48 

56.25 

49.79 

97.64 

32.72 

38.29 

114.26 

48.72 

Scenario 

number 

357 

869 

148 

305 

630 

171 

287 

636 

179 

Phase  II 

revenue 

32.20 

76.41 

56.21 

48.95 

96.33 

32.07 

38.03 

109.25 

48.42 

Algorithm 

number 

312 

802 

161 

271 

582 

176 

252 

597 

183 

Phase  I  revenue  (LP  value) 

32.46 

79.50 

56.29 

49.79 

97.65 

32.77 

38.29 

114.40 

50.27 

Phase  II  revenue 

32.18 

78.03 

47.76 

49.56 

95.83 

26.21 

37.35 

106.80 

40.54 

Table  4.4:  Effectiveness  of  the  Phase  II  Packing  Algorithm 

much  of  the  gap  is  due  to  the  quality  of  the  relaxation  and  how  much  is  due  to 
the  quality  of  the  heuristic. 


4.3.2  Evaluation  of  the  Call  Preemption 


In  this  section  we  present  an  evaluation  of  call  preemption  algorithm.  Suppose 
after  the  work  of  the  two-phase  algorithm  was  completed  a  set  of  rerouted  calls 
}Cr  was  obtained  together  with  their  routes.  Clearly,  /C,.  C  1C.  The  set  of  the 
new  routes  can  be  described  by  introducing  notation  similar  to  (4.1): 


for  all  k  <G  KLr.  e  <G  A:  aek  = 


1  if  call  k  is  carried  by  arc  e 
0  otherwise 


(4.31) 
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Now,  consider  the  following  IP  formulation: 


Maximize 

'y  '  bh(k)bkUk  y  ' 

(4.32) 

k£  Kr  m£. M 

subject  to 

y  '  ^ekbkVk  y  'j  bembmZm  A  Be 

for  all  e  e  A, 

(4.33) 

kElCr  meM 

0  <  Vk  <  1  and  integer 

for  all  k  G  ICr , 

(4.34) 

0  <  zm  <  1  and  integer 

for  all  m  G  A4. 

(4.35) 

It  can  be  easily  seen  that  solution  to  (4.32-4.35)  is  a  feasible  suboptimal  solution 
to  (IPR)  corresponding  to  the  call  rerouting  and  preemption  given  by  the  two- 
phase  algorithm.  The  output  of  this  integer  program  (4.32-4.35)  consists  of 
decisions  regarding  call  refusal  (the  variables)  and  call  preemption  (the  zm 
variables).  Thus,  in  a  certain  sense,  it  gives  the  optimal  solution  to  to  the  call 
setup  process  which  is  carried  in  a  distributed  fashion  by  the  network  nodes. 
Then  the  quality  of  the  two-phase  algorithm  can  be  estimated  by  how  close 
solution  of  (4.32-4.35)  is  to  that  of  (IPR).  As  it  was  noted  above,  the  latter 
problem  can  not  be  solved  explicitly  but  it  is  possible  to  estimate  the  value  of  its 
solution  from  above  by  solving  its  LP  relaxation  which  is  essentially  (CA).  The 
results  of  the  comparison  for  the  previously  described  set  of  calls  are  presented 
in  the  table  4.5.  We  consider  these  results  to  be  quite  encouraging. 

4.4  Conclusions 

In  this  chapter  we  have  presented  models  and  algorithms  which  address  certain 
rerouting  problems  which  arise  in  multi-class  traffic  environments.  The  approach 
has  potential  applicability  in  several  contexts.  Our  computational  experiments 
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Calls  Rerouted 

Mix  1 

Mix  2 

Mix  3 

30% 

50% 

65% 

30% 

50% 

65% 

30% 

50% 

65% 

Upper  Bound 

32.46 

79.50 

113.24 

49.79 

97.65 

121.64 

38.29 

114.40 

136.27 

Lower 

value 

31.58 

74.94 

109.34 

47.49 

93.89 

112.59 

37.54 

106.40 

132.93 

Bound 

number 

318 

835 

1245 

276 

588 

1084 

254 

604 

1033 

Gap 

2.8% 

6.1% 

3.6% 

4.8% 

7.5% 

2.5% 

Table  4.5:  Preemption  Evaluation 

have  demonstrated  that  these  methods  can  greatly  improve  the  response  of  a 
telecommunications  system  to  link  failures  over  more  direct  approaches.  In  order 
to  carry  out  the  implementation  our  methods,  specific  protocols  must  be  designed 
to  transfer  the  required  information.  Furthermore,  refinements  must  be  made  to 
insure  appropriate  response  in  a  real-time  setting. 
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Chapter  5 


Virtual  Path  Layout  in  an  ATM 
Environment 

Introduction 

The  Virtual  Path  (VP)  is  an  important  concept  in  ATM  multiservice  networks. 
In  our  research  framework  we  consider  VPs  as  logical  macrolinks  in  telecommu¬ 
nication  networks.  That  is,  in  general,  a  VP  may  traverse  several  physical  links 
in  the  network  to  provide  a  faster  connection  between  a  pair  of  remote  network 
nodes.  This  efficiency  results  from  a  greatly  reduced  processing  time  in  the  in¬ 
termediate  switching  nodes  of  the  VP.  However,  there  is  a  tradeoff  involved  since 
capacity  is  allocated  for  each  VP  individually.  In  various  ATM  network  models, 
the  allocated  capacity  is  treated  differently  but  most  often  it  is  considered  not 
to  be  available  for  use  by  the  other  VPs  on  the  link.  This  in  turn  reduces  the 
gain  of  statistical  call  multiplexing. 

In  the  literature,  there  exist  several  different  approaches  to  modeling  and 
solving  the  VP  layout  problem.  In  [5,  18],  the  routing  problem  in  VP-based 
ATM  networks  is  considered.  Dynamic  routing  policies  suggested  by  the  authors 
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attempt  to  minimize  the  VC  blocking  rate.  Simulations  were  used  to  evaluate 
the  proposed  routing  policies. 

References  [9,  30]  consider  the  problem  of  minimizing  the  total  VC  blocking 
rate  where  VC  demand  between  each  node  pair  is  assumed  randomly  arriving 
with  a  Poisson  distribution.  There  are  two  special  cases  of  the  problem  solved  by 
the  authors:  (1)  each  VC  uses  exactly  one  VP  and  (2)  each  VP  involves  exactly 
one  physical  link.  The  problem  was  formulated  as  a  nonlinear  non-differentiable 
combinatorial  optimization  problem.  Suggested  approaches  include  Lagrangean 
relaxation,  penalty  functions  and  subgradient  optimization  techniques. 

Reference  [3]  explored  various  issues  of  ATM  network  restoration,  particularly 
rerouting  strategies  emphasizing  restoration  speed  and  capacity  allocation  and 
redundant  capacity  allocation  methods  that  take  into  account  that  a  cut  of  a 
single  physical  link  may  result  in  multiple  service  failures.  However,  analytical 
models  for  the  call  rerouting  or  spare  capacity  allocation  were  not  proposed. 

References  [15,  14]  deal  with  optimizing  dynamic  virtual  path  bandwidth 
allocation  and  restoration  schemes.  It  is  assumed  that  VP-level  protocols  provide 
the  optimal  logical  assignment  of  VP  bandwidth  for  an  upcoming  control  period. 
The  objective  was  to  minimize  the  total  rejected  bandwidth  demand  subject 
to  QOS  constraints  (cell  loss  and  delay)  and  100%  restorability.  The  latter 
requirement  implied  that  there  should  be  a  certain  amount  of  spare  capacity 
kept  to  respond  to  link/node  failures.  The  authors  did  not  consider  explicit 
backup  routing  algorithms. 

In  reference  [33] ,  a  proposed  approach  was  described  to  minimize  the  amount 
of  expected  traffic  flow  loss  when  a  link  failure  occurs.  It  was  assumed  that 
the  traffic  is  restored  using  a  set  of  routes  bypassing  the  failed  link.  In  this 
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approach  the  bypassing  routes  begin  and  terminate  at  the  end  nodes  of  the  failed 
link.  Statistical  multiplexing  effects  which  would  affect  the  required  amount  of 
capacity  when  the  VPs  are  split  were  not  considered. 

The  VP-layout  methods  proposed  in  our  paper  allow  fast  restoration  for  ATM 
networks  based  on  the  concept  of  backup  VPs.  This  approach  was  suggested 
in  series  of  papers  [21,  20,  22],  We  follow  the  proposed  framework  and  build 
an  optimization  model  for  VP  layout  satisfying  certain  reliability  constraints. 
Essentially,  we  take  advantage  of  the  following  ATM  feature.  It  is  known  that 
the  flexibility  of  the  ATM  architecture  makes  it  possible  to  establish  a  zero 
capacity  backup  VP  and  then,  in  the  event  of  failure,  the  full  capacity  can 
be  captured  sufficiently  fast  to  provide  a  real-time  service  restoration.  This 
approach  has  the  advantage  that  the  time-consuming  process  of  setting  routing 
tables  at  the  switches  and  defining  VP  identifiers  (VPI)  can  be  carried  out  in 
an  off-line  mode  and,  when  a  failure  occurs,  an  affected  primary  VP  can  be 
switched  to  a  backup  without  noticeable  interruption  of  real-time  services.  The 
restoration  algorithm  is  very  simple  [22],  Basically,  it  consists  of  failure  detection 
followed  by  a  restoration  message  sent  along  the  VP  to  reserve  a  proper  amount 
of  capacity  (this  algorithm  is  outlined  in  more  detail  in  section  2.2). 

We  consider  four  variants  of  the  fault-tolerant  VP  layout  based  on  combina¬ 
tions  of  the  following  two  alternatives: 

1.  (a)  Exactly  one  link-disjoint  backup  VP  is  established  for  any  primary 

VP; 
or 

(b)  a  set  of  backup  VPs  is  established  for  a  primary  VP  so  that  each 
backup  VP  provides  a  backup  in  the  event  of  the  failure  of  a  corre- 
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Figure  5.1:  An  Example  of  Two  Rerouting  Strategies 
sponding  link  of  the  primary  VP. 

2.  (a)  The  amount  of  bandwidth  captured  on  the  backup  VP  during  the 

restoration  is  exactly  the  same  as  that  of  the  primary  VP; 
or 

(b)  the  backup  VP  is  allocated  less  than  the  total  primary  VP  bandwidth. 

(In  the  latter  case  a  set  of  VCs  must  be  dropped  when  the  failure  occurs.) 

The  approach  which  is  the  simplest,  at  least  conceptually,  is  to  establish  one 
backup  VP  for  each  primary  VP  (1(a)).  This  backup  VP  would  be  used  if  any 
link  in  the  primary  VP  fails.  Thus,  it  should  be  link  disjoint  from  the  primary 
VP.  Alternatively,  the  backup  route  may  be  set  up  to  detour  a  certain  link  in 
a  primary  VP  (1(b)).  Then  there  are  several  backups  for  each  primary  VP, 
generally  as  many  as  there  are  links  traversed  by  the  primary  VP.  Backups  of 
this  kind  serve  as  bypasses  around  faults.  An  example  of  these  two  strategies  is 
shown  in  Figure  5.1.  Here,  a  primary  VP  is  set  over  the  route  ABC.  Under  the 
first  strategy,  a  link-disjoint  backup  VP  p  is  used  if  any  of  the  links  AB  or  BC 
fails.  Under  the  second  strategy,  a  backup  VP  consisting  of  path  q  and  link  BC 
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is  used  when  link  AB  fails  and  a  backup  VP  consisting  of  links  AB  and  path  r 
is  used  when  link  BC  fails. 

Naturally,  one  should  expect  less  redundant  capacity  to  be  required  in  the 
network  when  implementing  the  second  strategy.  However,  it  has  a  certain 
disadvantage.  Each  link  in  the  ATM  network  may  carry  a  number  of  VPs  not 
exceeding  4096  (=  212)  which  is  limited  by  the  12  bits  reserved  for  the  VPI  in 
the  cell  header.  In  the  case  of  large  networks  and  a  large  number  of  long  VPs, 
there  might  not  be  enough  VPI  values  to  set  up  all  link-dependent  backup  VPs. 
Also  in  this  case  the  control  overhead  could  be  significantly  higher. 

In  this  chapter  we  consider  a  mathematical  model  for  case  la/2a.  Our  basic 
approach  was  motivated  by  the  solution  to  a  similar  problem  for  packet  net¬ 
works  proposed  by  Gavish  et  al  [12],  However,  their  model  contains  a  group  of 
non-linear  constraints  which  causes  certain  difficulties  when  a  linear  relaxation 
technique  is  applied.  We  show  how  these  constraints  can  be  linearized. 

Considering  the  second  pair  of  alternatives,  a  backup  VP  need  not  necessarily 
have  the  same  capacity  as  the  corresponding  primary  VP.  Lesser  capacity  could 
be  allocated  in  the  case  of  limited  network  resources.  When  such  a  backup 
VP  is  used,  some  lower-priority  calls  would  have  to  be  taken  down  or  rerouted 
separately.  In  order  to  accomplish  this,  the  network  protocols  should  have  the 
ability  to  take  down  or  reroute  individual  calls  or  VCs. 

As  was  stated  in  chapter  2,  we  also  suggest  a  more  general  model  for  fault- 
tolerant,  link-disjoint  VP  layout  where  a  VP  can  serve  as  both  primary  and 
backup  VP.  In  this  case  we  allow  traffic  flow  to  both  VPs  and,  when  failure  oc¬ 
curs,  traffic  is  switched  from  one  VP  to  the  other.  The  advantage  of  this  model 
is  additional  routing  flexibility  and  higher  throughput  as  will  be  shown  by  the 
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example  of  section  5.2.5.  In  chapter  6,  we  develop  a  two-step  heuristic  solution 
approach  for  this  model.  On  the  first  step,  a  VP  pair  is  routed  between  every 
node  pair  (using  link-disjoint  routes)  based  on  the  solution  of  a  multicommod¬ 
ity  flow  problem;  on  the  second  step,  a  sort  of  dimensioning  problem  is  solved 
assigning  appropriate  capacities  to  VPs  and  reserving  backup  capacities  on  the 
network  physical  links. 

Another  performance  measure  which  we  consider  in  our  research  is  adapt¬ 
ability  of  a  particular  VP  layout  to  changes  in  traffic  demand.  Quantitatively, 
we  measure  adaptability  in  terms  of  call  blocking  probability.  Our  conjecture 
is  that  longer,  lower  capacity  VPs  have  lower  adaptability  than  shorter,  higher 
capacity  VPs.  In  order  to  validate  this  conjecture,  we  simulate  the  random  gen¬ 
eration  of  calls  which  are  routed  through  the  network  using  previously  computed 
VPs.  Our  results  show  an  increase  in  call  blocking  probability  when  VPs  become 
longer  but  more  narrow. 

5.1  Statistical  Multiplexing 

We  start  our  study  of  the  VP  layout  by  building  a  very  simple  analytical  model 
for  statistical  multiplexing  of  VBR  calls  which  captures  the  main  properties  of 
statistical  multiplexing  effect  and  gives  a  linear  or  piece-wise  linear  approxima¬ 
tion  suitable  for  our  purposes. 

We  consider  the  two  most  common  VBR  traffic  types  in  telecommunication 
networks:  video  and  data.  From  the  literature  [35,  45],  it  is  known  that  the 
distribution  of  video  traffic  transmission  rates  can  be  approximated  by  a  normal 
distribution.  Data  traffic  is  usually  considered  as  generated  by  on/off  sources 
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with  alternating  periods  of  silence  and  transmission  at  peak  rate. 

We  assume  that  buffer  sizes  in  the  network  are  minimal.  This  is  the  case  for 
most  high  speed  networks  aimed  at  providing  real-time  service.  In  fact,  queues 
are  used  only  to  avoid  cell  loss  when  arrivals  from  several  multiplexed  sources 
occur  simultaneously.  For  example,  even  when  streams  from  CBR  sources  are 
multiplexed  together  at  a  node  and  the  outgoing  link  transmission  rate  is  greater 
than  the  total  cell  arrival  rate,  there  is  a  chance  that  instantaneous  incoming  rate 
may  exceed  the  transmission  rate  when  current  cells  on  all  links  arrive  within 
a  short  time  interval.  However,  even  in  this  case,  there  is  no  real  cell  queuing- 
delay  since  the  cells  in  the  queue  are  transmitted  before  the  arrival  of  the  next 
group  of  cells. 

Notice  that  this  model  presents  an  essential  simplification  when  compared 
with  other  models  (c.f.  [39])  where  buffer  size  is  an  essential  parameter.  Taking 
buffer  size  out  of  consideration  makes  our  model  insensitive  to  the  distribution 
of  the  active  period  lengths.  For  example,  adopting  this  model  we  can  not 
distinguish  between  Poisson  and  self-similar  types  of  traffic.  Though  it  is  an 
essential  simplification  of  the  model,  it  allows  us  to  restrict  our  consideration  to 
the  probability  distribution  of  the  transmission  rate. 

We  define  our  capacity  allocation  rule  as  providing  the  lowest  capacity  suf¬ 
ficient  to  yield  a  predefined  CLR,  e,  when  a  set  of  VBR  connections  1C  is  mul¬ 
tiplexed  together.  Let  Ak  be  a  random  variable  describing  transmission  rate  of 
call  k  e  1C  and  A  =  JCkeic  A k  be  the  total  transmission  rate.  Then  capacity 
allocation  rule  providing  CLR,  e,  can  be  formally  defined  as: 

C(K)  =  inf{C  :  Pr{A  >  C)  <  s}  (5.1) 

In  other  words,  we  want  the  probability  of  exceeding  the  allocated  capacity  to 
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be  less  than  predefined  CLR.  Based  on  this  definition,  we  will  derive  piece-wise 
linear  approximations  for  typical  transmission  rate  distributions  encountered  in 
telecommunication  networks. 

5.1.1  Video  Traffic  Capacity  Allocation 

In  this  section  we  derive  a  piece-wise  linear  approximation  for  the  capacity  al¬ 
location  rule  in  the  case  of  video  traffic  which  is  assumed  to  have  a  normally 
distributed  transmission  rate  [45]. 

Suppose  n  video  connections  are  multiplexed  together.  Then  their  total  trans¬ 
mission  rate  A  is  a  normally  distributed  r.v.  with  mean  //,  =  J2j Li  Ah  and  variance 
a2  =  Ya= 1  af  where  //.,  and  of  are  mean  and  variance,  respectively,  of  the  i-th 
connection.  Taking  into  account  properties  of  the  normal  distribution,  we  obtain: 

Pr{A  >C)  =  Pr{Z  >  (C  -  fi) /a)  (5.2) 

where  Z  is  a  random  variable  having  a  standard  normal  distribution. 

From  (5.1)  and  (5.2) 

C  =  n  +  ze<7  (5.3) 

Here,  z£  is  1—s  quantile  for  standard  normal  distribution.  It  can  be  approximated 
[35]  as  z£  ~  \J2  ln(l/t)  —  ln(27r). 

Thus,  in  case  of  n  video  calls  with  independent  identically  distributed  trans¬ 
mission  rates  with  ^  =  fi0  and  cq  =  <j0,  we  have  the  following  capacity  allocation 
rule 

C(n)  =  no  n  +  z£a0\/n.  (5.4) 

For  the  purposes  of  our  approximation  we  assume  that 
A<o  =  kuo 
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where  k  is  a  constant  between  3  and  4  chosen  such  that  the  probability  to  have 
zero  traffic  would  be  low  but  positive.  Then,  it  is  natural  to  assume  that  peak 
rate  R  =  //0  +  k  a0,  or  R  =  2 fiQ. 

The  function  on  the  right-hand  side  of  (5.4)  is  concave  but  with  very  low  cur¬ 
vature.  We  are  going  to  implement  a  conservative  over-dimensioning  approach 
by  approximating  the  concave  function  by  its  tangent  line.  It  is  easy  to  see  that 
the  linear  approximation  using  the  line  tangent  at  n'  is  the  following: 


C{n)  «  -z£a0 


(R o 


z£cr  o 

2 \fn* 


) 71 


or  introducing  notation: 


(5.5) 


r'  —  Ho  + 


Ze^0_ 

2y/rH 


we  obtain  the  following  fixed  charge  linear  approximation: 

|  0  if  n  =  0, 

C(n)  ~  <  (5.6) 

I  R'  +  r'  n  if  n  >  1 

We  will  call  R'  the  “fixed  charge”;  r'  will  be  called  the  “slope  coefficient”. 

For  example,  if  n'  =  20  and  k  =  4, 


R'  =  1.17  R 
r'  =  1.23r 


Here,  we  assume  that  for  a  typical  normal  approximation  R  =  '2/i.  is  the  peak 
rate,  r  =  /i  is  the  mean  rate  where  ji  is  the  mean  of  the  normal  distribution. 

This  approximation  is  plotted  in  Figure  5.2.  In  this  figure  ’o’  marks  corre¬ 
spond  to  quantiles  of  n-fold  convolution  of  normal  distribution  which  are  taken 
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at  the  highest  permitted  CLR  value  (s  =  10'  -8);’x’  marks  show  required  band¬ 
width  per  call;  solid  line  plots  formula  (5.4);  dashed  line  plots  approximation 
(5.6)  with  n!  =  20. 

Figure  5.3  plots  relative  approximation  error  (expressed  in  per  cent)  for  dif¬ 
ferent  choices  of  tangent  point  n'  =  3,  10,  20,  50.  When  implementing  this 
approach,  the  actual  choice  of  the  tangent  point  would  be  chosen  based  on  the 
number  of  video  calls  that  are  typically  multiplexed  together  on  a  link  and/or 
the  required  precision. 

Our  approximation  method  can  be  easily  generalized  for  the  case  when  there 
is  a  mix  of  independent  video  calls  h  €  R  having  different  distribution  types.  In 
this  case  expression  (5.3)  becomes: 

C  =  Vh  nh  +  ,./E^  (5.7) 

hen  Y  hen 

where  rih  is  the  number  of  calls  of  type  h.  Next,  we  make  another  approximation 
by  neglecting  terms  with  smaller  standard  deviation  values  in  the  summation 
under  square  root  in  (5.7).  This  gives  us  the  following  linear  approximation 
similar  to  (5.6): 


C 


(5.8) 


0  if  E hen  nh  =  0, 
maxftgi  R'h  +  E heL  r'h  nh  if  E hen  nh  >  1 
where  R'h  and  r'h  are  parameters  of  call  type  h  approximation. 

Using  alternative  notation,  let  K,p  be  the  set  of  calls  multiplexed  on  the  link 
(or  VP)  p  then: 


C(JCP) 


0  if  Kp  =  0, 
ilia R't  +  EkeK.  r'k  Otherwise. 


(5.9) 


We  call  the  value  ma xke/cpR'k  in  (5.9)  the  “variable  fixed”  cost  for  its  analogy 
to  the  fixed  cost  term  and  its  variable  nature. 
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This  approximation  is  illustrated  in  Figure  5.4.  In  this  example  every  forth 
call  has  mean  and  standard  deviation  3  times  higher  than  the  rest  of  the  calls. 
We  have  the  same  notation  as  in  Figure  5.2  for  quantiles  of  n-fold  convolutions 
and  per  call  bandwidths.  Solid  line  plots  formula  (5.7);  dotted  line  plots  the  same 
formula  but  with  averaged  mean  (1.5 pm  =  |(3 /i  +  (3 p))n)  and  standard  devi¬ 
ation  (\/3 an  =  \J\( 3ct2  +  (3a)2 )n)  where  p  and  a  are  parameters  of  the  lower 
bandwidth  calls;  dashed  line  plots  linear  approximation  (5.8)  (or,  its  equivalent 
(  5.9)). 

5.1.2  Data  Traffic  Capacity  Allocation 

In  this  section  we  derive  a  piece- wise  linear  approximation  for  the  capacity  allo¬ 
cation  rule  in  case  of  data  traffic.  We  assume  that  the  distribution  of  the  data 
transmission  rate  is  discrete,  with  only  two  non-zero  probabilities:  at  zero  and 
peak  traffic  values.  If  several  non-correlated  sources  are  multiplexed  together 
on  a  link,  their  traffic  distribution  is  calculated  as  convolution  of  their  original 
distributions.  The  minimal  required  bandwidth  is  determined  as  the  quantile  of 
the  negative  cumulative  density  function  (c.d.f.)  taken  at  maximal  permitted 
CLR. 

Let  p  ( q )  be  the  probability  the  source  is  transmitting  (silent).  Since  the 
source  is  considered  as  “on/off”,  p  +  q  =  1.  Suppose  all  sources  are  i.i.d.  If 
statistical  multiplexing  is  not  employed,  the  bandwidth  required  to  transmit 
data  from  n  sources  with  peak  rate  R  is  simply  nR.  If  statistical  multiplexing 
is  used,  the  required  bandwidth  will  be  (n  —  g)R.  We  will  call  g  statistical 
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multiplexing  gain  for  data  transmission.  Clearly, 

g{n)  =  sup{#  :  ^  pn  kqk  <  s}  (5.10) 

k=0  V  ^  ) 

where  £  is  the  highest  permitted  CLR. 

It  is  easy  to  see  that  if  n  is  small  enough  so  that  pn  >  e,  there  is  no  statistical 
multiplexing  gain,  i.e.  g(n)  =  0.  In  this  case  the  minimal  required  bandwidth  is 
equal  to  the  total  peak  rate:  bo(n)  =  Rn.  On  the  other  hand,  as  n  gets  larger, 
the  binomial  distribution  in  (5.10)  can  be  approximated  by  normal  with  mean 
np  and  standard  deviation  y/npq.  Then,  following  procedure  in  section  5.1.1,  we 
approximate  minimal  required  bandwidth  as: 

C(n)  npR  +  zeRy/npq  (5.11) 


As  we  did  in  section  5.1.1,  we  further  approximate  5.11  with  a  linear  function 
bo  ~  R  +  r'  n  where  R  =  [lne/ln p]qR-  Thus  the  approximating  function  is 
piece-wise  linear  and  consists  of  two  pieces: 


Rn  if  n  <  [lnt/lnp], 
R  +  r'  n  otherwise. 


(5,12) 


Notice,  that  the  first  equation  in  (5.12)  is  exact.  This  approximation  is  illustrated 
in  Figure  5.5.  Here,  we  consider  negative  c.d.f.  for  the  n-fold  convolutions 
of  discrete  distribution  with  Prob{transmission  rate  =  0 Mb/s}  =  0.9  and 
Prob{transmission  rate  =  1  Mb/s}  =  0.1  where  n  =  1, . . . ,  64.  Quantiles, 
taken  at  the  highest  permitted  CLR  value  (a  =  10— 8 ) ,  are  plotted  in  Figure  5.5 
where  the  ’o’  marks  correspond  to  the  total  required  bandwidth,  ’ x’  marks  show 
required  bandwidth  per  call,  solid  line  plots  approximation  (5.11),  and  dashed 
line  plots  approximation  (5.12).  As  it  is  seen  in  the  figure,  the  approximated 
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function  for  the  minimal  required  bandwidth  is  step-wise.  This  is  explained  by 
the  fact  that  the  minimal  required  bandwidth  is  a  multiple  of  peak  rate. 

Now  consider  statistical  multiplexing  of  data  calls  of  different  types  h  G  % 
with  different  “on/off”  probabilities  ph  +  (p,  =  1  and  different  peak  rates  R/,. 
Then  following  our  models  for  i.i.d.  data  calls  and  video  calls  of  different  types, 
we  obtain  that: 


C 


E hen  Rh  nh  if  -  E hen  nh  In  Ph  <  ~  In  fir, 
maxft^  Rh  +  E hen  A  nh  otherwise. 

Notice,  that  the  first  equation  in  (5.13)  is  exact. 


(5.13) 


5.1.3  Mixing  Heterogeneous  Calls 


The  following  model  is  proposed  for  the  case  when  there  are  calls  of  different 
types  multiplexed  together.  Let  'HVideo  and  'Hdata  be  the  sets  of  different  video 
and  data  call  types,  respectively;  R  =  'Hvideo  data-  Then 


a, 


C(tCVideo)  +  T.hendata  Rh'nh  if  E/ig'H  nh  lnp/(,  T 


[  max(maxfte%fe  Rh.  maxhenvideo  R’h)  +  Ehen  A  nh  otherwise, 
where  C'(/C,!ulifi0)  is  obtained  from  (5.9). 

Obviously,  the  second  part  of  the  formula  (5.14)  is  a  slight  modification 
of  (5.8)  where  the  mix  of  data  calls  is  approximated  as  normally  distributed. 
Concerning  the  first  part,  it  is  easy  to  see  that  this  implements  a  conservative 
approach  by  ignoring  statistical  multiplexing  between  data  and  video  calls.  How¬ 
ever,  since  in  this  case  both  bounds  are  rather  tight,  the  loss  of  accuracy  is  not 
very  high.  In  Figure  5.6,  we  plot  the  minimal  bandwidth  to  meet  CLR  require¬ 
ments.  In  this  example  we  consider  several  data  calls  multiplexed  together  with 
one  video  call.  Data  calls  are  i.i.d.  having  the  same  parameters  as  in  the  example 
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Minimal  Required  Link  Capacity: 


^^xxxxxxxxxxxx 


><XX><XXXX^^ 


10  20  30  40  50 

Number  of  Calls  Multiplexed  together  on  the 


pictured  in  Figure  5.5,  section  5.1.2  ( R  =  1,  p  =  0.1)  and  the  video  call  has  peak 
rate  10  and  mean  rate  5.  The  numbers  on  the  horizontal  axis  show  the  number 
of  data  calls  and  the  dashed  lines  are  a  linear  approximation  (5.13).  As  can  be 
seen  from  the  plot,  relative  error  of  the  approximation  does  not  exceed  20%. 

5.1.4  A  Comparison  with  Related  Results 

Gaussian  Self-Similar  Traffic 

The  network  traffic  stochastic  time  dependence  is  widely  regarded  to  be  statisti¬ 
cally  self-similar  meaning  that  it  shows  fractal  behavior  over  various  time  scales. 
This  concept  is  based  on  traffic  measures  taken  at  Bellcore  [27,  26].  It  implies 
the  following  capacity  allocation  rule  [8]: 

C  =  m  +  <I>  '(c)n'  s2//%  "’/"m1"2"’  (5.15) 

where 

m  is  total  mean  transmission  rate; 
a  is  transmission  rate  variance; 
x  is  buffer  size; 

H  is  Hurst  parameter,  1/2  <  H  <  1. 

It  can  be  seen  that  expression  for  C  in  (5.15)  consists  of  two  terms:  linear 
and  sublinear  on  mean  transmission  rate  m.  The  latter  is  dominated  by  the  first 
term.  Thus,  applying  steps  outlined  in  section  5.1.1,  we  can  approximate  C  as 
a  function  of  m  and  derive  a  fixed  charge  linear  approximation  similar  to  (5.4). 
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Two-Piece  Linear  Approximation 


In  references  [28,  29],  the  authors  showed  that  the  minimum  capacity  required 
for  a  number  of  “on-off”  VBR  connections  multiplexed  together  can  be  approx¬ 
imated  as  a  two-piece  linear  function  similar  to  (5.12).  In  [28],  the  authors 
proposed  a  simple  analytical  model  which  results  in  a  two-piece  linear  approx¬ 
imation.  In  [29],  in  addition  to  two  linear  pieces,  the  authors  considered  a 
“transient”  curve  between  the  two  approximating  lines. 

In  both  of  these  approaches,  the  first  linear  piece  of  the  approximation  corre¬ 
sponded  to  the  total  peak  rate  capacity  allocation  for  all  multiplexed  calls;  the 
second  piece  had  the  slope  approximately  equal  to  the  mean  rate  of  the  calls. 
Thus,  this  approximation  is  similar  to  (5.12)  derived  in  section  5.1.2. 

5.2  Problem  Formulation 

The  problem  formulated  in  this  section  can  be  described  as  follows.  Suppose 
there  is  a  certain  number  of  various  call  connections  to  be  set  up  between  each 
node  pair  in  a  network.  Each  call  is  characterized  by  its  mean  and  peak  band- 
widths  and  its  largest  tolerable  delay.  In  order  to  simplify  routing  and  reduce 
call  delay  and  network  control  overhead,  we  wish  to  establish  multilink  VPs  be¬ 
tween  some  pairs  of  nodes.  These  VPs  are  considered  as  logical  macro-links. 
The  problem  is  to  choose  which  node  pairs  are  to  be  connected  by  VPs,  what 
route  each  VP  should  take  and  the  amount  of  capacity  to  be  allocated  to  each 
VP. 

The  problem  constraints  are  dictated  by  the  following  two  QOS  requirements: 

1.  maximum  call  delay; 
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2.  minimum  cell  loss  probability. 


The  essential  property  of  cell  switching  in  ATM  networks  is  that  queuing  delay 
is  practically  completely  eliminated  at  the  intermediate  switching  nodes  of  a  VP 
through  which  the  cell  is  routed  [25,  35,  41].  Thus,  in  our  model,  we  assume  an 
end-to-end  delay  proportional  to  the  number  of  traversed  VPs. 

A  cell  loss  probability  constraint  translates  into  a  capacity  allocation  rule 
which  provides  a  guaranteed  cell  loss  probability.  As  was  shown  earlier,  the 
capacity  allocation  rule  is  essentially  non-linear.  For  the  purposes  of  our  model 
we  will  be  using  its  piece- wise  linear  approximations  derived  in  section  5.1. 

5.2.1  Bi-Criteria  Objective  Function 

A  good  VP  layout  involves  a  fundamental  tradeoff  between  the  network  through¬ 
put  and  its  grade  of  service.  Consider  the  quality  of  service  parameter,  transmis¬ 
sion  delay.  It  is  the  sum  of  propagation  and  queuing/switching  delays.  The  first 
part  is  negligibly  small  in  an  ATM  network  unless  it  contains  satellite  links.  The 
second  part  is  determined  by  the  number  of  VP-to-VP  switching  points  since, 
as  was  assumed  earlier,  delay  at  the  intermediate  VP  nodes  is  relatively  small. 
Hence,  transmission  delay  is  substantially  decreased  by  reducing  the  number  of 
VP-to-VP  switching  nodes. 

Thus,  establishing  “long”  VPs  decreases  transmission  delay  and  network  con¬ 
trol  overhead  and  simplifies  routing.  On  the  other  hand,  establishing  too  many 
low-capacity  VPs  may  take  away  statistical  multiplexing  gain  obtained  by  shar¬ 
ing  VP  capacity  among  multiple  variable-bit  rate  (VBR)  calls. 

However,  using  too  many  “short”  VPs  in  order  to  take  advantage  of  statistical 
multiplexing  generally  results  in  an  unacceptably  high  call  delay  and  network 
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Short, 

Long, 

High-capacity  VPs 

Low-capacity  VPs 

Throughput 

larger 

smaller 

Control  Overhead 

larger 

smaller 

Delay 

larger 

smaller 

Table  5.1:  VP  Layout:  Fundamental  Tradeoff 


control  overhead. 

This  concept  is  summarized  in  the  Table  5.1. 

In  order  to  capture  the  VP  layout  tradeoff,  we  formulate  the  problem  as 
bi-criteria  optimization: 

•  maximize  the  minimum  over  all  Origin-Destination  (O-D)  pairs  and  traf¬ 
fic  classes  of  the  ratio  of  admitted  traffic  to  traffic  demanded; 

•  minimize  the  total  mean  traffic  flow  through  the  VP  switching  nodes 
(except  at  origin  and  destination). 

The  ratio  in  the  first  part  of  the  objective  function  will  be  called  “relative 
throughput” .  Maximizing  the  minimum  of  the  relative  throughput  corresponds 
to  the  concept  of  “fairness”  in  the  ATM  network  management.  In  our  model, 
relative  throughput  may  take  values  greater  than  one  which  means  that  network 
can  accommodate  higher  throughput  than  projected  demand. 

The  second  part  of  the  objective  function  represents  the  amount  of  traffic 
flow  going  through  VP-to-VP  switches.  Since  these  nodes  are  intermediate  for 
the  virtual  channels,  this  part  will  be  called  “flow  through  intermediate  nodes” 
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or  “intermediate  flow”,  for  short.  Based  on  our  assumption,  this  flow  amount  is 
a  measure  of  overall  transmission  delay  and  control  overhead.  Note  that  if  there 
was  a  direct  VP  between  every  O-D  pair  then  this  value  would  be  zero. 

The  solution  to  this  type  of  optimization  problem  is  a  set  of  individual  so¬ 
lutions  to  optimization  problems  with  objective  function  composed  of  the  two 
criteria  taken  with  certain  weights.  The  whole  set  (also  known  as  Pareto-optimal 
solution  set)  is  obtained  by  varying  these  weights.  It  is  convenient  to  represent 
the  solution  set  as  points  on  the  relative  throughput  -  intermediate  flow  plane. 
It  will  be  also  called  a  tradeoff  curve. 

5.2.2  Basic  VP  Layout  Formulation 

In  this  section,  we  give  a  formulation  of  the  bi-criteria  optimal  VP  layout  problem 
which  we  consider  to  be  basic  in  the  sense  that  we  do  not  include  fault  tolerant 
constraints  into  this  formulation. 

Let  us  represent  a  telecommunication  network  as  a  graph  Q  =  (V,  A)  where 
J\f  and  A  are  the  sets  of  nodes  and  one-directional  links,  respectively.  Denote 
the  set  of  all  possible  routes  between  nodes  i  and  j  by  Vij  and  let  V  =  (J  Vij 
be  the  set  of  all  possible  virtual  paths  in  the  network.  Introduce  the  following 
notation: 

Vp  ,  capacity  assigned  to  VP  p  (unknown  variable); 

Be  ,  capacity  of  physical  link  e  (fixed). 

{1  if  VP  p  passes  through  the  link  e, 

0  otherwise; 

OUT(i)  ( IN(i )  )  ,  set  of  candidate  VPs  beginning  (ending)  at  node  i  G  J\f. 
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Let  1C  be  the  set  of  calls  to  be  (re)routed.  This  call  set  can  be  viewed  as  a 
known  network  traffic  pattern  or  a  set  of  calls  which  exists  in  the  network  at  a 
particular  time,  i.e.  a  “snapshot”.  With  each  call  A:  we  have  associated  its  origin 
and  destination  nodes  0(k )  and  D(k ),  respectively,  its  maximum  tolerable  delay 
dk  and  mean  and  peak  bandwidths  rk  and  Rk,  respectively.  We  introduce  two 
groups  of  integer  zero-one  call  variables: 

{1  if  call  A:  uses  VP  p. 

0  otherwise; 

{1  if  call  A:  is  accepted, 

0  if  it  is  blocked. 

Thus,  any  subset  of  “core”  call  set  ICcore  =  {k  G  /C  :  yk  =  1}  is  guaranteed  to  be 
routed. 

We  assume  that  call  set  1C  consists  of  two  subsets:  JCmdeo  and  ]Cdcda .  containing 
video  and  data  calls,  respectively.  Since  voice  calls  are  CBR,  for  the  purposes  of 
our  modeling,  we  consider  them  as  a  special  case  of  video  calls  with  zero  fixed 
charge.  In  general,  we  consider  a  set  R.  of  all  call  classes. 

Let  A';';  be  the  set  of  class  h  e  K  calls  with  the  origin  i  and  destination  j  and 
rk  be  the  mean  transmission  rate  of  call  k.  Then  the  ratio  of  the  total  traffic- 
routed  between  nodes  i  and  j  to  the  total  traffic  demand  between  these  nodes  is 

5Z  rkVk/  5Z  rk  ■ 

k-  V  kcK% 

We  would  like  to  maximize,  (>.  the  lowest  of  these  ratios  over  all  traffic  classes 
and  node  pairs  and  at  the  same  time  to  minimize  the  traffic  flow  through  the 
VP-to-VP  switches.  This  gives  us  the  following  expression  for  the  objective 
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function: 


Maximize 
subject  to 


P~vYrk  Y  XkP 

keK.  Pev-ouT(o(k)) 


(5.16) 


Y  rkVk  >  p  Y  rk  for  a11  '>■  j  eAf,hen  (5.17) 

keic^  *e/c*. 

Here,  coefficient  r)  is  used  to  put  certain  weights  on  the  two  criteria  which  are 
optimized  at  the  same  time. 

As  was  shown  in  section  5.1.2,  the  amount  of  capacity  required  to  carry  data 
traffic  is  a  concave  function  of  the  number  of  calls  carried.  We  approximate  this 
function  by  a  two-piece  linear  function.  The  first  piece  corresponds  to  the  case 
when  no  statistical  multiplexing  gain  is  observed  and  the  bandwidth  allocated  per 
call  is  exactly  its  peak  rate.  The  second  linear  piece  corresponds  to  the  case  when 
statistical  multiplexing  gain  is  in  effect.  We  will  call  these  cases  “non-saturated” 
and  “saturated”,  respectively.  Since  the  linear  approximation  is  concave,  we 
need  one  additional  zero-one  variable  per  link  and  two  continuous  variables  per 
call  and  link  pair  to  provide  an  adequate  linear  programming  formulation. 


1  if  the  data  traffic  is  saturated  on  VP  p: 


sp  ~  j 

0  otherwise; 

C/sp  is  the  slope  of  the  “non-saturated”  linear  piece,  when  traffic  is  not  saturated, 
zero,  otherwise. 


£kp  is  the  slope  of  the  “saturated”  linear  piece,  when  traffic  is  saturated,  zero, 
otherwise. 


The  next  group  of  constraints  forces  sp  to  one  when  traffic  is  saturated,  that  is, 
according  to  (5.14),  when  Y.keicdaia~^n  Pkxkp  A  —  Ins,  otherwise  sp  is  forced  to 
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zero. 


(1  -  sp)  (- lnt)  <  J2  -  In  PkXkp  ~  (-  Int)  < 

keKdata 

sp  ^2  ~  ^Pk  for  all  p  E  V,  (5.18) 

keKdata 

Coefficient  (—  lne)  in  the  left  hand  side  of  (5.18)  is  chosen  as  a  lower  bound  on 
the  expression  in  the  middle  part  of  (5.18);  coefficient  ^2  ~^nPk  in  the  right 

keKdaia 

hand  side  part  is  chosen  as  an  upper  bound  on  the  expression  in  the  middle  part. 

The  next  three  groups  of  constraints  define  a  fixed  charge,  Adata,  and  slopes 
for  two-piece  linear  approximation  to  capacity  allocation  for  data  traffic. 


Adata  >  Rki^kp  —  (1  —  sp))  for  all  p  EV,  k  E  JCdata,  (5.19) 

Rk  {'■I'kp  ~  sp)  <  (kp  <  (1  -  sp)Rk  for  all  p  E  V.  k  E  1Cdata ,  (5.20) 

rk{xkP  ~  (1  -  sp))  <  E,kp  <  rksp  for  all  p  E  V,  k  E  JCdata,  (5.21) 

0  <  sp  <  1  and  integer  for  all  p  E  V ,  (5.22) 

CkP,  tkp  >  0  for  all  pEV,  k  E  1C,  (5.23) 


The  following  two  groups  of  constraints  define  the  total  fixed  charge  as  the 
maximum  of  fixed  charges  over  all  call  classes. 


Ap  V  RkXkp 

for  all  p  E  V,  k  E  JCmde0, 

(5.24) 

A  ">  A  dodo 
d-*p  E  t-ip 

for  all  p  E  V , 

(5.25) 

The  next  two  constraint  groups  define  capacity  allocated  to  VP  p 

and  place 

capacity  constraints  on  physical  link  e,  respectively. 

Vp  A  Ap  + 

'y  '  1  k'Xkp  T  y  '  (C kp  +  Ckp) 

for  all  p  E  V, 

(5.26) 

kefCvideo  ketcdata 

SpeVP  <  Be 

for  all  e  E  A, 

(5.27) 

pev 

vp  >  0 

for  all  p  EV . 

(5.28) 
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Flow  conservation  constraints  introduced  next  and  integrality  constraints 
provide  a  single  path  between  origin  and  destination  of  the  accepted  calls. 


~Vk  if  D{k)  =  i 

^  '  %kp 

^  '  %kp  —  * 

yk  if  O(k)  =  i 

for  all  i  G  A",  k  G  /C, 

(5.29) 

peouT(i) 

pelN(i) 

0  otherwise 

VI 

VI 

o 

and  integer 

for  all  k  G  V. 

(5.30) 

Xkp  >  0  and  integer 

for  all  p  G  V,  k  G  1C, 

(5.31) 

Finally,  the  basic  formulation  of  the  VP  layout  problem  is  given  by 
BAS:  Maximize  (5.16)  subject  to  (5.17-5.31) 

It  can  be  seen  from  the  problem  formulation  that  the  best  way  to  reduce  the 
intermediate  flow  is  to  establish  VPs  between  each  pair  of  nodes.  However,  if 
this  is  done,  then  we  might  not  get  enough  statistical  multiplexing  gain  and,  as 
a  result,  would  not  have  enough  capacity  to  route  all  calls.  On  the  other  hand, 
establishing  VPs  using  exactly  one  physical  link  maximizes  statistical  multiplex¬ 
ing  gain  but  may  make  intermediate  flow  rather  high.  Short  VPs  also  tend  to 
increase  the  network  control  overhead.  In  order  to  favor  longer  VPs,  objective 
function  coefficient  i)  should  be  made  larger. 

Below,  we  study  the  properties  of  the  formulation  by  considering  two  marginal 
cases. 

Suppose,  we  perform  a  pure  relative  throughput  maximization  ignoring  the 
second  part  of  the  objective  function  (5.16),  or  equivalently,  coefficient  p  is  set 
to  zero.  Then,  it  is  easy  to  show  that  in  the  optimal  solution  each  VP  will 
traverse  one  physical  link  and  one  physical  link  will  contain  at  most  one  VP. 
Indeed,  if  there  are  VPs  traversing  more  than  one  physical  link,  then  they  can 
be  partitioned  into  sequences  of  VPs  associated  with  the  traversed  physical  links 
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and  the  objective  function  will  not  be  affected.  After  that,  any  existing  VP  will 
share  exactly  one  physical  link  with  some  other  VPs.  If,  now,  all  VPs  sharing  the 
same  link  are  merged  together  it  will  loosen  constraint  (5.27)  and  the  solution 
will  stay  feasible.  We  can  formulate  this  statement  as  the  following 

Proposition  1  If  in  the  problem  BAS  g  =  0,  then  there  is  an  optimal  solution 
having  VPs  equivalent  to  physical  links. 

Now,  assume  there  exists  a  VP  layout  connecting  each  O-D  pair  by  a  set  of 
VPs  and  admitting  at  least  one  call  in  each  class  between  this  pair.  Clearly,  with 
this  layout  the  second  part  of  the  objective  function  becomes  zero  and  solution 
value  becomes  positive.  On  the  other  hand,  for  any  feasible  solution  with  non¬ 
zero  second  criterion  (i.e.  non-zero  intermediate  flow),  coefficient  g  can  be  chosen 
large  enough  so  that  the  solution  value  be  negative.  If  such  end-to-end  layout 
does  not  exist,  then  zero  solution  will  still  be  feasible.  This  argument  implies 
the  following 

Proposition  2  If  coefficient  rj  in  objective  function  (5.16)  is  chosen  sufficiently 
large,  for  any  optimal  solution  the  second  term  of  the  objective  function  will  equal 
zero. 

5.2.3  Fault  Tolerant  Routing  Constraints 

In  this  section  we  expand  our  model  to  capture  link  failure  case.  In  this  case  we 
wish  to  reroute  calls  as  fast  as  possible  taking  advantage  of  the  two  level  VC/VP 
ATM  network  architecture.  That  is,  together  with  primary  VPs,  we  want  to 
establish  backup  VPs  so  that  in  the  event  of  failure  only  VP  level  rerouting  is 
required.  Since  this  kind  of  rerouting  does  not  touch  the  VC  level,  particularly 
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VC  identifier  tables,  it  is  expected  to  work  very  fast,  changing  only  VP  identifier 
tables.  We  propose  two  different  strategies: 

1.  Backup  VPs  use  paths  which  are  link-disjoint  from  the  primary  VPs.  When 
a  link  fails,  all  affected  primary  VPs  are  replaced  by  their  backups. 

2.  A  backup  VP  is  designed  for  each  VP  link  pair,  so  that  potentially,  there 
is  a  different  replacement  for  each  link  on  the  primary  VP. 

We  start  with  an  assumption  that  between  each  node  pair  there  are  at  most 
one  primary  and  one  backup  VP  established.  If  there  is  a  need  to  establish  more 
than  one  VP  between  the  same  end-points  (e.g.,  to  accommodate  multiple  traffic 
types),  then  our  problem  formulation  can  be  generalized  quite  easily. 


Link-Disjoint  Backup  Paths 

In  addition  to  the  variables  and  notation  introduced  earlier  in  this  chapter,  we 
define  the  following  zero-one  and  continuous  decision  variables  for  backup  VPs. 

{1  if  route  p  E  Vij  is  used  as  a  primary  VP, 

0  otherwise; 

{1  if  route  r  E  V,  ,  is  used  as  a  backup  VP, 

0  otherwise; 

Ur  =  capacity  assigned  to  backup  VP  r. 

If  vp  (ur)  is  equal  to  zero,  then  there  is  no  capacity  assigned  to  the  primary 
(backup)  VP,  which  means  that  Vp  =  0  (Ur  =  0).  If  any  of  the  variables  vp  or  ur 
is  one,  then  capacity  of  the  corresponding  VP  is  non-zero  but  does  not  exceed 
the  smallest  link  capacity  on  the  VP  route.  This  relation  can  be  formulated  as 
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follows: 


Vv  <  vpmm  Be  for  all  p  E  V  (5.32) 

y  '  e&p 

Ur  <  ur  min  Be  for  all  r  E  V  (5.33) 

eer 

Since  there  is  at  most  one  primary  and  backup  VP  between  each  pair  of  nodes, 
E  vv  =  E  «r  <  1  for  all  Vij.  (5.34) 

PC  T;  ,  r  Vij 

Since  the  backup  VP  should  have  the  same  assigned  capacity, 

E  Vp  =  E  Ur  for  all  Vij.  (5.35) 

peVij  reVij 

Finally,  the  link-disjoint  condition  is: 

E  Vp^pe  +  E  ur^re  <  1  for  all  Vij,  e  E  A  (5.36) 

peVij  rev^ 

Before  we  formulate  the  link  capacity  constraint,  let  us  introduce  one  more  group 
of  variables: 

Zpr,  the  amount  of  flow  transferred  from  the  primary  VP  p  to  the  backup  VP  r, 
when  the  primary  VP  is  affected  by  the  link  failure. 

In  our  model,  we  assume  all  flow  on  the  primary  VP  is  transferred  to  the  backup 
VP  so  that  Zpr  =  Vpur  =  vpUr  since  only  in  case  of  vp  =  ur  =  1  is  capacity 
transferred  from  a  primary  VP  p  to  its  backup  r.  However,  expressing  Zpr  in 
this  way  makes  the  problem  non-linear  which  is  not  desirable.  Instead,  the 
following  group  of  linear  constraints  gives  the  equivalent  expression: 

for  all  p  E  Vij  and  all  Vij,  (5.37) 

for  all  r  E  Vij  and  all  Vij,  (5.38) 


E  Zpr  —  Vp 


r(  Pi 


E  Zpr  ~  Ur 

P&'Pij 
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Suppose  link  e  fails.  Consider  the  change  of  flow  through  an  arbitrary  link 

e'  ^  e.  The  initial  flow  through  that  link  was  ^  8pe'Vp.  Since  link  e  fails,  all 

per 

VPs  p  satisfying  Spe  =  1  are  taken  down  and  the  total  flow  decrease  on  link  e! 

becomes  ^  dpe5pe>\rp.  The  increase  in  the  flow  on  link  e'  caused  by  new  backup 

Pev 

VPs  is  E  E  8re>5peZpr.  Finally,  capacity  constraints  for  the  links  become: 

all  Vij  p-rev-ij 

E(!  -  +  E  E  Sre'SpgZpr  <  Bei  for  all  e'  ±  e  <G  45-39) 

p  T  all  Vij  P'r<=vij 

The  problem  for  the  link-disjoint  path  is  given  by  (5.16-5.39)  and  zero-one  con¬ 
straints  for  vp  and  ur  as  well  as  non-negativity  constraints  for  Ur  and  Zpr: 


0  <  vp 

<  1 

and  integer 

for  all  p  G  V. 

(5.40) 

0  <  ur 

<  1 

and  integer 

for  all  r  <E  V. 

(5.41) 

Ur  >  0 

for  all  r  G  V. 

(5.42) 

Zpr  >  o 

for  all  p ,  r  G  V. 

(5.43) 

Link-Dependent  Backup  Paths 

In  order  to  implement  the  second  strategy,  he.  link-dependent  rerouting,  we 
treat  backup  route  r  as  paired  with  a  potentially  failed  link  e.  Then,  variables 
ur  and  Ur  are  replaced  with  link-specific  variables  ure  and  Ure,  respectively,  which 
correspond  to  the  routes  chosen  as  backup  VPs  and  activated  in  the  case  of  link 
e  failure. 

The  problem  formulation  for  the  link  dependent  rerouting  is  similar  to  link 
disjoint  formulation.  Here  constraints  (5.36)  forcing  the  link-disjoint  paths  be¬ 
come 


oU'pp  — 


0  for  all  r  E  V,  all  e  G  A, 


(5.44) 


which  bar  backup  VPs  from  passing  through  the  link  e. 

If  primary  VP  does  not  pass  through  the  link  e,  then  the  corresponding 
backup  VP  ure  is  not  setup.  This  requirement,  together  with  a  condition  limiting 
the  number  of  VPs  with  the  same  end  points  to  one,  replaces  constraints  (5.34- 
5.35)  as  follows: 


E 

Vp<  1 

for 

all 

Vij 

(5.45) 

p&'p-ij 

E 

ure  - 

—  ^  '  fipe^p 

for 

all 

Piji 

all 

e  G 

A 

(5.46) 

reVij 

P&'Pij 

E 

Ure  - 

=  E  Vp 

for 

all 

A?  5 

all 

e  G 

A 

(5.47) 

reVij 

p&Vij 

Constraints  (5.46-5.47)  do  not  let  backup  VP  ure  be  set  up  if  the  primary  VP 
with  the  same  end  points  does  not  go  through  the  link  e.  Constraints  (5.38) 
become: 

Ure  E  AAe  for  all  r  G  Tjj,  all  V-,j,  and  all  e  G  A.  (5.48) 

P&'Pij 

Finally,  we  have  the  following  formulation  of  the  problem. 

(5.16) 

subject  to 

(5.17-5.33),  (5.44-5.47),  (5.37),  (5.48),  (5.39-5.43) 

The  solution  to  the  problem  is  a  set  of  link  specific  backup  VPs  which  are 
activated  when  the  corresponding  link  fails. 

5.2.4  VP  contraction 

In  the  previous  section  it  was  assumed  that  any  primary  VP  should  be  backed  up 
by  a  backup  VP  with  the  same  capacity  since  the  currently  routed  calls  should 
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not  be  affected  by  a  link  failure.  This  may  cause  a  substantial  reduction  in  the 
number  of  calls  routed  over  a  VP  in  order  to  guarantee  their  complete  tolerance 
to  the  link  failure.  Alternately,  if  it  is  possible  to  take  down  some  calls  in  case 
of  link  failure,  the  number  of  calls  routed  under  normal  conditions  may  be  much 
higher.  In  this  case  the  backup  VP  may  have  lower  assigned  capacity  than  the 
primary.  Notice  that  we  do  not  try  to  reroute  calls  to  different  paths  due  to 
large  potential  number  of  the  affected  calls  and,  in  many  cases,  unreasonable 
complications  in  the  network  management  protocol.  We  assume  that  it  is  quite 
easy  to  take  down  some  calls.  We  study  this  model  by  modifying  the  formulation 
presented  in  section  5.2.3. 


Link-Disjoint  Backup  Paths 

In  the  case  of  VP  contraction,  in  addition  to  call-to-VP  assignment  variables  xkp, 
we  must  introduce  variables  indicating  which  calls  (VCs)  are  rerouted  together 
with  a  VP  when  it  is  affected  by  a  failure. 

{1  if  call  k  uses  VP  p  and  is  not  dropped  when  failure  occurs 
0  otherwise; 

Variables  xkp  and  x'kp  define  two  VP  sets  Vk  and  Vk ,  respectively.  The  first  set 
is  the  set  of  VPs  used  for  call  k  routing,  the  second  is  the  set  of  VPs  where  call 
k  is  backed  up.  Clearly,  V'k  C  Vk.  or 


x'kp  <  xkp  for  all  k  e  1C  and  p  €  V  (5.49) 

Let  fp  be  the  probability  that  VP  p  fails  during  the  next  time  period.  Since  any 
two  VPs  in  the  set  Vk  use  different  links  (pi,  p->  G  Vk ,  PiC\p2  =  0)  and  links 
fail  independently,  failures  of  VPs  pi  and  p2  are  also  independent.  Since  call  k 
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is  taken  down  if  and  only  if  VP  p  £  Vk\P'k  fails,  then 


Prob{ call  k  is  taken  down}  =  1—  (1  —  fp)  = 

pevk \Vk 

i~  n  (!-/?)-  e  fp=j2fp(xkp-x'kP)- 

p-xkp-x'kp=i  p-xkp-x'kp= i  pep 

The  above  approximation  is  valid  when  we  neglect  the  multiple  failures  case  and 
assume  fp  is  small. 

The  expected  loss  of  traffic  flow  due  to  a  link  failure  then  becomes 

E  rk  E  fpixkp  -  Xkp)  (5.50) 

keic  pev 

Since  the  backup  VP  in  the  current  model  may  have  a  bandwidth  lower  than 
the  primary,  the  corresponding  constraints  are  to  be  modified.  Constraints  (5.35) 
become 

E  Vp>  E  Ur  for  all  Tij,  (5.51) 

peVij  reVij 

and  (5.37)  becomes 

E  Zpr  <  Vp  for  all  p  £  Vjj  and  all  Vij ,  (5.52) 

rev i, 

whereas  (5.38)  stays  the  same.  A  failed  primary  VP  is  replaced  by  a  lower 

capacity  VP  with  new  capacity  ^  Zpr  where  V(p)  is  a  group  of  routes  between 

reV(p) 

nodes  i  and  j  which  are  the  end  points  for  VP  p.  Hence  the  new  constraints 
similar  to  (5.26-5.24)  are  added  for  the  remaining  calls  with  guaranteed  rerouting. 

J2  zpr  >  \  +  E  rk'x'kP  for  a11  peP,  (5.53) 

rEV(p)  k£lC 

A p  >  Rkx'kp  for  all  p  £  V,  k  £  Af,  (5.54) 
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The  entire  problem  becomes: 


Minimize  P~vJ2  r*(  5Z  xkP  +  fp(xkP  ~  x'kP ))  (5-55) 

keK  pev-ouT(o(k ))  per 

subject  to 

(5.17-5.34),  (5.51),  (5.36),  (5.52),  (5.38-5.43),  (5.49),  (5.53) 
x'kp  >  0  and  integer  for  all  p  G  V  and  k  E  1C.  (5.56) 

In  this  problem,  the  objective  function  (5.55)  is  obtained  by  composing  the 
objective  function  (5.16)  and  expression  (5.50)  for  the  expected  traffic  loss. 


Link-Dependent  Backup  Paths 

In  order  to  formulate  the  problem  for  the  second  strategy,  i.e.  link  dependent 
rerouting,  it  has  to  be  assumed  that  lower  amount  of  the  capacity  transferred  to 
the  backup  VP  is  the  same  for  the  failure  of  any  link.  Otherwise  the  variables 
Zpr  would  depend  on  failed  link  e  and  the  problem  would  become  intractable. 

The  problem  formulation  can  be  obtained  by  replacing  constraints  (5.51) 
with 


E  U’  e  <  E  for  all  Vij,  all  e  €  A  (5.57) 

reVtj  pev-ij 

which  is  similar  to  (5.47).  Then  the  problem  formulation  becomes: 


(5.55) 

subject  to 

(5.17-5.33),  (5.44-5.46),  (5.57),  (5.52), 
(5.48),  (5.39-5.43),  (5.49),  (5.53),  (5.56) 
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5.2.5  Network  Dimensioning  Subproblem  for  the  Gener¬ 
alized  Fault-Tolerant  VP  Layout  Model 

In  this  section  we  discuss  a  dimensioning  subproblem  of  the  generalized  fault- 
tolerant  VP  layout  model  which  was  stated  in  the  introduction  to  this  chapter. 
Here  we  make  two  essential  modifications: 

•  the  set  of  VP  routes  between  a  pair  of  end  nodes  is  replaced  by  a  set  of 
link-disjoint  VP  pairs  with  the  same  end  nodes; 

•  if  the  decision  is  made  to  set  up  both  VPs  from  a  link-disjoint  pair,  then 
under  normal  conditions  call  traffic  is  admitted  to  both  VPs. 

The  first  modification  does  not  change  the  model.  Its  advantage  is  that  we  are 
able  to  omit  constraints  (5.36)  and  get  rid  of  zero-one  variables  associated  with 
VP  routes  which  dramatically  speeds  up  solution  process.  Although,  it  comes 
for  the  price  of  potential  increase  in  the  number  of  VP  variables. 

We  argue  that  the  second  modification  provides  a  better  solution  to  the  VP 
layout  problem  than  the  original  approach  when  a  backup  VP  has  zero  capacity. 
Indeed,  a  set  of  feasible  solutions  to  this  problem  includes  the  “zero  backup  VP” 
solution  as  a  special  case. 

On  the  other  hand,  the  second  modification  implies  that  in  case  of  a  VP 
failure,  VP  traffic  can  be  switched  to  the  other  VP  which  may  be  already  carrying- 
some  non-zero  traffic  of  its  own.  Thus,  the  network  protocol  should  be  able  to 
increase  VP  capacity  from  an  arbitrary  value.  This  procedure  may  require  more 
advanced  network  signaling  than  was  assumed  earlier  when  the  VP  capacity  had 
to  be  increased  from  a  zero  value. 
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It  is  not  practically  possible  to  consider  all  link-disjoint  pairs  of  VP  routes 
(without  resorting  to  a  column  generation  approach).  We  propose  an  approach 
where  a  limited  number  of  VP  route  pairs  are  selected  in  advance  using  spe¬ 
cial  routing  heuristics.  After  that  the  problem  becomes  a  type  of  dimensioning 
problem  which  allocates  capacity  to  the  VPs.  The  solution  obtained  is  in  general 
approximate  but,  as  will  be  shown  by  our  computational  experiments,  its  quality 
does  not  degrade  substantially. 

An  exact  formulation  of  this  model  will  be  presented  in  the  next  chapter.  The 
advantage  of  the  simplified  model  is  demonstrated  by  the  following  example. 

Consider  a  network  instance  shown  in  Figure  5.7.  Suppose  that  links  ei,  62.  e 3 
are  of  unit  capacity,  there  is  no  fixed  charge  for  the  traffic  between  node  pairs 
(A,  A'),  ( B,B '),  (C,  C"),  and  all  node  pairs  have  unit  mean  traffic  demand.  Let 
us  compare  the  two  strategies: 

•  If  traffic  between  each  node  pair  is  restricted  to  one  VP,  then  in  an  opti¬ 
mal  solution  each  link  carries  one  primary  VP  and  has  an  equal  capacity 
reserved  for  a  backup  VP.  Hence,  capacity  assigned  to  each  VP  is  equal  to 
1/2  and  the  total  traffic  routed  between  any  node  pair  is  1/2. 

•  If  traffic  between  each  node  pair  is  allowed  to  be  spread  between  both  pri¬ 
mary  and  backup  VPs,  then  there  is  a  solution  providing  better  through¬ 
put.  Indeed,  allocate  1/3  capacity  units  to  each  VP  and  reserve  1/3  capac¬ 
ity  capacity  units  on  each  link  for  backup  purposes.  Then,  it  is  easy  to  see 
that,  first,  link  capacities  are  not  exceeded  and,  second,  there  is  sufficient 
amount  of  capacity  to  reroute  traffic  affected  by  any  single  link  failure. 
The  total  traffic  between  any  node  in  this  case  is  2/3  which  is  higher  than 
1/2  in  the  previous  case. 
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Chapter  6 


Solution  Methods  for  the  VP  Layout 
Problem 


In  the  previous  chapter  we  gave  an  exact  integer  programming  formulation  of 
the  fault  tolerant  VP  layout  problem  and  considered  different  versions.  Theoret¬ 
ically,  if  this  formulation  is  used  as  input  into  a  general  IP  solver,  then  an  exact 
solution  can  be  obtained.  However,  because  of  the  complex  problem  structure, 
the  running  time  of  this  solution  process  would  be  prohibitively  long.  In  this 
section,  we  propose  an  approach  to  solving  the  problem  which  would  provide  us 
with  a  high  quality  approximate  solution. 

First,  let  us  outline  the  main  ideas. 

•  Linear  Relaxation.  In  order  to  speed  up  computations,  we  drop  the 
integrality  constraints  on  variables  %kp  and  y^. 

•  Aggregation.  Similar  to  the  approach  presented  in  Chapter  4,  we  suggest 
aggregating  ly,.  and  Xkp  variables  by  merging  the  same  class  calls  originating 
at  the  same  node  together.  This  approach  is  motivated  by  the  fact  that 
we  anticipate  a  large  number  of  calls  between  any  pair  of  nodes  and  we  are 
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more  interested  in  overall  demand  than  in  the  particular  call  routing. 


•  Heuristics.  We  obtain  integral  zero-one  values  for  variables  by  rounding- 
off  their  values  obtained  via  linear  programming  relaxation  using  various 
thresholds. 

•  Valid  Inequalities.  In  order  to  tighten  the  linear  relaxation  of  the  prob¬ 
lem,  we  introduce  several  groups  of  valid  inequalities.  This  enables  us  to 
improve  the  quality  of  our  approximate  solution. 


6.1  Aggregation  of  the  Call  Variables 


The  first  step  in  obtaining  an  approximate  to  the  solution  of  the  original  problem 
is  to  aggregate  call  variables.  Earlier,  we  defined  a  traffic  class  as  a  group  of  calls 
having  the  same  peak  and  mean  bit  rates.  Here  we  also  assume  that  video  and 
data  calls  belong  to  different  classes  since,  as  it  follows  from  our  earlier  discussion, 
they  require  different  statistical  multiplexing  models.  Similar  to  the  approach 
used  in  Chapter  4,  call  variables  are  aggregated  within  the  groups  corresponding 
to  common  originating  node  and  to  the  same  call  class.  Let  h(k)  denote  a  class 
call  k  belongs  to.  Then  the  new  aggregated  variables  are  defined  as  follows: 


Zip  =  E  rkZkp 

(6.1) 

k:  0(k)=i,  h(k)=h 

Vij  =  E  rkVk 

(6.2) 

k:  0(k)=i,D(k)=j ,  h(k)=h 

Next,  we  define  another  set  of  zero-one  variables  which  would  indicate  ad¬ 
mission  or  rejection  of  a  call  class  to  a  VP. 

[  1  if  call  class  l  is  admitted  to  VP  p 


Zph 


0  otherwise; 
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Obviously, 


E  *ip  —  miii{ili,-',  B(p)}zph  for  all  p  €  V,  h  <E  R  (6.3) 

ie  A' 

where  M-1  =  Y  r^.  is  the  maximal  amount  of  class  h  traffic  which 

k:  0(k)=i,  h(k)=h 

can  originate  at  node  i,  B(p )  =  mineGp  Be  is  minimal  capacity  of  a  link  traversed 
by  VP  p. 

Let  Rh  and  h  denote  peak  and  mean  bit  rate,  respectively,  of  calls  within 
class  h.  Also,  let  Rmde0  and  Rdata  be  subsets  of  R  consisting  of  video  and 
data  call  classes,  respectively.  As  was  mentioned  in  the  discussion  of  statistical 
multiplexing  mechanism  in  section  5.1,  data  traffic  is  approximated  by  a  two- 
piece  linear  function.  The  first  part  corresponds  to  what  is  called  the  “non- 
saturated”  case  and  it  expresses  direct  proportionality  dependence  on  the  number 
of  calls.  The  second  part  implies  a  fixed  charge  cost.  Similar  to  model  (BAS), 
we  introduce  zero-one  variable  sp  indicating  whether  the  data  traffic  on  VP  p 
is  saturated.  After  we  replace  Yxkp  in  (5.18)  by  Y.  yy  and  make  appropriate 
adjustments,  we  have: 

(1  -sp)LB{  J2  E(L“ln  Phlj1  ~  f-lnt])}  < 

E  E(L“ln  Ph\y-  -  r-lnel)  < 

ie'Hdata  *eA/”  h 

spUB{  Y  E(L-lnfVj^-  r-lnel)}  for  all  p  G  V,  (6.4) 

ieTidata  iejj  rh 

where  LB{  }  and  UB{  }  are  lower  and  upper  bounds,  respectively. 

If  we  also  use  two  continuous  variables  Q,h  and  £ph  for  each  data  class  h  and 
VP  p  to  describe  the  slopes  of  the  two-piece  linear  function,  then  the  expression 
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for  the  data  traffic  becomes: 


o  < 

C ph 

—  (1  sp )  t/t 

(6.5) 

o  < 

^  ip 

ieAT  7'h 

( ph 

A  spTPh 

(6.6) 

= 

Cph  T  fph 

(6.7) 

j r data 

vp 

= 

'y  ]  {RhCph  T  )  hfph) 

(6.8) 

hendata 


Here,  tp  =  refine/ In  p/j]  is  the  saturation  threshold  and  Tpp  is  the  highest  class 
h  traffic  level  on  link  p.  Notice  that  only  one  of  the  variables  (ph  and  £ ph  is 
non-zero  at  a  time. 

The  following  two  propositions  show  that  the  inequalities  chosen  above  are 
quite  “tight”  meaning  that  LP  relaxation  extreme  points  are  expected  to  be  close 
to  the  IP  solutions. 

Proposition  3  If  variable  x  has  negative  lower  and  positive  upper  bounds:  LB(x )  < 
0  <  UB(x),  then  the  two-dimensional  polyhedron  {{x,y)  :  (1  —  y)LB(x)  <  x  < 
yUB(x );  0  <  y  <  1}  has  integer  y-coordinates  of  its  extreme  points. 

Proof: 

As  can  be  seen  from  Figure  6.1(a),  the  polyhedron  has  the  following  ex¬ 
treme  points:  (LB(x),  0),  (0,1),  (UB(x),  1),  (0,0),  all  of  which  have  integer 
//-coordinates.  □ 

Corollary  1  Constraints  (6. 4)  together  with  0  <  sp  <  1  define  feasible  sets  with 
zero- one  sp  coordinates  of  their  extreme  points. 

Proposition  4  The  three-dimensional  polyhedron  defined  by  the  constraints  (6.5- 
6.6)  and  0  <  sp  <  1  in  the  space  (( ph,£ph,sp)  has  integer  sp- coordinates  of  its 
extreme  points. 
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Figure  6.1:  Polyhedra 


Proof: 

As  can  be  seen  from  Figure  6.1(b),  the  polyhedron  is  a  tetrahedron  with 
extreme  points:  (0,0,0),  (4,,  0,0),  (0,X/t,l),  (0,  0, 1),  all  of  which  have  integer 
Sp-coordinates.  □ 

Next,  we  introduce  an  expression  for  the  fixed  charge  Ap  corresponding  to 
VP  p: 

i  Rhzph  if  h  e  Rvideo 

Ap  >  <  for  all  h  €  R,  p  €  V  (6.9) 

Rh{zph  -  (1  -  sp))  if  he  Rdata 


Ap  >  0  for  all  p  €  V  (6.10) 

Suppose  now  the  traffic  classes  are  ordered  such  that, 

0  =  Ro  <  Ri  <  R2  <  ■  ■  ■  <  R\n\  <  1 

That  is,  we  arrange  all  (video  and  data)  traffic  classes  in  increasing  order  of  the 
fixed  charge  term  in  their  linear  approximation.  Class  0  with  no  fixed  charge  is 
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used  for  CBR  traffic.  Impose  the  following  “order”  constraints  on  ^-variables: 


1  >  zpi  >  zP2  >  ■  ■  ■  >  zp \u\  >  0  for  all  p  G  V  (6.11) 

Conceptually,  this  means  that  if  a  call  class  is  admitted  to  VP,  then  all  call 
classes  with  a  lower  fixed  charge  are  admitted  as  well. 

Now,  consider  an  expression  for  the  slope  of  the  linear  approximation  for 
video  traffic  function: 

v;ide°=  E  E  4  (6-i2) 

he'Hvideo  ie Af 

Then  the  entire  aggregated  approximation  of  the  problem  becomes: 

AGGR: 

Maximize  p  —  t)  E  E  E  4  (6-i3) 

ie M  pev-ouT(i)  hen 

subject  to 

Vij  >  P  E  rk  for  all  i  ^  j  G  N,  h  G  %  (6.14) 

keK'y 

~Vij  /  i 

E  £ip~  E  *rp  =  j  v-  rh  v  i,j  eNjien  (6.15) 

peOUT(j)  /;•  /  V ;  ?)  /  Vim  1  •/  1 


{.  meA /”—{*} 

(6.3-6.12) 

E  (^;ide°  +  Vpdata  +  Ap)<^  <  for  all  e  G  A.  (6.16) 
per 

zph  G  {0,  1}  for  all  p  G  V,  h  G  (6.17) 

Sp  G  {0,  1}  for  all  p  G  V  (6.18) 

p>  0  for  all  i,j  G  AT,  p  G  V.  h  G  H  (6.19) 


In  this  formulation,  the  objective  function  (6.13)  and  constraints  (6.14), 
(6.15)  are  obtained  from  the  corresponding  objective  function  (5.16)  and  con¬ 
straints  (5.17),  (5.29)  by  substituting  the  expressions  for  %kp-  and  -variables 
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from  (6.1)  and  (6.2),  respectively.  Notice,  that  since  constraint  (5.30)  from  the 
basic  MIP  formulation  is  omitted  in  the  aggregated  formulation,  the  value  of  p 
may  exceed  one.  If  this  is  the  case,  then  the  network  is  capable  of  carrying  an 
amount  of  traffic  increased  by  the  factor  of  p.  Thus  p  can  be  interpreted  as  the 
measure  of  network  over-dimensioning. 

Notice  that  if  any  constraint  from  the  group  (6.14)  is  satisfied  as  a  strict 
inequality,  then  the  corresponding  variable  can  be  decreased  which  would 
cause  the  decrease  of  flow  between  the  nodes  i,  j  and,  as  a  consequence,  non¬ 
decrease  of  the  objective  function  (6.13).  This  observation  can  be  formulated  as 
the  following: 

Proposition  5  There  is  always  an  optimal  solution  to  (AGGH)  in  which  con¬ 
straints  (6.14)  are  satisfied  as  equations. 

Consequently,  in  order  to  simplify  the  formulation  of  the  problem,  we  can  elim¬ 
inate  '(/-variables  by  a  substitution: 

Vij  :=  cl^jp  for  all  i  /  j  E  A”,  h  ETL 

where  dfj  =  Ylkeich  rk  is  the  class  h  demand  between  nodes  i  and  j. 

After  the  substitution  the  problem  is  reformulated  equivalently  as: 

AGGR: 


Maximize 
subject  to 


(6.13) 


E 


rh  -  V  Xh 

JIP  /  J  IP 

peouT(j)  pem(j) 


-dijp  if  j  /  i 


E  dimP  if  3  = 

(6.3-6.12)  and  (6.16-6.18) 


V  i,j  E  J\f ,  h  E  TL  (6.20) 


%ip,  P  >  0 


for  all  i  E  fif,  p  EV,  h  ETl 


(6.21) 
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6.1.1  Computational  Model  for  the  Fault-Tolerant  VP 
Rerouting 

In  this  section  we  continue  our  work  on  defining  a  fault-tolerant  model  which 
was  started  in  section  5.2.5. 

One  more  comment  is  to  be  made  before  we  add  fault-tolerant  constraints  to 
the  problem  formulation  (AGGR).  According  to  our  linear  approximation  for 
the  capacity  allocation  rule,  when  multiclass  traffic  is  switched  from  one  VP  to 
another,  the  fixed  charge  is  equal  to  the  highest  fixed  charge  of  all  classes.  In 
other  words,  it  is  the  maximal  fixed  charge  of  the  two  VPs. 

We  first  introduce  necessary  notation: 

r(p)  ,  backup  for  the  VP  p.  (Since  each  VP  in  a  pair  is  backup  for  the  other, 
r(r(p ))  =  p) 

Qe  ,  amount  of  the  capacity  reserved  on  link  e  to  be  used  for  traffic  restoration 
in  case  of  any  link  failure; 

Zp  ,  increase  in  fixed  charge  on  the  backup  VP  r(p)  when  the  VP  p  fails. 

It  is  easy  to  see  that 

Zp  =  max{  Y,  ARh  (4  ~  4(p ■) )»  °)  (6-22) 

hen 

where  we  define  ARh  =  Rh  —  Rh-i  for  all  h  G  R  —  {0}  and  A R0  =  0. 

Finally,  the  aggregated  model  with  fault  tolerant  constraints  is  formulated 
as  follows: 
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AGGR-FT: 


Maximize  (6.13) 

subject  to 

(6.20) 


E  (EE 

pev-.eep  ieM  hen 

xip  E  ^-RhZp  )  +  Qe  V  Be 

hen 

for  all  e  E  A 

(6.23) 

Qe>  E  EE4U 

p:e'  ep,  eer(p)  ieAf  hen 

for  all  e  /  e'  G  4 

(6.24) 

zr>Y.mj4-P,„) 

hen 

for  all  p  E  V 

(6.25) 

Zp>  0 

for  all  p  E  V 

(6.26) 

(6.11) 

(6.17-6.21) 

In  this  formulation,  we  reserve  amount  Qe  of  link  e  capacity  (constraint  (6.23)) 
sufficient  to  restore  failed  VPs  whose  backups  traverse  link  e  (constraint  (6.24)). 
The  fixed  charge  increment  Zp  for  VP  p  backup  is  determined  by  constraints 
(6.25-6.26)  which  are  equivalent  to  expression  (6.22). 

6.2  Valid  Inequalities  for  the  Fixed-Charge  Lin¬ 
ear  Approximation 

In  this  section  we  present  a  polyhedral  approach  to  solving  the  VP  layout  prob¬ 
lem  in  which  the  capacity  allocation  rule  for  all  types  of  traffic  is  assumed  to  be 
approximated  by  a  fixed-charge  model.  As  was  pointed  out  earlier,  this  is  a  case 
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when  all  traffic  is  video  or  saturated  data.  Then  the  constraint  set  is  as  follows: 


(6.20) 

Y  (Y  Y  xiP  +  Y  ARhz^)  <  for  a11  e  e  A  (6.27) 

pev-.eep  ieJ\f  hen  hen 

(6.11) 

(6.17-6.21) 

6.2.1  Band  Inequalities 

Our  approach  to  formulating  valid  inequalities  in  this  section  follows  methodol¬ 
ogy  proposed  in  [42,  11], 

We  start  by  assuming  that  lower  bound  on  p  is  known  a  priori.  That  is 
P>  Po 

This  constraint  can  be  either  explicitly  present  in  the  problem  formulation  or  a 
lower  bound  might  be  known  in  advance. 

Let  us  introduce  some  notation: 

P  =  (Pe  >  0,  e  G  £),  vector  of  arbitrarily  assigned  link  lengths; 

7t,L  ,  shortest  path  between  nodes  i  and  j  with  respect  to  link  lengths  /7: 

Vp'  =  {p  :  E eep  Pe  >  0},  set  of  positive  length  VPs  with  respect  to  link  lengths 

jj. 

In  addition,  let  Be  be  a  link  e  capacity  and  D ^  be  a  flow  demand  between  nodes 
i  and  j.  Then  a  multicommodity  flow  satisfying  the  demand  can  be  routed 
through  the  network  if  and  only  if 

PeBe  >  ::i.i 1  for  a11  P  -  0  (6-28) 

ee£  i,je  Jf 
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This  characterization  of  the  capacity /demand  for  a  network  carrying  multicom¬ 
modity  flow  was  first  stated  in  [19,  34]  and  can  be  proved  using  linear  program¬ 
ming  duality. 

Notice  that  in  the  trivial  case  when  fj,e  =  1  for  links  e  forming  a  (multi-)cut 
and  zero  otherwise,  condition  (6.28)  becomes  equivalent  to  minimum  (multi-)cut 
necessary  conditions  on  multicommodity  flow  feasibility. 

In  our  case  after  substituting 

Be  =  Be-  £  Y,  AR>*4 
p-  eep  hen 

A. j  =  Z  My 

hen 

into  (6.28)  and  regrouping  the  parts,  we  obtain 

E  He)ARhz%  <J2Be-  My)  for  a11  P  >  0  (6-29) 

pfz-pp  hen  eep  ees  i,jeN  hen 

In  Figure  6.2  we  picture  our  approach  to  deriving  a  valid  inequality  based  on 
(6.29)  for  a  fixed  /X  Each  column  of  the  grid  corresponds  to  a  VP,  p  G  and 


Figure  6.2:  Band  Inequality 
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^  a, 


is  Y,eep  /G  wide;  each  row  corresponds  to  traffic  class  h  G  %  and  is  A/?/,  high. 
The  shaded  area  is  equal  to  Yeep  tleRhZpp  and  chosen  to  exceed  the  right-hand 
side  of  (6.29)  (denoted  further  as  RHS(/2)).  This  means  that  in  order  to  satisfy 
flow  feasibility  conditions  (6.28),  at  least  one  of  the  variables  z^p  has  to  be  zero. 
We  call  this  valid  inequality  “band  inequality”  following  [11]  and  formulate  it  as 
a  proposition: 

Proposition  6  If  for  any  p  G  VJ  hp  is  chosen  such  that 

£  £(£/<«)«»,  >£b,-  £  4(£M‘) 

pern1  hen  eep  ees  ije Af  hen 

then 

E  4P  <  1^1  - 1  (6-3°) 

pevP 

is  a  valid  inequality  for  (AGGH). 

Separation  Problem 

We  formulate  a  separation  problem  for  generating  violated  constraints  (6.30)  as 
a  multi-constraint  knapsack  problem  finding  the  most  violated  constraint  if  there 
is  any.  Though  the  knapsack  problem  is  NP-hard,  it  allows  pseudo-polynomial 
solution  algorithms.  Also,  the  rounding  heuristics  based  on  LP-relaxations  are 
normally  quite  effective. 

Variables: 

|  1  if  h  =  hv 

<  for  all  p  G  V11  and  h  G  H 

0,  otherwise 
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SP: 


Maximize  E  PevpT.henzPfP  (6-31) 

Subject  to 

T,hen  tp  =  1  for  all  p  G  (6.32) 

ZPev(T,  eePVe)Rhthp>  RHS*  (6.33) 

e  {0, 1}  for  all  p  e  V ,L  and  h  e  %  (6.34) 


If  the  problem  (SP)  is  not  feasible  or  the  value  of  the  optimal  solution  does 
not  exceed  [P^\  —  1  (right-hand  side  of  (6.30)),  then  all  valid  inequalities  for 
the  given  values  of  p0  and  p  are  satisfied;  if  the  optimal  solution  value  exceeds 
\V^\  —  1,  a  new  valid  inequality  is  added  to  the  formulation.  The  problem  solution 
method  can  be  easily  modified  to  allow  for  generating  more  than  one  violated 
constraint. 


Choice  of  fi 

As  was  mentioned  earlier,  any  cut  induces  a  vector  p  such  that 


pe 


1  if  arc  e  belongs  to  the  cut 
0  otherwise 


(6.35) 


In  our  computational  experiments,  we  consider  a  set  of  cuts  C  =  {(jV  — 
{?'},  {«})|  i  G  A/"}.  For  each  cut  in  C,  vector  p  is  defined  by  (6.35)  and  used  to 
solve  the  separation  problem  (SP)  which  generates  violated  band  inequalities  in 
the  form  (6.30). 
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6.2.2  Connectivity  Inequalities 

Consider  a  cut  C  =  ( V,  W  =  J\f  —  V)  in  the  network  Q  =  (A/",  V)  where  V  C  Af . 
Suppose  the  total  demand  between  the  node  subsets  defining  the  cut  is  non-zero. 
That  is,  3 h  G  H  :  DyW  =  Y,ij-.  ievjew  ^ij  >  0-  Let  ho  be  the  highest  among 
such  h.  (Here,  we  omit  indices  ij  at  h0).  Then,  in  order  to  provide  non-zero 
network  throughput,  at  least  one  arc  from  the  cut  arc  set  should  be  open  for 
class  ho  traffic.  If  we  denote  this  arc  cut  set  by  Vvw  =  U ievjew^Vo  then  the 
above  statement  can  be  formulated  as  the  following  proposition. 

Proposition  7  If  for  a  cut  (V,  W)  there  exists  h0  such  that 

ho  =  max{/i  :  D 'yW  >  0},  (6.36) 

then 

£  %*”>!  (6.37) 

P&'Pvw 

is  a  valid  inequality  for  (AGGR). 

Separation  Problem 

Clearly,  there  are  in  general  exponentially  many  inequalities  (6.37)  for  a  given 
problem  instance.  However,  once  the  linear  programming  relaxation  of  (AGGR) 
is  solved,  the  violated  connectivity  constraints  (6.37)  can  be  found  in  polynomial 
time  using  the  following  approach. 

For  every  node  pair  (/,  j)  such  that  there  exists  h0  =  max{/i  :  dff  >  0}, 
solve  a  single  commodity  flow  problem  sending  a  unit  amount  of  flow  from  i  to 
j  over  arc  set  V  with  capacities  z}f° .  If  the  flow  problem  is  not  feasible,  add 
valid  inequalities  (6.37)  for  the  min  cut  (VW)  to  the  constraint  set  of  problem 

(AGGR). 
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6.2.3  Capacity  Inequalities 


Capacity  inequalities  can  be  considered  as  strengthened  connectivity  inequalities 
introduced  in  Section  6.2.2  in  a  special  case  of  equal  capacity  links.  This  type  of 
valid  inequality  was  considered  in  the  literature  [7,  31]  as  applied  to  a  network 
design  problem. 

Similar  to  the  way  it  was  introduced  in  section  6.2.1,  we  assume  that  a  lower 
bound  on  variable  p,  denoted  by  po,  is  known. 

Consider  a  cut  (V,  W).  If  all  physical  links  associated  with  the  cut  are  of  the 
same  capacity  B  and  a  fixed  charge  associated  with  class  h  is  Rh,  then  in  order 
to  accommodate  demand  po -Dyyy  across  the  cut,  at  least  \poDyW/(B  —  Rh)] 
links/VPs  should  be  open.  This  can  be  formulated  as  the  following  proposition. 


Proposition  8  If  all  links  associated  with  cut  (V,  W)  have  the  same  capacity 

B, 


E  4>\PoD™ 


B  —  Rt 


pE'Pyw  ^  ^ 

is  a  valid  inequality  for  fAGGRj. 


(6.38) 


Separation  Problem 


Using  approach  outlined  in  [31],  we  propose  the  following  heuristic  to  detect 
some  of  the  violated  capacity  inequalities  (6.38). 

STEP  1.  Test  valid  inequality  (6.38)  for  all  h  G  R  and  all  V  =  V7;  =  { i } . 
i  G  A f .  If  a  violated  inequality  is  detected,  add  it  to  the  problem  and  stop, 
otherwise  goto  STEP  2. 

STEP  2.  Grow  node  sets  Vp 


V*  :=  Vi  [J{.f  =  arg  min 


Po  D 


ViWi 


jeMhen  B  —  Rh 
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where  W,  =  AT  -  V,.  Goto  STEP  3. 

STEP  3.  If  there  are  any  violated  inequalities  found  in  STEP  2,  add  them 
to  the  problem  and  stop,  otherwise  repeat  STEP  2  while  \Vt  I  <  l-A/1- 

The  idea  here  is  to  apply  a  kind  of  “greedy”  approach  in  order  to  find  as 
many  violated  inequalities  as  possible. 

6.2.4  Bundle  Inequalities 

As  it  was  discussed  earlier,  the  problem  under  the  scope  is  essentially  a  bi¬ 
criteria  optimization  problem.  In  the  marginal  case,  we  wish  to  maximize  the 
throughput  when  the  flow  through  intermediate  nodes  is  equal  to  zero.  That 
means  that  at  least  one  direct  end-to-end  VP  is  set  up  between  each  pair  of 
nodes  with  non-zero  traffic  demand. 

Proposition  9  If  flow  through  intermediate  nodes  is  constrained  to  be  zero,  that 
is, 

E  E  E  4  = 0 

ieM per,  pgouT(i)  hen 

and  there  exists  ho  =  max}//  :  dh  >  0} 
then 

£  4°  >  1  (6-39) 

peVij 

Notice  that  in  case  of  a  single  VP  candidate  between  each  node  pair,  constraint 
(6.39)  uniquely  defines  the  solution. 
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6.3  Algorithms  and  Heuristics 

Linear  programming  relaxation  of  the  problem  (AGGR-FT)  together  with  var¬ 
ious  valid  inequalities  presented  in  sections  6. 1-6. 2  allows  us  to  obtain  an  upper 
bound  on  the  solution  value.  However,  it  does  not  provide  us  with  a  feasible 
solution.  Also,  the  number  of  VP  candidate  paths  between  node  pair  (i,j)  con¬ 
tained  in  a  VP  candidate  set  Vij  can  be  exponentially  large.  In  this  section  we 
discuss  various  aspects  of  obtaining  feasible  solutions  using  a  limited  number  of 
VP  candidate  routes  and  outline  a  general  procedure  for  solving  the  problem. 


6.3.1  Feasible  Solution  Heuristic 


Problem  formulation  (AGGR-FT)  contains  zero-one  ^-variables.  Since  the 
problem  size  is  quite  large,  it  does  not  appear  practically  possible  to  find  an 
exact  optimal  solution.  Instead,  we  propose  a  heuristic  method  to  obtain  a 
feasible  solution  from  a  solution  to  the  linear  programming  (LP)  relaxation  with 
little  additional  computations. 

Consider  solution  values  z ^  of  the  LP  relaxation  variables  z%.  Suppose,  we 
are  given  a  threshold  X  such  that  0  <  T  <  1.  If  variables  are  fixed  to  their 
rounded  values  : 


0  if  V  <  T 
-h  _  J  p  ^ 

Zp 

1  otherwise 


(6.40) 


then  it  is  easy  to  see  that  the  rounded  values  satisfy  order  constraints  (6.11). 
Now,  if  problem  (AGGR-FT)  is  solved  again  as  an  LP  with  its  ^-variables  fixed 
to  the  z^  values  and  has  a  feasible  solution,  this  new  solution  will  be  IP-feasible. 
The  remaining  question  is  how  to  choose  a  rounding  threshold  T.  If  the 
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chosen  threshold  is  too  large,  too  few  ^-variables  will  be  fixed  to  one  and  the 
network  may  become  disconnected.  On  the  other  hand,  if  the  chosen  threshold 
is  too  small,  too  many  ^-variables  will  be  fixed  to  one  and  the  throughput  will  be 
affected.  Our  computational  experiments  showed  that  the  solution  value  is  not 
very  sensitive  to  the  particular  choice  of  the  threshold.  In  order  to  find  the  best 
solution  by  rounding,  we  try  several  thresholds  starting  with  initial  threshold 
value  Ti  and  sequentially  decreasing  the  threshold  value  by  the  constant  factor. 


6.3.2  Selection  of  Initial  VP  Route  Set 


As  was  stated  in  Proposition  1,  the  optimal  solution  to  VP  layout  problem  with 
objective  function  weight  coefficient  r)  =  0  (pure  throughput  maximization)  and 
without  fault-tolerant  constraints,  does  not  involve  VPs  longer  than  one  hop. 
Thus,  the  optimal  solution  can  be  obtained  without  creating  VP  candidate  route 
sets  Vij.  We  use  this  property  to  obtain  an  initial  solution  which  is  later  used 
to  generate  an  initial  set  of  VP  routes.  Note  that  for  every  VP  we  must  find  a 
link-disjoint  backup  VP. 

The  heuristic  can  be  outlined  as  follows.  Consider  a  solution  to  the  VP 
layout  problem  with  pure  throughput  maximization  and  with  the  set  Vtl  defined 
as  follows: 


V,  =  < 


( '• ./ )  c  A 


(6.41) 


0  otherwise 

This  solution  provides  flow  values  from  any  node  i  obtained  by  summing  the 
variables  x\e  over  all  traffic  classes.  For  any  node  pair  (i,j)  consider  a  highest 
capacity  flow  path  ptj .  This  path  is  defined  as  a  path  maximizing  the  lowest 
flow  value  on  its  arcs  and  can  be  found  by  a  simple  modification  of  the  shortest 
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path  algorithm.  In  our  implementation  we  use  the  Bellman-Ford  algorithm  since 
it  guarantees  that  the  resulting  path  will  have  the  lowest  number  of  hops.  The 
second  link-disjoint  backup  path  is  obtained  by  deleting  the  links  on  the  first 
path  and  reapplying  the  method. 

6.3.3  Outline  of  the  Solution  Process 

In  this  section  we  give  a  general  outline  of  our  algorithm  by  providing  an  ap¬ 
proximate  solution  to  (AGGR-FT). 

1.  Solve  problem  (AGGR)  with  rj  =  0 

2.  Use  VP  routing  heuristic  to  obtain  a  set  of  link-disjoint  path  pairs 

3.  for  rj  =  770,?7i, . . .  ,r)N  do 

4.  Solve  an  LP  relaxation  of  (AGGR-FT)  with  valid  inequalities 

5.  for  T  =  Ti,  T2, . . . , Tm  do 

6.  Round  ^-variables  values  using  threshold  T 

7.  Solve  problem  (AGGR-FT)  with  ^-variables  fixed 

8.  Choose  the  best  solution  to  (AGGR-FT)  with  fixed  variables 

9.  Choose  a  subset  of  non-dominated  solutions  to  (AGGR-FT) 

In  this  algorithm,  values  of  the  T]  coefficients  were  chosen  empirically  such  that 
rj o  =  0,  7]n  was  large  enough  to  provide  end-to-end  VP  layout  (Proposition  2). 
As  was  mentioned  earlier,  different  threshold  values  Ti,T2,  . . . ,  TM  are  obtained 
by  sequentially  decreasing  its  initial  value  by  a  constant  factor.  Our  particular 
choice  was  the  initial  value  of  1/2  and  the  decrement  factor  of  2.  Thus,  Tm  =  2~m. 
Details  of  the  LP  computations  will  be  described  in  the  next  section. 
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6.4  Computational  Experiments 


In  this  section  we  present  results  of  the  computational  experiments  and  software 
implementation  of  our  approach  to  solving  problem  (AGGR-FT).  Our  compu¬ 
tational  experiments  were  designed  to  demonstrate  telecommunication  as  well 
as  combinatorial  optimization  aspects  of  the  problem.  We  also  present  instances 
of  the  test  problems. 

6.4.1  Implementation 

The  algorithm  was  implemented  on  SUN  SPARC  10  workstation,  coded  in  C 
and  used  CPLEX  3.0  LP  solver  callable  library.  Below  we  give  some  detail  of 
the  implementation. 

Since  our  problem  is  a  generalization  of  a  multicommodity  flow  problem, 
we  developed  an  approach  which  would  take  advantage  of  the  embedded  net¬ 
work  structure.  An  initial  basis  was  obtained  by  solving  a  problem  with  single- 
commodity  network  flow  constraints  (i.e.  a  subproblem  of  the  original  problem 
obtained  by  omitting  coupling  constraints)  by  using  CPLEX’s  built-in  network 
optimizer.  The  resulting  basis  was  not  feasible  but  nevertheless  provided  a  good 
starting  point  for  solving  the  LP  relaxation.  At  this  stage  we  used  the  primal 
simplex-method  providing,  as  was  empirically  shown,  a  faster  way  to  solve  the 
problem  than  dual  simplex.  The  resulting  optimal  solution  was  used  to  obtain 
rounded  2-variable  values  and  its  optimal  basis  was  used  as  initial  dual-feasible 
basis  to  solve  the  problem  with  fixed  2-variables.  The  latter  problem  was  solved 
several  times  for  various  rounding  thresholds.  Since  its  initial  basis  was  dual- 
feasible  we  applied  the  dual  simplex-method.  To  initiate  a  solution  process  for 
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a  linear  programming  relaxation  with  a  new  ij  objective  function  coefficient,  the 
optimal  basis  from  the  previous  iteration  was  used.  Since  this  basis  was  primal- 
feasible,  the  primal  simplex  was  used. 

6.4.2  Test  Data 

Our  computational  experiments  were  conducted  for  two  network  configurations 
presented  in  Figures  6.3  and  6.4.  The  first  is  an  example  of  a  grid  network  with  8 
nodes  and  10  undirected  links,  the  other  represents  a  geographically  distributed 
network  with  12  nodes  and  20  undirected  links.  All  links  in  these  two  networks 
have  the  same  capacity  of  250  Mb/s.  Traffic  demand  between  each  node  pair 
was  generated  randomly  using  the  same  uniform  distribution. 


0 - 0 - 0 - 0 


6 - 6 - 6 - 6 


Figure  6.3:  Grid  Network  Configuration 


6.4.3  Telecommunication  Aspects 

We  present  solutions  to  the  bi-criteria  optimization  VP  layout  problem  in  the 
form  of  a  curve  connecting  approximate  Pareto-optimal  solution  points.  This 
curve  demonstrates  a  tradeoff  between  relative  throughput  and  VP-to-VP  switch¬ 
ing  (intermediate  flow).  We  present  these  results  in  Figures  6.5  and  6.6  for  the 
“grid”  and  “distributed”  networks,  respectively.  The  curves  are  shown  by  solid 
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Figure  6.4:  Geographically  Distributed  Network  Configuration 

lines.  As  can  be  seen,  the  throughput  increase  can  be  achieved  only  at  the  price 
of  increased  flow  through  intermediate  nodes.  In  the  accompanying  Tables  6.1, 
we  provide  additional  statistics  for  the  VP  layout  such  as  the  number  of  VPs, 
average  hop  length  and  allocated  capacity.  In  accordance  with  the  theory,  the 
VPs  tend  to  become  “slimmer”,  “longer”  and  more  numerous  as  we  increase  the 
weight  coefficient  rj. 

VP  Layout  Adaptability 

As  was  mentioned  in  the  introduction  to  Chapter  5,  we  define  adaptability  as 
the  ability  of  a  particular  telecommunication  network  (and  its  VP  layout)  to 
function  well  when  the  traffic  demand  varies  significantly  from  the  projected 
traffic  (under  which  its  design  was  based). 

Our  conjecture  was  that  layout  using  “fat”  and  “short”  VPs  should  have 
better  adaptability  than  a  layout  using  narrower  and  longer  VPs. 
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Grid”  Network 


V 

Number 

of  VPs 

Average  Number 

of  Hops  on  VP 

Average  VP 

Capacity,  Mb/s 

Relative 

Gap,  % 

0.0000 

47 

2.15 

39.784 

2.17 

0.0001 

85 

2.48 

16.093 

7.84 

0.0002 

88 

2.55 

15.941 

8.28 

0.0004 

93 

2.59 

12.862 

10.53 

0.0006 

97 

2.59 

11.545 

1.86 

“Distributed”  Network 

Number 

Average  Number 

Average  VP 

Relative 

V 

of  VPs 

of  Hops  on  VP 

Capacity,  Mb/s 

Gap,  % 

0.000 

153 

1.95 

29.608 

3.18 

0.00005 

164 

1.97 

26.449 

4.38 

0.001 

190 

2.16 

18.259 

7.11 

0.002 

188 

2.14 

17.418 

8.35 

0.003 

200 

2.18 

16.181 

9.74 

0.004 

208 

2.21 

14.941 

1.90 

Table  6.1:  VP  layout  Statistics 
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50  100  150  200  250  300  350  400  450  500  550 

Total  Flow  through  Intermediate  Nodes 

Figure  6.5:  Results  for  “Grid”  Network  -  Tradeoff  Curve  (solid  line)  and  Band¬ 
width  Blocking  Probabilities  for  Various  VP  Layouts  (dashed  line). 

We  used  call  blocking  probability  as  a  measure  of  network  adaptability.  Since 
in  our  model  different  calls  have  different  bandwidth,  we  are  using  weighted  call 
blocking  probability,  or  bandwidth  blocking  probability  defined  as  the  ratio  of 
the  total  rejected  bandwidth  to  the  total  bandwidth  requested. 

In  order  to  confirm  our  conjecture,  we  conducted  a  series  of  computational 
experiments  simulating  call  requests  -  setup/blocking  -  takedown  using  a  pre¬ 
computed  VP  layout.  We  considered  several  sets  of  call  requests  obtained  by 
increasing  projected  traffic  demand  in  the  VP  layout  by  1,  10  and  30  %.  After 
that  the  requested  connections  were  routed/rejected  using  methods  described 
in  Chapter  3.  The  networks  used  consisted  of  VPs  laid  out  according  to  the 


119 


Figure  6.6:  Results  for  “Distributed”  Network  -  Tradeoff  Curve  (solid  line)  and 
Bandwidth  Blocking  Probabilities  for  Various  VP  Layouts  (dashed  line). 

solution  of  problem  (AGGR)  with  varying  bi-criteria  objective  function  weight 
coefficients  for  the  throughput  -  intermediate  flow  criteria. 

The  results  are  shown  in  Figures  6.5  and  6.6  for  the  “grid”  and  “distributed” 
networks,  respectively.  In  these  figures,  dotted  lines  connect  points  representing 
different  bandwidth  blocking  probabilities  (multiplied  by  5  for  better  represen¬ 
tation  purposes)  for  different  VP  layouts.  The  three  lines  correspond  to  three 
traffic  demand  increments  mentioned  above. 

It  can  be  seen  from  the  figures  that  “fat”  and  “short”  VPs  have  better  adapt¬ 
ability  than  “slim”  and  “long”  VPs. 
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6.4.4  Combinatorial  Optimization  Aspects 

As  was  mentioned  earlier,  it  practically  impossible  to  find  an  optimal  solution  to 
the  problem.  Instead,  we  produce  an  approximate  solution  and  estimate  quality 
of  the  approximation. 

From  the  previous  discussion  of  the  solution  process,  it  can  be  concluded  that 
there  are  two  factors  which  may  cause  suboptimality  of  the  solution: 

1.  instead  of  generating  all  paths  between  each  node  pair  within  the  linear 
program,  we  obtained  a  single  pair  using  a  multicommodity  flow  heuristic; 

2.  instead  of  solving  IP  to  optimality,  we  chose  a  feasible  solution  using  a 
certain  rounding  heuristic. 

Since  it  is  difficult  to  give  a  quantitative  estimation  to  the  loss  of  accuracy 
caused  by  the  first  factor,  we  argue  that  potential  loss  of  accuracy  is  not  sub¬ 
stantially  high  by  comparing  layouts  for  one-  and  two-route  per  node  pair  VP 
candidate  sets.  The  comparison  results  obtained  for  the  solution  to  the  (AGGR) 
problem  are  shown  in  Figures  6.7  and  6.8  for  “grid”  and  “distributed”  networks, 
respectively. 

It  can  be  seen  that  the  curve  shapes  are  very  similar  and  adding  the  second 
path  to  a  VP  candidate  set  does  not  substantially  improve  the  solution. 

An  accuracy  loss  caused  by  the  second  factor  which  is  also  known  as  integral¬ 
ity  gap  can  be  estimated  by  comparing  solution  values  of  the  LP  relaxation  and 
feasible  solutions.  Quantitative  results  are  provided  in  Tables  6.1  (last  column) 
which  contains  the  ratio  (per  cent)  of  the  integrality  gap  (difference  in  value  of 
integer  and  LP  solutions)  and  the  value  of  the  linear  relaxation  solution. 

It  can  be  seen  from  the  Tables  6.1  that  the  integrality  gap  does  not  exceed 
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2.9 


Figure  6.7:  Pareto-optimal  Solution  Sets  for  One-  and  Two-  VP  Candidate  Route 
Approaches  (“Grid”  Network) 

10%  and  it  becomes  quite  low  for  the  last  solution  corresponding  to  end-to-end 
VP  layout  where  we  can  apply  rather  strong  “bundle”  valid  inequality  (6.39). 
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Total  Flow  through  Intermediate  Nodes 


Figure  6.8:  Pareto-optimal  Solution  Sets  for  One-  and  Two-  VP  Candidate  Route 
Approaches  (“Distributed”  Network) 

6.5  Conclusions 

In  this  chapter  we  have  presented  models  and  algorithms  which  address  the  fault- 
tolerant  VP  layout  problem.  Using  VPs  as  a  backup  mechanism  enables  ATM 
networks  to  respond  in  real  time  to  single-link  failures  by  switching  traffic  to  pre¬ 
computed  backup  routes.  The  amount  of  time  required  to  perform  such  rerouting 
is  negligibly  low  and  determined  by  the  speed  with  which  VPI  switching  table 
is  changed  in  the  beginning  node  of  each  affected  VP. 

Our  problem  formulation  involves  bi-criteria  optimization  for  throughput 
maximization  and  control  overhead/delay  minimization  objectives.  This  ap- 
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proach  results  in  output  consisting  of  a  solution  set  also  called  a  tradeoff  curve 
which  provides  a  network  manager  with  a  useful  decision  making  tool. 

In  order  to  carry  out  the  implementation  of  our  methods,  specific  protocols 
must  be  designed  enabling  VP  switching  as  soon  as  link  failures  are  detected. 
Furthermore,  refinements  to  our  algorithms  should  be  made  to  insure  efficient 
interaction  with  network  operation  and  management. 
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Chapter  7 


Conclusions 


We,  now,  summarize  the  contributions  of  the  dissertation,  and  describe  poten¬ 
tial  directions  for  the  future  research.  Since  our  work  was  essentially  interdis¬ 
ciplinary,  involving  telecommunication  as  well  as  mathematical  programming 
issues,  we  consider  our  contributions  to  both  of  these  areas. 

Modern  telecommunication  networks  are  designed  to  carry  multiclass  traffic 
and  often  have  a  geographically  distributed,  complex  structure.  In  this  situa¬ 
tion  it  becomes  very  important  to  develop  reliable  mechanisms  working  in  both 
reactive  and  proactive  modes  for  responding  to  network  failures  without  inter¬ 
ruption  of  real-time  network  services.  The  primary  contribution  of  our  work 
in  the  area  of  fault  management  of  telecommunication  networks  is  the  devel¬ 
opment  of  tools  capable  of  restoring  traffic  under  a  wide  variety  of  conditions 
and  network  architectures.  In  chapters  4  and  5  we  considered  two  approaches  to 
rerouting:  reactive  and  proactive,  performed  on  VC  and  VP  levels,  respectively. 
The  first  approach  was  designed  to  respond  to  multiple  link  failures  and/or  for 
use  in  telecommunication  networks  lacking  two  layer  VC/VP  architecture.  The 
other  approach  which  is  less  universal  but  more  effective  was  designed  to  respond 
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to  single-link  failures  -  the  most  common  failure  type  in  ATM  networks.  It  is 
proactive  and  explicitly  predefines  a  set  of  VP  switching  operations  to  fix  any 
potential  link  failure. 

Certain  multiclass  telecommunication  networks  allow  for  preemption  of  lower 
priority  services  when  higher  priority  services  are  affected  by  a  network  failure 
and  need  to  be  rerouted.  This  feature  was  incorporated  in  to  our  model  to 
compute  VC-level  rerouting.  Another  feature  of  modern  telecommunication  net¬ 
works  -  their  ability  to  perform  distributed  computations  -  was  also  used  in  our 
model. 

Our  approach  to  VP  layout  can  be  considered  as  an  advanced  approach  to 
resource  allocation.  The  novel  features  here  are  enhancement  of  existing  VP 
backup  schemes  allowing  for  non-zero  traffic  over  the  backup  VP  and  bi-criteria 
optimization  involving  throughput  maximization  and  control  overhead/delay 
minimization  objectives.  This  approach  results  in  output  consisting  of  a  so¬ 
lution  set  also  called  a  tradeoff  curve  which  provides  a  network  manager  with  a 
useful  decision  making  tool.  It  requires  a  minimal  set  of  on-line  operations  and 
because  of  that  can  be  implemented  over  broad  spectrum  of  network  protocols. 

Development  of  these  two  approaches  can  be  considered  as  a  first  step  towards 
an  integrated  fault  management  system.  The  future  research  directions  would 
involve  refinements  of  the  proposed  models  and  design  of  the  protocols  which 
could  take  advantage  of  the  suggested  fault  recovery  techniques. 

Our  contribution  in  the  area  of  mathematical  programming  can  be  high¬ 
lighted  as  follows: 

•  Generalization  of  the  multicommodity  flow  problem  with  zero/one  vari¬ 
ables  by  including  preemption  variables  (chapter  4). 
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•  Design  of  new  solution  methods  oriented  toward  obtaining  a  high-quality 
approximate  solution  in  real  time. 

•  Development  of  a  two-phase  approach  taking  advantage  of  distributed  com¬ 
putational  resources. 

•  Piece-wise  linear  approximation  for  capacity  allocation  rule  in  case  of  sta¬ 
tistical  multiplexing. 

•  Integer  programming  formulation  for  VP  layout  problem. 

•  Design  of  solution  methods  for  VP  layout  problem  based  on  valid  inequal¬ 
ities  and  rounding  schemes. 

Finally,  we  list  the  most  relevant  future  research  directions.  The  first  two 
deal  with  telecommunication  aspects  of  the  VP  layout  problem;  the  last  two  are 
related  to  mathematical  programming  issues. 

1.  VP  layout  for  point-to- multipoint  and  multipoint-to- multipoint  connec¬ 
tions.  This  type  of  connection  is  expected  to  be  quite  common  since  it 
supports  video  conferencing  and  data  broadcasting  through  the  network. 
The  challenging  issue  here  is  large  amount  of  required  bandwidth  and  high 
burst  iness. 

2.  VP  layout  capable  of  supporting  rare  high-bandwidth  connections.  This 
connection  type  can  add  a  lot  of  variance  to  the  traffic  pattern  and  require 
high  adaptability  of  the  VP  layout. 

3.  Dynamic  VP  route  generation.  Since  our  VP  route  set  is  limited  to  only 
one  pair  of  link-disjoint  routes  between  any  node  pair,  there  is  a  need 
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to  develop  techniques  capable  of  dynamically  generating  additional  link- 
disjoint  VP  route  couples  during  the  solution  process.  Conceptually,  this 
process  should  be  similar  to  column  generation. 

4.  Narrowing  the  integrality  gap  for  the  approximate  solutions.  This  is  a 
classical  problem  faced  by  everyone  designing  approach  to  solving  an  IP. 
Since  the  formulated  problem  is  quite  large  and  has  a  complex  structure,  it 
does  not  seem  promising  to  apply  any  sort  of  branching  technique.  Instead, 
this  problem  should  be  solved  by  improving  the  quality  of  the  approximate 
solution  and  tightening  the  LP  relaxation  bounds. 
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